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Abstract 

Micro  ai  r  V  chicles  (  MAVs)  ar  e  i  ntended  f  or  f  uture  i  ntelligence,  s  urveillance,  an  d 
reconnaissance  use.  T  o  adequately  fulfill  a  cl  andestine  capacity,  MAVs  must  operate  in  close 
proximity  to  their  intended  target  without  eliciting  counter-observation.  This  obj  ective,  along 
with  D  area’s  constraint  of  a  sub-15  centimeter  span,  requires  future  MAVs  to  mimic  insect 
appearance  an  d  f  light  characteristics.  T  his  t  hesis  d  escribes  a  n  ex  perimental  m  ethod  f  or 
conducting  a  s  tructural  a  nalysis  of  a  Manduca  Sexta  (hawkmoth)  f  orewing.  G  eometry  i  s 
captured  via  computed  tomography  (CT),  and  frequency  data  is  collected  using  laser  vibrometry 
in  a  ir  and  va  cuum.  A  finite  e  lement  (EE)  mo  del  is  c  onstructed  using  qua  dratic  be  ams  a  nd 
general-purpose  shell  elements,  and  an  eigenanalysis  is  conducted.  A  preliminary  verification  of 
the  EE  model  is  carried  out  to  ensure  the  Manduca  Sexta  forewing  is  adequately  characterized, 
providing  a  ba  sis  f  or  f  uture  f  luid-structural  i  nteraction  c  omputations.  Included  i  s  a  s  tudy 
regarding  the  aeroelastic  effects  on  flapping-wing  insect  flight,  and  an  analysis  of  the  structural 
dynamic  anomalies  of  conventional,  fiat,  semi-rigid  flapping  wings. 

Experimental  tests  revealed  the  first  three  modes  of  a  cl  amped  Manduca  Sexta  wing  in 
vacuum  are  86  Hz,  106  Hz,  and  155  Hz;  tests  in  air  reveal  a  frequency  shift  of  26. 5%  from 
vacuum,  indicating  a  possible  aeroelastic  contribution  to  frequency  response.  The  finite  element 
model  produced  first  three  modes  of  84.6  Hz,  106.1  Hz,  and  317.7  Hz,  indicating  that  the  model 
is  limited  to  the  second  wing  mode  and  lower  frequencies.  Possible  sources  of  error  include  poor 
geometric  modeling  due  to  low  CT  resolution,  inadequate  modeling  of  camber,  and  inaccurate 
estimation  of  material  properties. 
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A  Structural  Dynamic  Analysis 


OF  A  Manduca  Sexta  Forewing 


I.  Introduction 


1.1.  Background 

Over  the  course  of  last  two  decades,  the  demand  for  unmanned  intelligence,  surveillance 
and  reconnaissance  (ISR)  assets  has  grown  at  an  extraordinary  pace.  Physical  separation 
of  an  aircraft  and  its  human  operator  has  provided  unprecedented  loiter  times,  range,  low-cost 
compared  toe  quivalent  m  anned  s  ystems,  and  pi  lot  s  afety.  T  o  supplement  future  ISR 
capabilities,  many  have  proposed  “micro”  air  vehicles  (MAVs)  that  would  revolutionize  the  field 
of  remote  sensing  with  their  low  cost^,  extreme  maneuverability,  and  inconspicuous  operation. 
As  pa  rt  of  t  heir  pr  oposed  m  ission,  t  hese  ve  hides  w  ill  ha  ve  a  w  ingspan  of  1  ess  t  ban  15 
centimeters,  c  omplete  with  f  lapping  wings,  to  min  imize  c  ounter-observation  [27].  If 
successfully  m  anufactured,  t  he  M  AV  of  t  he  f  uture  w  ill  be  c  apable  o  f  a  m  ultitude  of  c  lose- 
quarters  reconnaissance  tasks,  ranging  from  battlefield  operations  to  safety  inspections  of  civilian 
structures.  Table  1  outlines  the  Defense  A  dvanced  R  esearch  P  rojects  A  gency’s  (DARPA) 
vision  for  the  MAV  [38]. 

A  hr  ief  1  ook  a  t  any  c  onventional  flight-line  r  eveals  t  he  s  tark  co  ntrast  b  etween  current 
aircraft  design  and  that  of  the  proposed  MAV.  Today,  all  conventional,  powered  aircraft  feature 
prominent  engines  with  unnatural-looking  wings.  M  oreover,  a  primary  goal  in  wing  design  is 
maintaining  s  tructural  i  ntegrity  b  y  reducing  th  e  ma  gnitude  o  f  s  tructural  o  scillations.  E  ven 

*  Current  specifications  are  for  MAVs  to  cost  less  than  approximately  $1,500  (2009  USD)  [27]. 

^  Minimizing  counter-observation  is  not  only  visual.  D  etecting  an  appropriately  designed  MAV  out  of  a  crowd  of 
insects  with  a  radar  would  be  an  extraordinary  feat. 
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bounded  ha  rmonic  r  esponses  -  such  as  limit  c  ycle  o  scillations  -  are  avoided  toe  ircumvent 
material  fatigue  [23]. 


Table  1.  MAV  Design  Requirements  [38] 


Specification  Requirements  Details 


Size 

<  15.24  cm 

Weight 

~100g 

Range 

1  to  1 0  km 

Endurance 

60  min 

Altitude 

<  150  m 

Speed 

15  m/s 

Payload 

20  g 

Cost 

$1,500 

Maximum  Dimension 
Objective  GTOW 
Operational  range 
Loiter  time  on  station 
Operational  ceiling 
Maximum  flight  speed 
Mission  dependent 
Maximum  cost,  2009  USD 


Turning  to  nature  for  insight  into  future  MAV  designs  reveals  a  sharp  departure  from 
over  a  c  entury  of  d  esign  pr  actice:  not  onl  y  are  i  nsect  w  ing  os  cillations  p  ermitted,  th  ey  are 
abundant,  highly-efficient,  flight  mechanisms  [28].  Following  in  these  biological  footsteps,  the 
micro  a  ir  v  chicle  o  f  t  he  f  uture  w  ill  be  pu  rposely  d  esigned  w  ith  f  lapping  wings,  c  losely 
mimicking  insect  flight.  This  fundamental  shift  from  quasi-rigid  wings  -  be  they  fixed  or  rotary 
-  with  separate  power  plants,  to  a  flexible  airfoil  that  integrates  the  generation  of  lift  and  thrust, 
constitutes  a  dr  amatic  c  hange  i  n  de  sign  phi  losophy.  Through  r  esearch,  t  he  ubi  quitous  flying 
insects  serve  as  a  blueprint  from  which  structural  engineers  might  learn  important  design  criteria; 
among  these  insects  is  the  Manduca  Sexta,  or  hawkmoth  (see  Figure  1)  [20].  One  of  the  most 
notable  advantages  of  mimicking  Manduca  Sextae  is  that  they  lie  within  the  bounds  of  DARPA’s 
MAV  requirements  in  Table  1  [29]. 
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Figure  1:  An  adult,  kmaXe  Manduca  Sexta  (Hawkmoth)  [21]. 

Many  biologists  and  entomologists  have  studied  the  mechanics  of  Manduca  Sexta  flight. 
During  these  studies  and  ensuing  debates,  several  fundamental,  structural  questions  arose.  T  o 
date,  published  research  has  not  conclusively  addressed  such  important  issues  as  the  magnitude 
of  aer  oelastic  contributions,  the  s  tructural  lin  earity  of  the  w  ings,  t  he  n  atural  f  requencies  an  d 
modeshapes  of  the  wings,  or  the  level  of  sophistication  required  to  manufacture  a  wing  of  similar 
capabilities  and  characteristics.  With  these  questions  unanswered,  it  is  impossible  to  produce  a 
validated  structural  model  that  might  be  useful  in  future  numerical  studies. 

1.2.  Research  Objectives 

The  Manduca  Sexta  has  been  identified  as  an  ideal  research  candidate  to  assist  the  study 
of  insect  flight.  Still,  there  are  many  questions  that  remain  unanswered  by  previous  work.  Does 
the  Manduca  Sexta  beat  its  wings  near  its  natural  frequency?  Is  the  hawkmoth’s  forewing  a 
linear  or  nonlinear  structure?  Is  the  response  dominated  by  inertial  forces  or  aeroelastic  effects  - 
or  bot  h?  W  hat  1  evel  of  m  odel  s  ophistication  i  s  ne  cessary  t  o  adequately  represent  bi  ological 
wings?  A  sked  another  way,  what  simplifications  could  a  manufacturer  of  synthetic  Manduca 
Sexta  forewings  ma  ke  to  a  chieve  s  tructural  c  haracteristics  s  imilar  to  the  b  iological  w  ing? 
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Answering  t  hese  que  stions  r  equires  num  erical  a  nalyses  of  t  he  s  tructural  a  nd  fluid 
meehanics  t  hat  d  ominate  Manduca  Sexta  wings.  F  urthermore,  f  luid-struetural  in  teractions 
cannot  be  adequately  modeled  without  first  producing  a  validated  structural  model.  T  herefore, 
the  obj  ective  of  this  work  is  to  produce  a  finite  element  model  of  a  Manduca  Sexta  forewing 
grounded  i  n  experimental  vi  bration  t  esting.  In  s  o  doi  ng,  m  any  of  t  he  a  hove  que  stions  a  re 
addressed  and  answered,  laying  the  foundation  for  future  work. 

In  an  effort  to  elucidate  the  above  issues,  liberated  wings  from  Manduca  Sextae  are  used 
to  pursue  two  research  endeavors.  First,  the  natural  modes  (frequencies)  and  modeshapes  of  the 
hawkmoth  f orewings  are  i dentified  v  ia  1  aser  vibrometry.  Testing  is  a ccomplished  in  air  and 
vacuum,  to  observe  aeroelastic  effects.  Second,  a  finite  element  model  capable  of  representing 
the  observed  modal  behavior  is  developed  and  analyzed. 

1.3.  Motivation 

The  pace  of  unmanned  aerial  vehicle  (UAV)  development  and  operation  is  currently  at  an 
all-time  hi  storical  hi  gh,  w  ith  no  a  pparent  e  nd  i  n  s  ight.  World  a  ir  fo  rces  o  f  a  11  s  izes,  p  rivate 
corporations,  a  nd  e  ven  terrorist  gr  oups  c  urrently  operate  U  AVs  f  or  p  urposes  r  anging  f  rom 
weather  observation  to  munitions  delivery  [36].  As  the  number  of  UAVs  in  existence  grows,  it  is 
only  natural  that  the  level  of  specialization  of  these  systems  also  increases.  Due  to  the  widely- 
recognized  potential  afforded  by  MAVs,  there  are  a  multitude  of  nations  with  ongoing  research 
and  development  programs,  as  shown  in  Figure  2  [29].  M  icro  air  vehicle  development  may 
prove  to  be  a  critical  step  in  keeping  America's  technological  edge  on  the  battlefield. 

To  d  ate,  t  he  v  ast  m  ajority  o  f  M  AV  r  esearch  is  c  oncentrated  i  n  t  he  fields  of  f  luid 
mechanics  and  controls,  with  a  small  minority  -  about  10%  -  in  structures  [29].  As  a  result,  little 
work  has  concentrated  in  structural  dynamics  of  MAVs,  with  most  of  the  work  focusing  on  fluid- 


4 


structural  interaction.  Even  less  research  has  focused  on  the  structural  dynamics  of  insect  wings, 
and  their  applicability  to  bio-inspired  designs. 
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Figure  2:  International  MAV  research  locations.  Major  contributors  are  noted  by  name  [29]. 


1.4.  Flapping  as  a  Flight  Mechanism 

In  their  work  at  the  Georgia  Institute  of  Technology,  Michelson  and  Reece  [28]  produced 
several  compelling  arguments  for  flapping  wing  flight  as  the  preferred  mechanism  for  "stealthy" 
MAVs.  Their  rationale  for  selecting  flapping  wings  over  fixed-wing  or  rotary  flight  was  based 
upon  three  criteria:  necessary  flight  speed,  energy  efficiency,  and  stealthiness.  The  following 
three  paragraphs  are  a  summary  of  their  conclusions: 

Since  MAVs  will  likely  be  operated  in  urban  (including  indoor)  environments,  reaction¬ 
time  and  automated  obstacle  avoidance  are  critical  to  operational  success.  F  ixed  wing  designs 
require  i  ncreased  f  orward  ve  locity  to  sustain  flight,  min  imizing  th  is  p  recious  r  eaction-time. 
Moreover,  f  ixed-wing  designs  ne  ither  appear  nor  be  have  1  ike  s  imilarly-sized  hi  ological 
specimens,  increasing  the  probability  of  counter-observation.  T  hus,  fixed-wing  MAVs  are  ill- 
suited  for  their  intended  task. 
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Rotary-wing  aircraft  b  enefit  from  greater  maneuverability  at  low  flight-speeds,  and  are 
able  to  hover.  Several  commercial  helicopters  slightly  larger  than  the  MAV  class  already  exist 
with  high-resolution  cameras  (lOSOp  ).  These  devices  transmit  at  frequencies  ranging  from  80 
MHz  (FM)  to  10  GHz  (X-Band);  one  example  is  shown  in  Figure  3  [17]. 


Figure  3:  The  D  raganFlyer  X6  co  mmercial  U  AV  with  a  h  igh- 
resolution  (1080p)  camera,  has  an  85  cm  span  [3]. 


Flapping  wings  are  in  many  ways  a  compromise  between  fixed  and  rotary  wing  designs. 
As  e  videnced  i  n  n  ature,  f  lapping  w  ings  are  cap  able  o  f  1  ow-speed  f  light  and  hove  ring, 
maximizing  reaction  time  and  observation  c  apability.  Like  fixed-wing  aircraft,  flapping  w ing 
vehicles  have  chordwise  r  igidity,  allowing  t hem t  o  e  xpend  1  ess  e nergy  by  gliding.  Finally, 
flapping  wings  use  varying  wing  stroke  amplitudes,  producing  an  overall  average  audible  energy 
less  than  rotary  wings. 

Although  not  addressed  by  Michelson  and  Reece,  it  should  be  noted  that  several  hybrid 
designs  have  been  tested  to  address  the  issues  listed  above.  Figure  4  shows  one  such  example, 
the  Sikorsky  S  -72  X-Wing  h  ybrid  c  onfiguration  aircraft  [42].  A  s  ignificant  b  enefit  o  f  th  e 
stopped  rotor,  or  X-Wing,  configuration  is  its  ability  to  capitalize  on  t  he  maneuverability  of  a 

^  1080p  is  a  video  resolution  standard  with  1920  x  1080  pixels,  updated  sequentially  in  each  frame  [1]. 
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helicopter  and  efficiency  of  a  fixed  wing  aircraft.  Although  the  S-72  is  a  manned  aircraft  with  an 
18.9  meter  span,  its  capabilities  could  be  ported  to  the  MAV  class. 


Figure  4:  Sikorsky  S-72  X-Wing  hybrid,  manned  rotorcraft,  1986  [4]. 


Despite  its  design  advantages,  X-Wing  hybrids  are  still  comprised  of  unnatural-looking 
components,  and  emit  a  constant,  audible  whine.  Additionally,  hybrid  designs  such  as  the  S-72 
and  t  he  V  -22  O  sprey  s uffer  f rom  a  dded  w  eight  due  t  o  t  he  r  equired  s  ystem  c  omplexity  [  5],  a 
challenge  not  easily  overcome  by  MAVs. 


1.5.  Species  Selection 

Having  settled  on  a  flapping  wing  design,  with  the  intention  of  producing  an  MAV  easily 
mistaken  for  an  insect,  one  begins  the  tedious  task  of  selecting  an  insect  to  model.  Certainly 
there  are  a  plethora  of  species  that  lie  well  within  the  sub- 15  centimeter  wingspan  dictated  by 
DARPA,  e  ach  w  ith  i  ts  own  uni  que  set  of  s  tructural  a  nd  m  otive  be  haviors.  The  eh  allenge  of 
selection  is  ensuring  that  the  species  is  sufficiently  “typical”  -  else  any  knowledge  gained  during 
the  experimental  and  modeling  processes  might  be  too  narrowly  focused  for  broad  design  and 
application. 
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University  o  f  C  ambridge  r  esearchers  Willmott  a  nd  E  llington  reported  the  f  ollowing 
regarding  the  Manduca  Sexta  [40]; 

The  wingbeat  kinematics  contain  the  major  elements  of  which  have  been  consistently 
identified  in  other  insect  species,  but  without  any  of  the  marked  idiosyncrasies  which 
have  been  seen  in  certain  groups  ...The  stroke  amplitude  is  close  to  the  value  of  120°  that 
is  often  cited  as  typical  for  insects  \l>9\..The  hawkmoth,  therefore,  offers  a  reasonably 
generalized  wingbeat. 

Another  i  mportant  f  acet  of  the  Manduca  Sexta  is  t  he  c  onsistency  of  i  ts  w  ing  s  tructure  and 

deflection  from  specimen  to  specimen.  Willmott  and  Ellington  go  on  to  state  [40], 

The  Manduca  S  exta  wingbeat  proved  to  be  remarkably  consistent.  Significant 
intraspecific  variation  has  been  reported  for  many  insect  groups  including 
butterflies  [1]  and  flies  [8].  However,  the  variation  in  kinematics  between 
successive  wingbeats,  and  also  between  individuals,  was  small  in  the  current 
study. 

Einally,  the  Manduca  Sexta  is  consistent  in  the  frequency  of  its  wingbeat.  During  wind-tunnel 
testing,  W  illmott  a  nd  E  llington  f  ound  m  ultiple  s  pecimens  toe  xhibit  f  lapping  f  requencies 
between  24.8  and  26.5  Hz,  a  narrow  band  independent  of  flight  speed  [40]. 

Although  the  Manduca  Sexta  is  a  four-winged  synchronous  flapper,  Jantzen  and  Eisner 
experimentally  determined  that  th e  h awkmoth  c  an  f ly  w ithout  its  rear  wings  [24].  S  imilarly, 
Zhao  and  D  eng  [43]  used  a  s  ix-component  force  b  alance  t  o  the  c  ompute  1  ift  a  nd  dr  ag 
coefficients  0  f  Manr/wca  Sexta  wings.  U  sing  these  values,  Zhao  and  Deng  verified  that  the 
majority  of  the  in-flight  power  distribution  is  contained  in  the  forewings.  Based  on  the  above, 
sufficient  ev  idence  ex  ists  t  o  s  elect  t  he  Manduca  Sexta  forewing  as  a  prototypical  b  asis  f  or 
studying  flapping  wing  flight.  The  design  is  evidently  viable,  since  the  hawkmoth  is  capable  of 
flight  -  even  without  rear  wings. 
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1.6.  System  Identification 


After  s  uccessfully  id  entifying  th  e  Manduca  Sexta  as  an  i  nsect  w  ith  r  easonably  t  ypical 

wingbeat  characteristics,  Willmott  and  Ellington  posited  the  following  [40]; 

The  flight  musculature  must  function  at  the  natural  resonant  frequency  of  the 
mechanical  system  involving  the  thorax  and  wings.... Insects  appear,  therefore, 
to  fly  within  the  narrow  range  of  frequencies  which  maximizes  the  mechanical 
performance  of  their  flight  muscles. 

This  hypothesis  -  that  Manduca  Sextae  operate  near  or  at  a  resonant  frequency  of  their  wings 
where  displacement  is  maximized  -  is  based  upon  two  factors.  First,  the  Manduca  Sexta  has 
muscles  in  its  thorax,  but  none  in  the  wings.  This  configuration  produces  a  response  similar  to  a 
cantilever  b  earn  undergoing  ha  rmonic  base  m  otion.  S  econd,  multiple  s  tudies  d  emonstrate  a 
narrow  wingbeat  frequency  range. 

While  t  hese  lemmata  are  correct,  i  t  i  s  un  clear  whether  h  awkmoths  are  highly  e  fficient 
flyers  that  achieve  maximum  deflection  with  minimal  muscular  input  -  that  is,  operate  near  a 
resonant  frequency.  Without  validation,  this  inherent  assumption  seems  to  be  generally  accepted 
in  the  field  of  entomology  [10,  15]. 

From  a  s  tructural  d  ynamics  s  tandpoint,  one  m  ust  c  onsider  t  he  c  urious  a  bility  f  or  t  he 
Manduca  Sexta  to  fly  without  its  rear  wings,  if  indeed  the  forewings  operate  close  to  a  natural 
frequency.  If  the  f  orewing  is  lightly  damped,  t  he  ha  wkmoth  w  ould  be  a  ble  t  o  generate 
significant  displacements  -  thus  lift  -  by  simply  modulating  its  beat  frequency.  W  ith  a  sudden, 
dramatic  i  ncrease  i  n  en  ergy,  t  he  Manduca  Sexta  would  t  hen  b  e  m  ore  t  han  cap  able  o  f  ev  asive 
maneuvers,  eliminating  the  need  for  rear  wings  altogether.  Since  the  Manduca  Sexta  retains  its 
rear  wings  for  evasive  procedures  [24],  it  seems  unlikely  that  the  hawkmoth  “must”  operate  near 
resonance,  a  s  pos  ited  b  y  W  illmott  a  nd  E  llington  [  40] .  If  1  arge  d  amping  i  s  t  he  r  eason  f  or  r  ear 
wing  retention,  the  magnitude  of  the  wing  displacement  would  be  small,  even  near  resonance  - 
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also  casting  doubt  on  the  basis  of  the  elaim.  Therefore,  if  damping  is  present,  it  must  not  lie  near 
either  of  the  extrema.  A  counter  argument  in  support  of  Willmott  and  Ellington  is  the  possibility 
that  aeroelasticity  is  a  moderate  damping  meehanism.  If  so,  this  moderate  damping  would  both 
allow  operation  near  resonance  and  require  the  presence  of  rear  wings. 

1. 7.  The  Aeroelastic  Question 

After  e  ompleting  a  r  eview  of  ove  r  700  pa  pers  a  rchived  b  y  t  he  A  meriean  Institute  of 
Aeronauties  and  Astronauties  (AIAA),  Norris  [29]  was  unable  to  find  any  eonclusive  evidenee 
supporting  or  di sproving  a ero elasticity  as  a  dom inant  f actor  i n  t he  r esponse  of  flapping  wing 
flight.  Instead,  he  reported  a  “controversy  among  seholars  in  the  field”  [29]  regarding  the  issue. 
As  a  r  esult,  many  researchers  eurrently  utilize  aeroelastic  f luid-structural  i nteraction  m ethods, 
either  out  of  uncertainty  or  foregone  conclusion. 

Among  the  most  prominent  voiees  in  the  debate  regarding  aeroelastieity  (specifically  in 
Manduca  Sextae)  are  Combes  and  Daniel  [10,  11,  12,  13].  Combes  and  Daniel  carried  out  their 
studies  on  t  he  ba  sis  that  unde  rstanding  t  he  fluid-structural  i  nteractions  s  urrounding  t  he 
hawkmoth’s  forewing  helps  guide  the  modeling  process.  If  aeroelastic  effects  are  negligible  - 
that  is,  there  is  little  or  no  feedbaek  between  wing  deformations  and  fluid  mechanics  -  they  may 
be  omitted  during  modeling  [12].  This  omission  eonstitutes  a  dramatic,  desirable  simplification 
but  cannot  be  made  without  substantiation. 

Examining  the  degree  to  which  hawkmoths  are  affeetedby  aeroelastic  forces,  C ombes 
and  Daniel  performed  static  loading  tests  on  forewings  of  Manduca  Sextae,  as  shown  in  Eigure  5 
and  Eigure  6  [11,  12].  During  these  tests.  Combes  and  Daniel  used  a  pin  to  apply  a  point  load  to 
various  locations  on  the  wing,  and  collect  deflection  data.  Erom  these  defleetion  data.  Combes 
and  D  aniel  obs  erved  that  w  ing  stiffness  va ries  with  s panwise  1  oeation,  a nd  pr odueed  a  finite 
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element  model  to  represent  their  findings.  Inherent  in  their  methodology  is  the  assumption  that 
structural  (static)  stiffness,  El,  is  a  function  of  span.  Since  /,  the  moment  of  inertia,  is  a  function 
of  spanwise  length,  so  is  static  stiffness  -  even  if  E,  Young’s  modulus,  were  constant."^ 


Figure  5:  Static  test  setup  used  by  Combes  Figured:  Finite  element  model  developed  by 

and  Daniel  [10].  Combes  and  Daniel.  Color  bands 

denote  stiffness  as  a  function  of 
spanwise  distance  [11]. 

Following  their  stiffness  measurements.  Combes  and  Daniel  went  on  to  perform  flapping 
tests  in  air  and  helium  [11].  During  these  experiments,  hawkmoth  forewings  were  recorded  with 
high  speed  video.  Observing  their  footage.  Combes  and  Daniel  concluded  that  aeroelasticity  is  a 
function  of  wing  size  and  the  spanwise  stiffness  distribution  they  previously  calculated.  Based 
on  their  static  stiffness  tests.  Combes  and  Daniel  developed  a  parameter  to  describe  the  degree  to 
which  aeroelasticity  affects  various  flapping  wing  designs. 

During  t  heir  load  analysis,  C  ombes  a  nd  D  aniel,  c  iting  G  ordon  [  19],  u  sed  a  n  E  uler- 
Bernoulli  be  am  f  ormulation  toe  ompute  s  tiffness  a  s  a  f  unction  of  w  ingspan  [10].  S  everal 
assumptions  unde  rlie  E  uler-Bemoulli  t  heory:  beams  a  re  p  rismatic,  i  sotropic,  horn  ogeneous 
structures,  loaded  in  a  plane  of  symmetry,  with  plane  sections  remaining  plane.  Eurthermore, 

In  this  case,  it  is  more  likely  that  i?  is  a  weight  modulus  based  on  the  moduli  and  volume  fractions  of  the  veins  and 
membrane.  Thus,  the  wing  may  be  thought  of  as  a  composite  material,  obeying  the  rule  of  mixtures  [16]. 
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beams  are  eharaeterized  by  one  dimension  that  is  greater  than  the  other  two  by  at  least  one  order 
of  magnitude,  and  undergo  small  deformations  [18].  U  nfortunately,  inspecting  Manduca  Sexta 
wings  r  eveals  a  v  ery  d  ifferent  p  icture.  T  he  f  orewings  ar  e  n  on-prismatic,  non  -homogeneous 
plates  with  a  venation  pattern  resembling  composite  fibers.  N  oting  this  contrast.  Combes  and 
Daniel  acknowledged  the  possibility  of  a  “systematic  error  in  the  reported  values,  which  could 
potentially  underlie  the  observed  scaling  relationships”  [10]. 

In  m  any  w  ays,  C  ombes  a  nd  D  aniel  t  ook  a  r  emarkable  f  irst  s  tep  i  n  s  tudying  t  he 
mechanisms  behind  insect  flight.  C  ertainly,  their  work  is  widely  cited,  highly  publicized,  and 
often  de bated.  T heir  w ork  ha s  achieved  prestige,  as  well  as  been  mocked  as  conjecture  [9]. 
Gogulapati,  Friedman  and  Shy  argued  in  support  of  Combes  and  Daniel,  based  on  their  finite 
element  an  alysis  o  f  a  flapping  box-beam  [  20] .  Chen,  C  hen  a  nd  C  hou  a  rgued  a  gainst  C  ombes 
and  Daniel  on  the  basis  of  their  in  vacuo  photonic  frequency  analysis  of  a  dragonfly  wing  [9]. 

Fortunately,  a  pplying  c  lassie  s  tructural  d  ynamics  t  heory  s  erves  t  o  s  ettle  s  uch  an 
impassioned  debate.  While  Combes  and  Daniel  correctly  identified  the  need  for  flapping  tests  to 
determine  aeroelastic  effects,  their  work  was  based  upon  p  revious  static  loading  computations. 
As  previously  discussed,  the  validity  of  these  findings  i  s  qu  estionable,  due  to  the  s  implifying 
assumptions  inherent  in  Euler-Bernoulli  beam  theory.  Second,  static  loading  does  not  adequately 
characterize  a  structure’s  response  to  oscillations  -  structural  dynamics  theory  is  required. 

The  field  of  analytical  structural  dynamics  is  predicated  upon  de  termining  equations  of 
motion,  applying  an  input  (or  forcing  function),  and  solving  for  resultant  singularities  in  those 
equations  of  m  otion  [  23,  26].  While  c  omputational  me  thods  exist,  t  heir  pur  pose  i  s  t  o 
approximate  singularities  via  numerical  convergence,  where  an alytical  methods  are  intractable 
[14].  Combes  a  nd  D  aniel  di  d  i  nvoke  t  he  f  inite  element  m  ethod  ( FEM),  but  onl  y  ve  rify  t  heir 
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load-deflection  tests  [10].  Gogulapati,  Friedman  and  Shy  also  utilized  FEM,  but  their  objective 
was  to  study  fluid-structural  i  nteractions  (FSI);  due  to  the  limitations  of  computational  fluid 
dynamics  (CFD)  solvers,  they  modeled  a  beam  [20].  The  applicability  of  findings  based  on  a 
box-beam  model  to  a  Manduca  Sexta  forewing  is  questionable. 

In  an  attempt  to  eliminate  concerns  arising  from  the  above  methods,  Chen,  Chen,  and 
Chou  conducted  a  frequency  analysis  of  a  dragonfly  [9].  Unfortunately,  their  photogrammetric 
analysis  technique  required  them  to  paint  metallic  dots  on  t  he  wing,  adding  nearly  30%  to  its 
mass  (Figure  7).  Moreover,  Chen,  Chen,  and  Chou  failed  to  adequately  refute  the  aeroelastic 
effects  reported  by  Combes  and  Daniel,  since  their  vacuum  testing  eliminated  the  possibility  of 
aerodynamic  contributions.  Still,  Chen,  Chen,  and  Chou  deserve  credit  for  introducing  structural 
dynamics  testing  to  an  aeroelasticity  debate  dominated  by  biologists  and  entomologists. 


Figure  7:  (a)  Freshly  liberated  dragonfly  wing;  (b)  painted  wing 
used  by  Chen,  Chen,  and  Chou  [9]. 


1. 8.  Selected  Approach 

Borrowing  from  the  best  of  approaches  of  Combes  and  Daniel  [12]  and  Chen,  Chen,  and 
Chou  [9],  a  motor  was  used  to  flap  a  liberated  Manduca  Sexta  forewing  at  random  frequencies 
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(technically,  pseudo-random).  A  Polytec®  s  canning  1  aser  vi  brometer  was  us  ed  t  o  determine 
modes  and  modeshapes  ,via  fast  Fourier  transform  (FFT),  from  frequencies  approaching  zero  up 
to  one-thousand  Hertz.  Unlike  the  method  used  by  Chen,  Chen,  and  Chou,  the  laser  vibrometer 
offers  the  advantage  of  noninvasive  measurement.  The  specimen  remains  unaltered,  contributing 
neither  mass  nor  stiffness. 

Again  combining  the  methods  of  Combes  and  Daniel  and  Chen,  Chen,  and  Chou,  tests 
were  conducted  in  air  at  room  temperature  and  pressure,  then  repeated  in  vacuum.  Frequency 
response  functions  (FRFs)  of  hot  h  cases  were  compared  t  o  ev  aluate  t  he  extent  t  o  w  hich 
aeroelastic  forces  affect  flapping  wings. 

Computed  t  omography  (CT)  s  cans  pr  ovided  de  tailed  X  -Ray  i  mages  o  f  t  he  Manduca 
Sexta  vein  structure,  and  were  used  as  the  geometric  basis  for  modeling.  A  finite  element  model 
was  constructed  using  SolidWorks®,  Abaqus®^,  and  CoordGen^  [35].  Finally,  the  results  of  the 
frequency  analysis  in  Abaqus  were  compared  to  the  eigenvalues  and  eigenvectors  produced  by 
the  Polytec  scanning  laser  vibrometry  studies.  A  preliminary  investigation  into  what  level  of 
geometric  s  implification  mig  ht  b  e  p  ossible  when  manufacturing  a  synthetic  Manduca  Sexta 
forewing  was  also  conducted. 

Such  simplifications  are  desirable  for  three  reasons.  First,  manufacturing  costs  scale  with 
design  complexity.  Since  cost  reduction  is  one  of  many  attractive  aspects  of  future  MAVs,  trade 
studies  of  simplification  versus  performance  are  valuable  tools.  Second,  many  fluid-structural 
interaction  (FSI)  codes  do  not  interact  well  -  or  at  all  -  with  non -beam  elements.  The  simpler 
the  finite  element  model,  the  more  feasible  future  CFD  analyses  become.  Third,  the  proposed 
clandestine  nature  of  the  MAV  warrants  designing  many  types  of  synthetic  insect  look-alikes, 

^  “SolidWorks”  and  “Abaqus”  are  registered  trademarks,  and  are  proprietary  works  of  Dassault  Systemes  SA. 

^  “CoordGen”  is  copyrighted,  open-source  work,  published  under  license  in  Appendix  A  -  CoordGen. 
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since  not  all  insects  are  indigenous  to  all  parts  of  the  globe.  In  this  sense,  less  speeific  wing 
designs  o  ffer  in  creased  f  lexibility  in  application,  potentially  reducing  f  uture  research  an  d 
development  costs. 

1.9.  Document  Overview 

The  following  f  our  c  hapters  ar  e  a  rranged  i  n  ch  ronological  o  rder.  Chapter  1 1  discusses 
experimental  w  ork  c  ondueted  a  1 1  he  A  ir  F  oree  Institute  of  T  echnology  (  AFIT)  a  nd  at  t  he 
Materials  Directorate  of  the  Air  F oree  Researeh  Laboratory  (AFRL).  The  experimental  data  in 
Chapter  1 1  are  t  he  f  oundation  upon  w  hich  s  ubsequent  c  hapters  a  re  bui  It.  Chapter  III  details 
several  finite  element  and  numerieal  studies  conducted  before  proceeding  with  the  monumental 
task  of  producing  a  final,  finite  element  model.  These  analyses  are  custom-tailored  problems 
that  represent  the  charaeteristics  of  the  final  model,  but  are  far  more  tractable.  Chapter  IV  is  the 
results  and  discussion  section.  Here,  such  topics  as  aeroelastic  effects,  frequency  analysis,  tested 
biologieal  a  nd  s  ynthetie  s  pecimens,  a  nd  wing  s  tress  ar  e  p  resented.  Chapter  V  reviews  t  he 
highlights  of  eaeh  of  the  preeeding  chapters  to  present  conclusions  and  opportunities  for  future 
endeavors. 

The  appendices  follow  Chapter  V,  and  contain  supplemental  material,  from  whieh  one 
might  gain  further  insight  into  the  methods  used,  souree  code,  instructions  for  those  wishing  to 
conduct  similar,  future  researeh,  and  general  information  outside  the  primary  scope  of  this  work. 
Finally,  a  bibliography  eontains  any  references  used  in  preparation  of  this  doeument. 
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II.  Experimental  Methods 


The  p  rimary  f  ocus  during  the  experimental  phase  of  studying  the  Manduca  Sexta  was 
recording  the  natural  frequencies  and  associated  modeshapes  of  the  forewings.  In  an  effort  to 
reduce  an  y  possible  an  atomical  v  ariations,  onl  y  females  were  s  tudied.  Laser  v  ibrometry  w  as 
selected  as  t  he  pr  eferred  m  ethod  of  f  requency  extraction  due  t  o  i  ts  noni  nvasive  na  ture, 
repeatability,  and  rapid  test  times.  T  he  wings  were  mounted  to  a  motor,  exciting  the  wing  at 
pseudo-random  frequencies  up  to  approximately  400  Hz.  Time  domain  data  was  acquired  with 
velocity  decoders  and  converted  to  the  frequency  domain  via  the  fast  Fourier  transform  (FFT). 

Tests  were  conducted  in  air  and  in  vacuum  (less  than  530  Pa.,  or  approximately  1/150*  of 
an  a  tmosphere)  to  in  vestigate  th  e  p  ossible  a  eroelastic  c  ontributions  d  uring  flight.  Clamped 
boundary  conditions  were  verified  by  comparing  laser  scans  of  beams  to  analytical  solutions. 
The  results  of  each  of  these  tests  are  found  in  Table  2  and  Table  3. 

2.1.  Experimental  Setup 

The  frequency  analysis  of  a  Manduca  Sexta  forewing  was  conducted  using  the  Polytec 
PSV-400  scanning  laser  vibrometer.  Wings  were  mounted  into  a  foam-lined  clamp,  which  was 
bolted  1 0  a  m  otor.  T  he  m  otor  w  as  c  onnected  t  o  a  n  a  mplifier  t  hat  r  eceived  s  ignals  f  rom  t  he 
Polytec  function  generator.  The  motor  and  wing  were  encased  in  a  transparent,  acrylic  vacuum, 
angled  to  prevent  back-scatter  to  the  laser  head  (see  Figure  8).  During  vacuum  tests,  the  pressure 
was  h  eld  co  nstant  at  r  oom  t  emperature  and  ap  proximately  530  Pa.  During  te  sts  in  a  ir,  th  e 
vacuum  chamber  was  closed,  but  kept  at  room  temperature  and  pressure. 
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Figure  8:  Schematic  of  experimental  setup. 


Two  1  aser  he  ads  w  ere  u  sed  w  hile  condueting  frequeney  extraction  t  ests.  A  r  eference 
laser  focused  on  the  bolt-head  holding  the  wing  clamp  in  place,  while  the  other  laser  scanned  the 
wing,  as  shown  in  Figure  9.  The  scanning  laser  followed  a  pre-set  pattern  of  scan  points  over  the 
surface  of  the  wing.  This  scanning  point  cloud  was  sufficiently  dense,  and  an  average  of  five 
scans  per  point  was  recorded  and  processed  during  the  FFT.  During  the  tests,  the  motor  flapped 
the  wing  at  very  low  amplitudes^,  to  minimize  any  aeroelastic  influence  that  might  not  be  present 
during  normal  operation  by  the  Manduca  Sexta.  A  pseudo-random  signal  drove  the  motor  at 
approximately  0.2  vol  ts,  a  t  f  requencies  up  t  o  a  pproximately  400  Hz.  All  s  ignal  p  recessing, 
including  c  oherence  pi  ots  a  nd  frequency  response  f  unction  pi  ots,  w  as  au  tomatically  achieved 
with  Polytec’s  proprietary  software  suite. 


^  While  not  quantified,  these  amplitudes  were  small  enough  that  they  were  barely  visible  with  the  naked  eye. 
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Figure  9:  Wing  motor  and  laser  configuration  used  during  testing. 


Vacuum  t  ests  w  ere  co  nducted  i  nside  an  A  bbess  Instruments  s  tainless  s  teel  /  acr  ylic 
vacuum.  T  he  vacuum  chamber  is  a  cube,  transparent  on  five  sides,  each  measuring  24  inches. 
To  ensure  no  a  mbient  air  1  eaked  i  nto  t  he  va  cuum  dur  ing  t  esting,  a  di  gital  pr  essure  gage  i  s 
attached  to  the  vacuum,  and  monitored  during  testing.  The  vacuum  was  manually  operated,  and 
was  not  configured  to  automatically  compensate  for  any  (small)  changes  in  pressure;  however, 
the  vacuum  is  sufficiently  sealed  to  provide  quasi-equilibrium,  at  no  point  during  vacuum  testing 
did  the  pressure  exceed  the  4  Torr  self-imposed  limit.  Tests  in  air  were  held  at  room  temperature 
and  pressure,  by  allowing  the  vacuum  relieving  valve  to  remain  open. 
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Figure  10:  Abbess  Instruments  vacuum  chamber  used  during  testing. 


2.2.  Wing  Preparation 

Before  any  vibration  testing  can  be  conducted,  the  forewings  must  be  liberated  from  the 
Manduca  Sexta.  In  an  effort  to  minimize  any  unnecessary,  potential  suffering,  the  insects  were 
cooled  i  n  a  r  efrigerator  f  or  a  pproximately  t  hirty  minutes,  t  hen  checked  f  or  r  esponsiveness. 
Using  a  scalpel,  the  wings  were  removed  by  an  incision  parallel  to  the  body,  and  immediately 
distal  of  the  wing-thorax  adjacency  (shown  in  Figure  11).  During  this  process,  great  care  was 
taken  to  avoid  damaging  the  wing  prior  to  testing. 
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Figure  11:  Cut  location  during  wing  liberation  [21], 


Observed  during  the  wing  liberation  proeess  was  a  natural  camber  in  the  Manduca  Sexta 
forewing.  In  order  to  measure  the  frequency  characteristics  of  the  wing,  it  must  be  placed  inside 
a  rigid  c lamp.  U  nfortunately,  c lamping  the  wing  flattens  its  arched  c enter,  thus  changing  the 
shape  of  the  wing  and  producing  internal  stress.  These  changes  in  geometry  and  internal  stress 
alter  t  he  n  atural  f  requencies  an  d  m  ode  s  hapes,  an  u  nacceptable  r  esult  i  n  an  an  alysis  m  ethod 
selected  for  its  noninvasive  characteristics.  T  o  circumvent  this  unnatural  flattening  of  the  wing, 
a  polymeric  foam  paint  was  applied  inside  the  clamp  on  both  sides  of  the  wing,  shown  in  Figure 
12. 


©  0 


Figure  12:  Side  view  of  a  clamped  Manduca  Sexta  forewing  (a),  and  the  foam-lined  clamp  of 
cross-section  A-A  (b).  The  four  symbols  shown  in  (a)  are  reference  points  for  future 
photo-modeling  use. 
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2.3.  Verification  of  Boundary  Conditions 


The  Manduca  Sexta  wing  was  secured  by  a  foam-coated  clamp,  sufficiently  eliminating 
artificial  changes  to  its  natural  shape  (camber).  During  frequency  analysis  and  subsequent  finite 
element  modeling,  the  nature  of  the  applied  boundary  conditions  is  of  critical  importance  [14, 
26].  While  a dequately  maintaining  the  ha wkmoth  wing’s  undisturbed  shape,  a dding  pol ymer 
foam  causes  one  to  question  whether  t  he  boun  dary  c  ondition  is  s  ufficiently  clamped,  i  n  t  he 
analytical  sense.  An  easy  way  to  answer  this  question  is  to  conduct  a  frequency  analysis  of  a 
well-known  pr  oblem,  f  or  w  hich  ana  nalytical  s  olution  i  s  a  vailable  -  namely  a  n  i  sotropic, 
homogeneous,  pr  ismatic  be  am.  If  t  he  r  esults  f  rom  the  m  odal  an  alysis  of  a  b  earn  h  eld  b  y  t  he 
foam-lined  clamp  match  the  analytical  solution,  the  clamp  may  be  considered  rigid. 

2. 3. 1.  Frequency  Analysis  of  a  Beam 

A  beam  was  co  nstructed  b  y  cutting  a  rectangular  strip  from  a  piece  of  standard  copier 
paper.  ^  The  b  earn  w  as  p  laced  i  nside  t  he  f  oam-lined  c  lamp  i  ntended  f  or  1  ater  us  e,  and  a 
frequency  an  alysis  w  as  co  nducted  u  sing  t  he  P  olytec  s  canning  1  aser  s  canner.  T  he  t  est  w  as 
conducted  i  n  va  cuum  c  onditions  ( less  t  han  530  P  a.,  or  a  pproximately  4  T  orr)  t  o  a  void  any 
aeroelastic  contributions,  as  well  as  verify  the  ability  of  the  scanning  laser  to  operate  through  a 
transparent,  acrylic  vacuum.  The  frequency  response  function  (FRF)  is  plotted  below,  in  Figure 
13. 


*  Paper  used  was  Hewlett-Packard  Office  Recycled  Paper,  201b,  92  brightness. 
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Figure  1 3 :  Frequency  response  function  of  a  paper,  fixed-free,  cantilevered  beam. 

The  s  pectral  c  oherence  is  pi  otted  ve  rsus  f  requency  in  Figure  14.  S  pectral  co  herence 


satisfies  linearity  and  ergodicity^,  and  is  described  by 


^xy\ 


r  r 

^xx  ^yy 


(1) 


for  two  signals,  x(t)  and  y(t),  where  Gxy  is  the  cross-spectral  density  between  x(t)  and  y(t),  and 
Gxx  and  Gxy  are  t  he  au  tospectral  d  ensities  o  f  xfifj  and  respectively  [5].  S  ince  co  herence 
describes  the  degree  to  which  the  signals  x(t)  and  y(t)  are  related  in  a  linear  system,  it  also  serves 
as  an  i  ndicator  o  f  ef  fectively  calibrated  in  strumentation:  a  ssuming  a  s  tructure  is  lin  ear,  its 
vibration  measurement  and  reference  signals  should  produce  a  coherence  equal  to  one.  As  one 
might  expect,  the  beam  used  here  is  a  linear  structure  that  produces  adequate  values  of  coherence 
up  to  approximately  250  Hz. 


^  A  signal  is  ergodic  if  it  contains  sufficient  data  to  be  considered  stochastic.  In  so  doing,  its  statistical  properties 
may  be  computed  from  a  single  signal  that  is  recorded  for  a  sufficiently  long  time  [5]. 
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Figure  14:  Spectral  coherence  for  the  fixed-free  beam  analysis. 


2.3.2.  Analytical  Fixed-Free  Beam  Theory 

Euler-Bemoulli  b earn  theory  has  many  useful  applications  for  structures  that  experience 
negligible  shear  through  the  thickness  [34].  A  vibrating,  cantilevered  beam  of  constant  Young’s 
modulus,  E,  satisfies  this  basic  requirement  and  -  as  will  be  demonstrated  in  the  paragraphs  that 
follow  -  has  a  p  redictable  m  odal  behavior.  Most  not  ably,  t  he  na  tural  frequency  r  atios  of  a  n 
Euler-Bemoulli  be  am  a  re  a  Iways  t  he  s  ame,  s  o  1  ong  as  t  he  b  earn  p  roperties  a  nd  bounda  ry 
conditions  are  unchanged.  Asa  result,  the  experimentally  obtained  mo  des  o  f  any  prismatic, 
homogeneous,  i  sotropic  b  earn  -  for  w  hich  t  he  m  aterial  pr  operties  a  re  unknow  n  -  can  b  e 
compared  t  o  t  he  analytical  s  olution  t  o  as  sess  the  ef  fectiveness  o  f  t  he  real-world  boundary 
condition.  M  atched  a  nalytical  a  nd  e  xperimental  f  requency  r  atios  a  re  onl  y  pos  sible  i  f  t  he 
boundary  conditions  in  the  experiment  satisfy  the  criteria  of  the  analytical  end  conditions. 

The  e  quation  of  E  uler-Bemoulli  be  am  be  nding  i  s  g  iven  i  n  E  q.  (2)  [34],  w  here  E  is 
Young’s  modulus,  /  is  the  moment  of  inertia,  F  is  the  forcing  function,  M  is  the  bending  moment. 
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X  is  the  horizontal  coordinate,  y  is  the  vertical  coordinate,  t  is  time,  and  L  is  the  length  of  the 
beam.  A  schematic  of  the  beam  described  by  Eq.  (2)  is  shown  in  Figure  15  [26]. 


dx"^ 


EI(x) 


d^yjx,  t) 
dx'^ 


+  F{x,  t)  =  M(x) 


d^yjx,  t) 


0  <  X  <  L 


(2) 


Since  the  beam  is  homogeneous  and  prismatic.  Young’s  modulus,  the  moment  of  inertia, 
and  the  beam’s  mass  are  constant:  EI(x)=EI,  and  m(x)=m. 


b. 


f[x,t)dx 

ttl 

M(x,t)  -  1  ^  ox 


y(x.t)  Q(x.t) 


Figure  15:  Euler-Bemoulli  beam  in  bending  vibration  (a),  and  the  ifee- 

body  diagram  for  a  planar,  infinitesimal  beam  element  (b)  [26]. 


Inserting  the  homogeneous,  prismatic  assumptions  into  Eq.  (2),  and  assuming  F(t)  is  a 
harmonic  function,  yields  Eq.  (3),  where  O)  is  the  frequency  of  oscillation. 


dx'^ 


-p^Yix)^0, 


0  <  X  <  L] 


a)^m 

El 


(3) 


Applying  the  rigidly  clamped  boundary  condition,  the  displacement  and  rotation  are  both 
zero  at  the  wall: 
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X  —  0 


(4) 


=  0 


Likewise,  the  free  end  is  not  subjeeted  to  a  bending  moment  or  shear  stress,  thus: 

(d^Y{x)  _  ^ 

,  X  —  L 


dx'^ 


(5) 


d^Y(x) 


=  0 


dx^ 

Finally,  i  nvoking  t  he  know  n,  g  eneral  s  olution  of  E  q.  (3),  a  pplying  t  he  bounda  ry 
conditions  i  n  E  qs.  (4)  and  (5),  a  nd  c  onsiderable  a  Igebraic  s  implification  yields  th  e  c  anonical 
bending  equation  for  a  fixed-free  cantilevered  beam: 


cos(l3L)  coshipi)  —  —1  (6) 

Meirovitch  g  ives  t  he  t  ranscendental  e  quation  gove  rning  t  he  m  odeshapes  i  n  R  ef.  [  26],  a  nd 
provides  the  associated  frequencies  shown  in  Eigure  16.  Included  in  Eigure  16  are  the  frequency 
values  corresponding  to  the  numerical  solutions  of  Eq.  (6). 

An  important  facet  of  Eq.  (6)  is  that  for  a  prismatic,  isotropic,  homogeneous  beam,  the 
ratio  of  the  natural  frequencies  is  always  the  same  (provided  the  fixed-free  boundary  conditions 
are  satisfied).  In  other  words,  one  is  not  required  to  know  the  material  properties  (i.e.,  Young’s 
modulus,  E),  measure  the  beam’s  length,  or  calculate  the  moment  of  inertia  to  assess  whether  a 
beam  is  s  atisfactorily  c  lamped.  S  imply  c  omparing  th  e  r  atios  o  f  a  b  cam’s  me  asured  n  atural 
frequencies  to  the  ratios  given  in  Eigure  16  provides  an  excellent  metric  for  comparison. 
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Figure  16:  The  first  two  natural  frequeneies  and  associated 

modeshapes  for  a  fixed-free,  cantilevered  beam  [26]. 


2.3.3.  Comparison  of  Analytical  and  Experimental  Values 

The  results  of  the  frequency  extraction  are  shown  below  in  Table  2.  Table  3  compares 
the  ratios  of  the  paper  beam’s  frequency  values  to  those  obtained  analytically. 

Table  2.  Natural  frequencies  of  a  fixed-free  paper  beam. 

Frequency,  Hz. 

uji  15.625 

^2  100.0 

Table  3.  Natural  frequency  ratios  of  analytical  and  experimental  beams. 

Analytical  Solution  Paper  Beam  Error 

c<j2/c<Ji  6.267  6.40  2.08  % 

The  values  given  in  Table  3  lead  one  to  the  conclusion  that  the  paper  beam  is  made  of  a 
homogeneous,  isotropic  material,  and  sufficiently  satisfies  the  clamped  boundary  condition In 
so  doing,  the  simple  exercise  of  conducting  a  frequency  analysis  of  a  strip  of  paper  produced 
important  r  esults.  O  ne  can  now  s  ay  w  ith  c  onfidence  t  hat  the  pol  ymer  foam  pa  inted  ont  o  t  he 

***  Prior  to  this  experiment,  it  was  unknown  whether  copy  paper  was  isotropic  and  quasi-homogeneous. 
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clamp  does  not  adversely  affeet  its  performanee  as  a  fixed  boundary  eondition,  without  ehanging 
the  g  eometry  of  a  mounted  Manduca  Sexta  wing.  Moreover,  i  t  i  s  evident  t  hat  t  he  P  olytec 
seanning  laser  vibrometer  is  eapable  of  testing  through  the  transparent  vacuum  used  in  the  setup 
described  above. 

2.4.  Test  Flow  and  Timeline 

With  t  he  e  lamped  boun  dary  e  ondition  v  erified,  a  nd  t  he  s  uecessful,  opt  ieal  f  requeney 
analysis  through  aerylie  eompleted,  the  proeess  of  testing  a  Manduca  Sexta  wing  begins.  Wing 
testing  is  a  multi-step,  time-eonsuming  proeess.  Unfortunately,  liberating  the  wing  exposes  the 
eireulatory  members  at  the  wing  root,  eausing  the  wing  to  dry.  As  the  wing  dries,  its  properties 
ehange,  nullifying  the  frequeney  data  and  any  eonelusions  drawn  from  them.  As  part  of  his  work 
at  AFIT,  Norris  eonducted  an  investigation  of  Manduca  Sexta  wing  drying  effeets.  During  this 
exereise,  h  e  aeq  uired  f  requeney  d  ata  f  rom  m  ultiple  lib  erated  w  ings,  a  nd  c  ompared  th  eir 
responses  to  the  liberation  timeline.  He  eoneluded  that,  under  normal,  indoor  conditions,  wings 
were  still  “healthy”  up  to  approximately  four  hours  from  separation  [30]. 

Operating  within  this  four-hour  window,  two  main  studies  were  eondueted.  C  omputed 
tomography  (CT)  scans  captured  the  geometry  of  the  wing  members.  Immediately  following  the 
CT  i  maging,  t  he  w  ing  w  as  p  laeed  i  nside  t  he  v  aeuum  eh  amber  an  d  a  nalyzed  u  sing  1  aser 
vibrometry. 

2.5.  Computed  Tomography  (CT)  Imaging 

The  Manduca  Sexta  forewing  is  a  highly  eomplex  strueture.  There  are  a  multitude  of 
veins  and  membranes,  many  of  whieh  intersect  multiple  times.  Beeause  of  this  eomplexity,  a 
high-fidelity  geometrie  mo  del  is  r  equired  to  p  roeeed  w  ith  a  ny  f  inite  e  lement  e  fforts. 
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Unfortunately,  t  he  s  eale  a  nd  i  ntricacy  of  t  he  w  ing  limit  th  e  n  umber  of  deviees  c  apable  o  f 
aceurately  resolving  its  components.  O  ne  apparatus  (among  ot  hers)  t  hat  is  cap  able  o  f 
characterizing  the  geometry  of  the  wing  is  a  computed  tomography  (CT)  “scanner.” 

Computed  tomography  scanners  rotate  an  X-Ray  imager  around  a  specimen,  while  taking 
a  s  erics  oft  wo-dimensional  X  -Ray  i  mages.  T  hese  t  wo-dimensional  i  mages  a  re  t  hen  pos  t- 
processed  t  0  bui  id  a  three-dimensional  digital  model  [22].  The  r  esolution  of  the  model  is 
strongly  dependent  on  t  he  number  of  images  taken,  however  increasing  the  number  of  X-Ray 
images  r  esults  in  lo  nger  s  can  time  s.  Fortunately,  C  T  i  s  a  popul  ar  i  maging  t  echnique  i  n  t  he 
medical  c  ommunity,  th  us  ma  ny  of  its  p  ost-processing  f  eatures  ar  e  optimized  for  e  xamining 
organic  tissue  -  including  moth  wings. 

2.5.1.  Imaging  a  Manduca  Sexta  Wing 

A  freshly  liberated  wing  was  clamped,  as  shown  in  Figure  12,  and  placed  inside  the  CT 
scanner  at  the  Materials  Directorate  of  the  Air  Force  Research  Laboratory  (AFRL/RX).  Previous 
tests  by  Norris  resulted  in  high-fidelity  CT  models  after  only  eight  minutes  [30].  Choosing  to  err 
on  the  side  of  caution,  the  Manduca  Sexta  test  specimen  was  scanned  for  20  m  inutes.  A  plan 
view  of  the  wing  X-Ray  is  shown  below,  in  Figure  17. 

The  nature  of  multi-image  composite  modeling  used  in  computed  tomography  allows  one 
to  vi  ew  c  ross-sectional  “  slices,”  t  hrough  va  rying  t  issue  de  nsities.  T  his  be  comes  a  n  i  mportant 
feature  w  hen  a  ttempting  t  o  di  mension  t  he  c  omponents  i  n  t  he  ha  wkmoth’s  w  ing.  P  reviously, 
others  ha  ve  a  ssumed  t  he  ve  ins  w  ere  c  ircular  i  n  c  ross-section  a  nd  measured  v  ein  w  idths 
(diameters)  from  above  [25].  A  quick  glance  at  Figure  17  reveals  the  difficulty  one  might  have 
in  such  a  process.  Where  does  the  vein  stop,  and  the  membrane  begin?  How  large  is  the  inner 
diameter  of  each  vein  (they  are,  after  all,  circulatory  members)? 
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Figure  17:  Unprocessed  CT  image  of  a  female  Manduca  Sexta  ’s  right  forewing. 
Changing  the  density  threshold  in  the  CT  model  helps  answer  these  questions. 


Figure  18  shows  the  same  wing  as  Figure  17,  capturing  only  high-density  tissue.  Note 
that  in  th  is  ima  ge,  th  e  in  ner  a  nd  o  uter  d  iameters  of  t  he  v  eins  a  re  clearly  visible  upon  c  lose 
inspection.  L  ikewise,  the  boundary  between  vein  and  membrane  is  much  easier  to  distinguish. 
Perhaps  one  of  the  most  obvious  artifacts  of  only  resolving  the  higher-density  structures  is  the 
noticeable  lack  of  membrane  tissue.  It  appears  that  the  wing  skin  is  of  equal  density  at  all  points, 
possibly  i  ndicating  horn  ogeneity  -  casting  doub  t  on  previous  assertions  o  f  v  arying  ma  terial 
properties  [10,11,  25]. 
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2.5.2.  Wing  Geometry 


Ultimately,  the  goal  of  utilizing  CT  imaging  techniques  is  capturing  the  geometry  of  the 
wing  CO  mponents.  O  nee  co  mpleted,  t  hese  m  easurements  ar  e  u  sed  toe  reate  a  finite  el  ement 
model  of  t  he  w  ing.  S  ince  t  he  C  T  mo  del  is  th  ree-dimensional,  it  is  p  ossible  to  me  asure  th  e 
dimensions  of  t  he  ve  ins  a  nd  m  embranes.  T  his  i  s  a  ccomplished  b  y  m  easuring  c  ross-sectional 
“slices”  of  the  wing,  from  root  to  tip.  Figure  19  shows  one  such  slice,  located  approximately  two 
millimeters  from  the  wing  root^^ 

Consisting  of  a  total  of  five  veins,  the  image  illustrates  the  complexity  of  the  Manduca 
Sexta  vein  structure  -  no  two  veins  are  alike,  though  all  are  approximately  circular.  Furthermore, 
the  vein  closest  to  the  trailing  edge  (right-most  in  Figure  19)  is  dramatically  different  than  the 
others  at  this  particular  cross-section  location  (Figure  18,  A- A).  Here,  two  veins  are  merging 
into  one,  larger  vein. 


Figure  19:  Cross-section  A-A  (see  Figure  18)  from  a  Manduca  Sexta  wing,  approximately  2 
mm.  from  the  root.  The  prominent  left  vein  is  the  leading  edge. 

With  the  vein  and  membrane  dimensions  clearly  visible,  measurements  were  taken  at  16 
locations,  or  approximately  every  six  percent  of  the  46  c  m.  spanwise  length  .  Cross-sectional 
slices  were  concentrated  at  areas  where  the  vein  diameter  changed  most  rapidly  (especially  near 
the  tip).  Although  the  vein  diameters  tapered  with  length,  their  cross-sectional  shapes  remained 


*’  For  readability’s  sake,  grey  pixels  within  3%  of  white  (denoted  in  hexadecimal  color  code  as  #FFFFFF)  were 
converted  to  white.  The  original,  unaltered  image  is  located  in  Appendix  B  -  CT  Images. 

The  exact  locations  are  given  in  Appendix  B  -  CT  Images. 
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relatively  constant  between  measurement  locations.  Only  at  the  locations  where  veins  merge  or 
branch  were  irregularities  observed. 

All  veins  were  observed  to  be  relatively  circular  in  shape,  with  one  exception.  N ear  the 
midpoint  of  the  leading  edge,  several  veins  merge  into  one  structure;  during  this  transition,  the 
two  leading  edge  veins  gradually  become  closer,  until  the  two  circular  structures  merge  into  a 
binocular-shaped  doubl  e  pi  pe,  then  finally  into  a  b  ox -beam.  At  a  pproximately  60%  of  t  he 
spanwise  length,  the  box -beam-shaped  leading-edge  vein  resumes  its  circular  shape.  Unlike  the 
veins,  however,  the  membrane  remained  a  nearly  constant  0.1 2  millimeter  thickness,  deviating 
only  ±0.02  millimeters  at  the  extrema. 

2.6,  Modal  Identification  of  a  Manduca  Sexta  Forewing  in  Vacuum 

A  Manduca  Sexta  wing  was  analyzed  using  a  scanning  laser  vibrometer  to  determine  its 
frequency  c  ontent.  As  described  p  reviously,  t  he  f  requency  d  ata  w  ere  collected  in  vacuo,  t  o 
eliminate  an  y  p  otential  aer  oelastic  ef  fects.  T  he  au  tomatically-computed  f  requency  r  esponse 
function  for  the  Manduca  Sexta  forewing  is  shown  below,  in  Figure  20.  Table  4  gives  the  first 
three  natural  frequencies  observed  in  this  experiment. 

Table  4.  Natural  frequencies  of  a  Manduca  Sexta  forewing  in  vacuum. 

Frequency,  Hz 

0)^  86 

0)2  106 

0)3  155 

Notice  t  hat  t  he  co  herence  d  ecreases  as  t  he  f  requency  ap  proaches  t  he  400  H  z  cu  toff. 
Since  t  his  a  nalysis  i  s  p  rimarily  concerned  w  ith  onl  y  t  he  first  t  hree  m  odes,  t  he  r  eduction  i  n 
coherence  values  at  high-frequencies  do  not  negatively  impact  the  fidelity  of  the  low-frequency 
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Coherence,  ND  Magnitude,  nnm-s'ViTim-s 


(less  than  250  Hz)  measurements.  These  high  coherence  values  seem  to  indicate  that  the  wing  is 
a  linear  structure,  since  any  structural  response  is  strongly  correlated  to  its  excitation. 


Frequency,  Hz 

Figure  20:  Frequency  r  esponse  function  ( magnitude) 
for  a  Manduca  Sexta  wing  in  vacuum. 


T 
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Frequency,  Hz 

Figure  21 :  Spectral  coherence  for  the  Manduca  Sexta  wing  analysis. 
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The  m  odeshapes  -  how  t  he  s  tructure  responds  t  o  a  na  tural  f  requency  i  nput  -  are 
important  parameters  in  describing  the  wing  .  Due  to  the  orthogonality  of  the  eigenvalues  and 
eigenvectors,  each  of  the  natural  frequencies  (modes)  listed  in  Table  4  has  a  unique,  associated 
modeshape  [26].  The  modeshapes  corresponding  to  the  first  three  natural  frequencies  are  shown 
below.  For  ease  of  understanding,  all  mode  shapes  are  presented  first  in  an  isometric  view  with 
the  trailing  edge  in  front  and  the  leading  edge  in  the  rear.  Each  modeshape  is  also  accompanied 
by  a  tip  view,  looking  inward  towards  the  wing  root. 

The  first  natural  frequency  occurs  at  86  Hz,  and  exhibits  a  modeshape  analogous  to  the 
first  bending  mode  of  a  cantilevered  b earn.  Similar  to  the  beam  shown  by  Meirovich,  the  first 
modeshape  of  the  Manduca  Sexta  forewing  has  a  maximum  displacement  at  the  tip,  and  flaps  in 
only  one  plane  of  motion. 


(a)  (b) 

Figure  22:  First  modeshape  (86  Hz)  of  a  Manduca  Sexta  wing,  shown 
from  an  isometrie  view  (a)  and  a  tip  view  (b). 


The  second  mode  of  the  Manduca  Sexta  forewing  occurs  at  106  H  z  with  a  modeshape 

similar  to  the  first  torsional  mode  of  a  plate  fixed  at  one  end.  The  maximum  displacement  is  at 

the  wing  tip,  with  the  leading  and  trailing  edges  deflecting  asynchronously.  Also  akin  to  the 

second  mode  of  a  plate  fixed  at  one  end  (first  torsion),  the  second  modeshape  of  the  hawkmoth 

wing  rotates  about  a  central  axis,  approximately  located  at  the  half-chord  line. 

It  s  hould  ben  oted  t  hat  t  he  term  “modeshape”  i  s  used  here  tod  enote  an  o  perational  d  eflection  shape  t  hat  i  s 
strongly  dominated  by  the  wing  response  at  resonance. 
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(a)  (b) 

Figure  23:  Second  modeshape  (106  Hz)  of  a  Manduca  Sexta  wing, 
shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 


The  t  bird  m  ode  oc  curs  at  1 55  H  z.  In  m  any  w  ays,  t  he  be  havior  of  t  he  w  ing  a  t  th  is 
frequency  i  s  s  imilar  tot  hat  of  i  ts  s  econd  m  ode,  w  ith  a  n  i  mportant  e  xception.  Instead  o  f 
asynchronous  flapping  a  bout  t  he  ha  If-chord  line,  t  his  mode  di  splays  s  ynchronous  (  in-phase) 
flapping  of  the  leading  and  trailing  edges  at  the  tip.  As  this  occurs,  the  half-chord  line  exhibits 
bending  similar  to  the  first-beam  bending  mode,  but  out  of  phase  with  the  flapping  in  the  leading 
and  trailing  edges.  The  combination  of  the  leading  /  trailing  tip  flapping  and  half-chord  bending 
creates  an  overall  deflection  shape  resembling  a  saddle.  This  “saddle  mode”  is  shown  below,  in 
Figure  24.  A  similar  saddle  mode  is  seen  in  the  fifth  modeshape  of  a  rectangular  plate  fixed  at 
one  end. 


(a)  (b) 

Figure  24:  Third  modeshape  (155  Hz)  of  a  Manduca  Sexta  wing, 
shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 


34 


Examining  t  he  f  requency  r  esponse  f  unction  r  eveals  t  hat  t  he  m  agnitudes  of  high- 
frequency  modes  are  small,  with  the  exception  of  368  Hz  mode.  Since  the  Manduca  Sexta  beats 
its  w  ings  ne  ar  25  Hz.  [40],  it  is  unlikely  that  the  h  igher-frequency  m  odes  ha  ve  m  ore  t  ban 
miniscule  effects  on  the  operational  deflection  shape  of  the  wing.  T  he  fourth  and  higher  mode 
shapes  are  shown  in  Appendix  C  -  Higher  Order  Modes,  but  are  omitted  from  further  analysis. 
Every  continuous  structure  has  an  infinite  number  of  natural  frequencies,  forcing  one  to  define 
an  uppe  r  bound  during  a  nalysis  [  26].  Considering  onl  y  t  he  first  t  hree  m  odes  provides  a  n 
adequate  basis  for  future  model  verification,  as  well  as  limits  the  scope  of  the  investigation. 

2. 7.  Modal  Identification  of  a  Manduca  Sexta  Forewing  in  Air 

The  Manduca  Sexta  frequency  identification  e  xperiment  was  r  epeated  using  t  he  s  ame 
setup  and  methods,  only  at  room  temperature  and  pressure.  To  eliminate  drying  effects,  a  fresh 
wing  was  used.  The  modeshapes  of  the  wing  were  the  same  as  those  given  above;  however  the 
frequencies  at  which  they  occurred  were  different.  Evidently,  the  fluid  interaction  with  the  wing 
serves  as  a  damping  mechanism,  reducing  the  frequency  content  by  as  much  as  26%  at  the  lower 
modes.  Eikewise,  th  e  ratios  o  f  th  e  n  atural  f  requencies  changed,  an  d  seemed  t  o  s  calc  with 
increasing  frequency. 

This  f  inding  provides  insight  i  nto  t  he  pr  eviously  m  entioned  -  and  m  uch  de  bated  - 
aeroelastic  question.  E urthermore,  future  studies  regarding  the  nature  of  the  damping  involved 
will  likely  prove  beneficial,  since  the  Manduca  Sexta  does  not  operate  in  vacuum  conditions  in 
nature.  Table  5  lists  t he  f requencies  as sociated  w ith  each  mode,  and  Table  6  compares  the 
frequency  ratios. 


35 


Tables.  Natural  frequencies  of  a  Manduca  Sexta  wing  in 


air  and  vacuum. 

Air,  Hz 

Vacuum,  Hz 

Change 

0)1 

58.75 

80 

26.56% 

0J2 

75 

98 

23.47% 

OJ3 

95 

118 

19.49% 

Table  6. 

Comparison  of  frequency  ratios  of  Manduca  Sexta 
and  wing  in  air  and  vacuum. 

Air 

Vacuum 

Change 

1.28 

1.23 

4.04% 

0)3/0)! 

1.62 

1.80 

8.78% 

2.8.  Paper  Wing  Tests 

In  1  ight  of  t  he  c  amber  obs  erved  during  t  he  e  omputed  t  omography  (CT)  s  eans,  one 
wonders  w  hether  e  amber  i  s  an  i  mportant  s  truetural  eh  araeteristic  i  n  t  he  h  awkmoth  wing. 
Certainly,  it  is  a  widely  aceepted  faet  that  eamber  inereases  low-Reynolds  number  aerodynamie 
efficiency  [33]  -  a  valuable  attribute  for  small  wings.  However,  to  date,  flat,  rigid  wing  designs 
with  homogeneous  m  aterial  account  f  or  a  1  arge  number  o  f  M  AV  r  esearch  [  29] .  These  w  ings 
provide  computational  simplicity,  and  are  cost  effective. 

Flat,  rigid  flapping  wings  are  often  relevant  in  the  scientific  endeavor.  In  fact,  several 
studies  ha  ve  f  ound  r  igid,  pi  ate-like  w  ings  on  t  he  Calliphora  (blowfly),  Tibicen  canicularis 
(cicada),  and  others  [25].  Since  the  Manduca  Sexta  falls  inside  the  D  ARPA  size,  weight,  and 
payload  r  equirements,  i  t  i  s  be  neticial  t  o  know  whether  a  flat,  horn  ogeneous  w  ing  r  easonably 
approximates  the  actual  forewing.  While  one  could  construct  a  finite  element  model,  a  simple 
process  i s  1 0  m ake  a p  aper  replica,  conduct  a  frequency  analysis,  and  compare  the  results  to 
Manduca  Sexta  forewing  data. 
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Previous  tests  showed  that  copier  paper  is  a  h  omogeneous  material,  suitable  for  testing 
simple,  semi-two-dimensional  cases.  T  o  explore  the  extent  to  which  the  veins  and  membranes 
are  quasi-isotropic,  and  the  degree  to  which  camber  effects  the  hawkmoth  wing  modeshapes,  a 
paper  wing  of  the  same  planform  was  constructed.  To  simplify  the  process,  a  Manduca  Sexta 
wing  was  placed  inside  a  photocopier;  the  copy  was  cut  out  in  the  shape  of  the  printed  wing  and 
placed  inside  the  vacuum  chamber  for  testing Figure  25  shows  the  wing  used,  and  the  results 
of  this  test  are  shown  in  Figure  26  and  Table  7. 


Clamped  section 


Figure  25:  Paper  wing 
used  in  frequency 
experiments. 


Frequency,  Hz 

Figure  26:  Frequency  response  function  of  the  paper  wing. 


14 


A  tacit  assumption  here  is  that  the  ink  on  the  copier  paper  does  not  significantly  affect  the  “wing”  response. 
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Table  7.  Natural  frequencies  of  a  paper  wing  in  vacuum. 


Frequency,  Hz 


20 

0)2 

96 

0)3 

162 

(U4 

472 

"5 

692 

Table  8  compares  the  ratios  of  the  natural  frequencies  obtained  during  the  paper  wing 
tests  to  those  of  the  Manduca  Sexta  in  vacuum.  Clearly,  the  two  are  unrelated  in  their  frequency 
content,  and  for  reasons  other  than  their  difference  in  materials  and  mass.  While  one  would 
expect  their  natural  frequencies  to  be  different,  the  hawkmoth  wing  and  paper  wing  would  have 
similar  frequency  ratios  if  material  properties  and  mass  were  their  only  differences.  It  is  likely 
that  the  Manduca  Sexta  wing  stiffening  due  to  the  presence  of  circulatory  members  is  affecting 
its  frequency  content. 


Table  8. 

Comparison  of  frequency  ratios  of  Manduca 
Sexta  and  paper  wings  in  vacuum. 

Manduca  Sexta  Wing 

Paper  Wing 

0)2l0)\ 

1.23 

4.8 

6J3/6J1 

1.80 

8.1 

0)^1 0)-^ 

2.17 

23.6 

6J5/6J1 

2.69 

34.6 

Examining  t  he  m  odeshapes  of  t  he  pa  per  w  ing  r  eveals  a  d  ifferent  p  icture.  In  f  act,  the 
modeshapes  are  similar  to  the  Manduca  Sexta,  even  though  their  frequency  content  is  radically 
different  -  suggesting  wing  geometry  is  a  dominant  factor  in  modeshape  appearance.  Although 
at  different  frequencies  and  a  different  frequency  ratio,  the  first  two  modeshapes  of  the  paper 
wing  are  nearly  identical  to  those  of  the  hawkmoth.  The  first  mode  is  a  flapping  mode,  similar  to 
a  homogeneous  fixed-free  beam’s  first  bending  mode  (see  Figure  27).  The  second  modeshape  is 
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similar  to  the  second  mode  of  a  flat,  rectangular  plate  fixed  at  one  end;  exhibiting  torsion  about 
the  central  axis,  from  root  to  tip  (see  Figure  28). 


Figure  27:  First  modeshape  (20  FIz)  of  a  paper  Manduca  Sexta  wing, 
shown  from  an  isometrie  view  (a)  and  a  tip  view  (b). 


(a)  (b) 

Figure  28:  Seeond  modeshape  (96  Hz)  of  a  paper  Manduca  Sexta 

wing,  shown  from  an  isometrie  view  (a)  and  a  tip  view  (b). 


The  third  modeshape,  shown  in  Figure  29,  appears  at  162  Hz.  Here,  the  wing  appears  to 
approximate  the  second  bending  mode  of  a  fixed-free  beam  at  the  half-chord  line.  The  leading 
and  trailing  edges  appear  to  lag  the  phase  of  the  half-chord  line,  but  mimic  its  behavior.  The 
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third  m  odeshape  is  uni  que  t  o  a  paper  w  ing,  suggesting  t  he  Manduca  Sexta  wing  h  as  m  ore 
chordwise  stiffness  than  its  paper  counterpart. 


(a)  (b) 

Figure  29:  Third  modeshape  ( 1 62  Hz)  of  a  paper  Manduca  Sexta 

wing,  shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 


As  a  CO  unter-argument,  one  m  ay  s  uggest  t  hat  t  he  t  bird  m  odeshape  i  n  t  he  pa  per  w  ing 
might  b  e  p  resent  in  the  Manduca  Sexta  wing,  but  out  of  order,  z.e.,  i  t  ap  pears  at  higher 
frequencies.  T  his  i  s  uni  ikely,  how  ever,  s  ince  t  he  hi  gher  f requency  m  odes  s  eem  toe  xcite  th  e 
entire  structure  less  with  increasing  frequency.  The  high-order  modeshapes  found  Appendix  C  - 
Higher  Order  Modes  paint  a  clear  picture:  high-frequency  content  is  seen  most  at  the  wing  tip 
and  trailing  edge,  near  the  half-span. 

The  fourth  paper  wing  mode  appears  at  472  Hz,  and  approximates  the  saddle  mode  (third 
mode)  of  the  hawkmoth  wing,  but  with  added  tip-deflection.  The  leading  and  trailing  edges  flap 
synchronously,  and  out  of  phase  with  the  vertical  motion  of  the  half-chord  line.  An  important 
distinction  be  tween  t  he  Manduca  Sexta  wing  and  t  he  pa  per  m  odel  i  s  t  he  t  ip  be  havior.  A 
comparison  of  Figure  24  (b)  and  Figure  30  (b)  highlights  the  increased  tip  motion  in  the  paper 
specimen,  relative  to  the  rest  of  the  wing.  The  third  modeshape  of  the  Manduca  Sexta  wing 
showed  asynchronous  deflections  in  the  trailing  edge  and  the  tip  of  the  wing,  with  only  slight 
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leading  edge  movement.  In  contrast,  the  fourth  modeshape  of  the  paper  wing  shows  prominent 
tip  deflection  in  phase  with  the  central  portion  of  the  wing.  While  this  characterization  is  purely 
qualitative  i  n  na  ture,  i  t  provides  us  eful  i  nsight  i  nto  w  hat  f  eatures  a  re  i  mportant  i  n  t  he  modal 
behavior  of  a  Manduca  Sexta  wing. 


(a)  (b) 

Figure  30:  Fourth  modeshape  (472  Hz)  of  a  paper  Manduca  Sexta 

wing,  shown  from  an  isometrie  view  (a)  and  a  tip  view  (b). 


Combined  with  the  flexibility  of  the  leading  edge,  the  vertical  deflection  of  the  wing  tip 
in  Figure  30  appears  to  be  much  more  dominant  than  that  observed  in  the  Manduca  Sexta  wing; 
this  flexibility  in  the  paper  specimen  likely  results  from  the  lack  of  veins  or  cambering  stiffness. 
Tests  w  ere  co  nducted  in  vacuo,  1  eading  one  t  o  pos  tulate  t  hat  t  he  pr  esence  of  ve  ins  i  ncreases 
kinematic  tip-stiffness. 

The  fifth  mode  is  given  in  Appendix  C  -  Higher  Order  Modes.  As  described  above,  the 
higher-order  modes,  while  useful  for  model  validation  play  a  miniscule  role  in  the  behavior  of 
the  Manduca  Sexta  wings.  This  is  even  more  the  case  for  the  paper  specimen,  since  the  fifth 
mode  occurs  at  692  Hz  -  28  times  the  beat  frequency  of  the  Manduca  Sexta  wing.  Inspecting  the 
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frequency  response  function  of  Figure  26  reveals  that  the  first  three  modes  are  dominant,  with  a 
significant  loss  in  power  in  the  fourth  and  higher  modes. 

Comparing  the  actual  Manduca  Sexta  wing  and  the  flat,  paper  specimen  gives  insight 
regarding  the  importance  of  the  wing  components.  The  hawkmoth  wing  appears  to  be  dominated 
by  plate-like  modeshapes  -  as  one  might  expect  given  its  dimensions.  By  contrast,  the  paper 
wing  s  eems  to  e  xhibit  a  c  omposite  b  earn  /  p  late  b  ehavior,  in  termittently  r  esembling  either 
structure.  The  absence  of  the  second  beam  bending  mode  from  the  Manduca  Sexta  wing  and  the 
lack  of  tip  motion  in  the  saddle  mode  suggests  the  veins  have  structural  importance  -  that  is,  they 
are  not  only  useful  as  circulatory  members.  At  this  point,  it  is  evident  that  mimicking  Manduca 
Sexta  flight  will  require  stiffening  members,  though  it  is  unclear  to  what  extent  camber  plays  in 
the  wing  response. 
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III.  Numerical  Methods 


The  primary  focus  during  the  numerical  phase  of  study  is  identifying  features  of  the  finite 
element  method  that  will  prove  useful  in  eventually  modeling  the  Manduca  Sexta  forewing.  The 
computed  t  omography  ( CT)  i  mages  pr  ovided  t  he  ge  ometric  b  asis  f  or  m  odeling,  a  nd  t  he  1  aser 
vibrometry  provided  the  wing  modes  and  modeshapes.  T  ogether,  these  methods  provide  a  vast 
amount  of  i  nformation  about  t  he  ha  wkmoth  w  ing  s  tructure.  I  ndeed,  r  everse  en  gineering  a 
Manduca  Sexta  wing  is  no  trivial  task,  and  requires  one  to  make  simplifications  in  the  modeling 
process.  This  chapter  explores  what  simplifications  are  permissible,  what  attributes  must  closely 
mimic  the  Manduca  Sexta,  and  how  one  might  proceed  with  the  modeling  process. 

Such  t  opics  as  el  ement  s  election,  g  eometric  simplification,  c  amber  effects,  me  sh 
convergence,  and  modal  effects  of  vein  diameter  are  discussed  in  the  sections  that  follow.  All 
modeling  is  carried  out  in  Abaqus  Standard,  Version  6.8.  From  these  analyses  it  is  clear  that  the 
veins  should  be  m  odeled  us  ing  qua  dratic  be  ams,  a  nd  t  he  w  ing  m  embrane  s  hould  be  m  odeled 
using  linear,  g  eneral-purpose  shell  el  ements.  A  c  onvergence  s  tudy  o  f  t  he  f  inal  m  odd  i  s 
presented  in  Chapter  IV. 

This  chapter  begins  with  a  discussion  of  model  construction  using  the  data  collected  from 
the  CT  scans.  A  discussion  of  element  selection  follows,  with  convergence  studies  for  each  case. 
Following  the  results  of  the  studies  involving  the  membrane  elements,  a  discussion  of  membrane 
element  t  heory  explains  t  he  f  ormulation  o  f  pi  ane-stress  elements,  a  nd  their  s  hortcomings  in 
modeling  Manduca  Sexta  wings.  Abaqus  is  validated  through  a  comparison  of  experimental  and 
finite  e  lement  results  i  nvolving  a  p  aper  wing.  Finally,  s  tudies  are  c  onducted  e  xamining  t  he 
effects  of  wing  camber  and  vein  dimensions  on  natural  frequency  in  the  Manduca  Sexta  wing. 
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3.1.  Model  Construction 


The  c  omputed  t  omography  ( CT)  s  cans  pr  ovided  a  n  e  xcellent  t  hree-dimensional  d  igital 
model  of  the  moth  wing.  Since  this  model  is  built  directly  from  X-Ray  images  of  the  wing,  it  has 
only  one  geometric  simplification.  The  wing  surfaces  are  exported  to  stereolithography  (STL) 
CAD  files,  with  continuous  surfaces  approximated  by  several  triangular  facets.  These  STL  files 
can  be  i  mported  i  nto  va  rious  C  AD  pr  ograms  a  nd  t  hen  i  nto  A  baqus  for  an  alysis.  There  ar  e 
several  p  otential  b  enefits  in  u  sing  th  is  method  for  f  inite  e  lement  a  nalysis:  th  e  STL  mo  del 
approximates  t  he  act  ual  s  tructure  w  ell,  nog  eometry  m  easurements  ar  e  n  ecessary,  an  d  m  odel 
construction  -  one  of  the  most  time-consuming  tasks  -  is  already  complete. 

Unfortunately,  an  analysis  using  the  STL  model  of  the  Manduca  Sexta  forewing  has  one 
major  downfall;  meshing.  Figure  31  shows  a  portion  of  the  hawkmoth  wing  root.  The  second 
and  third  veins  (from  the  leading  edge)  are  shown,  oriented  with  the  leading  edge  to  the  right^^. 
Inspecting  only  a  small  portion  of  the  wing  reveals  it  is  unlikely  that  even  a  highly  refined  mesh, 
consisting  only  of  tetrahedral  elements,  would  successfully  discretize  this  structure.  Moreover,  if 
such  a  mesh  were  successful,  it  would  hardly  be  desirable  due  to  the  excess  bending  stiffness 
associated  with  tetrahedrons  -  especially  in  thin  structures'^  [1,2,  44]. 

As  a  n  a  Iternative  to  a  ttempting  a  d  irect  me  sh  o  f  th  e  S  TL  mo  del,  o  ne  could  a  ttempt  a 
geometry  “cl  ean-up”  p  rocess.  A  reasonable  geometric  s  implification  could  be  achieved  with 
smooth  surfaces  instead  of  triangular  facets.  Likewise,  one  could  enforce  a  discrete  separation  of 
veins  and  membrane  instead  of  the  nebulous  transitions  shown  in  Figure  31.  Still,  this  approach 
is  undesirable  due  to  the  induced  stiffness  of  a  tetrahedral  mesh,  and  the  time  (possibly  weeks) 

A  planform  view  of  this  image  is  shown  in  Appendix  B  -  CT  Images. 

This  statement  is  regarding  volumetric  and  shear  locking  in  linear  tetrahedral  elements.  Q  uadratic  tetrahedrons 
offer  more  accuracy,  but  the  overall  complexity  would  force  the  model  to  have  large  quantities  (possibly  millions) 
of  unknown  degrees  of  freedom. 
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required  to  aehieve  an  adequately  simple  geometric  representation.  C  onversely,  the  thinness  of 
the  wing  membrane  and  tapering  of  the  generally  circular  veins  suggests  it  is  possible  to  achieve 
a  reasonably  accurate  solution  using  shell  and  beam  elements,  respectively.  T  hus,  the  selected 
approach  was  to  create  a  CAD  model  that  reasonably  approximates  the  Manduca  Sexta  forewing 
using  measurements  taken  from  the  CT  images. 


Figure  3 1 :  Stereolithography  (STL)  image  of  a  Manduca  Sexta  wing  root.  Shown 
are  the  second  and  third  veins  from  the  leading  edge  (leading  edge  is  to 
the  image  right). 


3.2.  Element  Selection 

The  Manduca  Sexta  wing  is  divided  into  two  major  components:  ve  in  and  membrane. 
The  ve  ins  a  re  generally  or  iented  i  n  t  he  s  panwise  di  rection,  a  nd  t  aper  f  rom  r  oot  t  o  t  ip.  T  he 
membrane  has  no  obvious  orientation  (from  CT  scans),  and  has  a  generally  constant  thickness  of 

1  ”7 

approximately  0.12  millimeters.  The  relative  thinness  of  the  wing  membrane  requires  one  to 
consider  what  element  type  is  appropriate  for  modeling.  The  veins,  however,  can  be  modeled  as 

The  term  “membrane”  is  used  here  in  the  tissue  sense  of  the  term,  not  as  an  element  type. 
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wire  b  earns,  with  appropriately  as  signed  cross-seetions.  S  till,  n  umerieal  studies  m  ust  be 
conducted  to  determine  to  what  extent  the  wing  frequencies  are  affected  by  the  vein  taper.  In 
other  words,  a  total  of  16  c  ross-sectional  measurements  were  taken  from  the  C  T  images;  how 
many  of  these  cross-sectional  assignments  are  required  to  adequately  model  the  behavior  of  the 
veins  as  wire  beams?  The  studies  that  follow  answer  these  questions. 

3.2.1.  Modeling  with  Membrane  Elements 

Physical  handling  of  the  Manduca  Sexta  wing  reveals  a  structure  that  is  highly  dominated 
by  vein  stiffness,  with  the  membrane  fulfilling  a  displacement  continuity  role.  It  appears  that  the 
“membrane”  -  i.e.,  t  he  t  issue  be  tween  t  he  v  eins  -  may  be  s  tructurally  eh  aracterized  b  y 
membrane  elements,  in  the  finite  element  sense  of  the  term.  By  convention,  membrane  elements 
are  sheets  that  exist  in  plane-stress  [1].  Any  stresses  are  carried  parallel  to  the  midsection  of  the 
element,  with  all  others  assumed  zero.  Membrane  elements  have  no  bending  or  transverse  shear 
stiffness.  As  a  result,  any  strains  normal  to  the  element  or  in  the  transverse  direction  (As,)  the 
result  of  Poisson  effects  in  the  constitutive  relations,  when  the  element  lies  in  the  plane  described 
by  Du  and  022  [44].  A  schematic  of  a  four-noded  membrane  element  and  its  positive  normal 
vector,  n,  are  shown  below,  in  Figure  32. 
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Before  attempting  to  create  a  model  that  represents  the  entirety  of  the  hawkmoth  wing,  a 
model  was  constructed  to  examine  the  behavior  of  the  vein-membrane  interaction.  T  hree  wire 
beams  of  circular  cross-section  are  connected  by  two  planar  sections  of  membrane  elements,  as 
shown  in  Figure  33.  The  three  beams  and  two  membrane  sections  were  tied  together,  forcing 
their  displacements  to  be  the  same  at  the  node  locations,  and  were  of  dimensions  representative 
of  the  Manduca  Sexta  wing.  Two  test  cases  were  run;  the  first  with  linear  membrane  elements, 
and  the  second  with  quadratic  elements.  For  simplicity,  aluminum  was  selected  as  the  material 
in  both  cases. 


bottom)  connected  by  two  membrane  sections. 


The  f  irst  f  requency  an  alysis  was  c  onducted  us  ing  two-noded  linear  b  earns  ( B3 1)  and 
four-noded  qu  adrilateral  membrane  el  ements  (M3D4).  Mesh  d  ensities  b  etween  3  04  el  ements 
(348  node  s)  a  nd  750,000  e  lements  (  756,504  no  des)  w  ere  c  onsidered  f  or  c  onvergence.  The 
model  is  planar  with  only  one  element  through  the  thickness,  thus  a  full-integration  method  was 
used  to  avoid  any  idiosyncrasies  associated  with  constant  strain,  reduced  integration  formulations 
[7].  The  first  m  odeshape  i  s  shown  below  in  Figure  34.  Most  not  ably,  t  he  c  entral  b  earn 
(analogous  to  a  vein)  is  in  motion,  affecting  only  its  adjacent  membranes.  There  is  transference 
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of  motion  among  the  membrane  elements  in  the  out-of-plane  direction;  in  its  present  state,  this 
beam-membrane  configuration  is  incapable  of  describing  the  eigenvectors  seen  experimentally  in 
the  Manduca  Sexta  wing. 


Figure  34:  Finite  element  results  of  three  quadratic  wire  beams  connected  by  two 


linear  membrane  sections.  Color  bands  denote  displacement. 

The  large  vertical  displacement  of  the  beam  relative  to  the  undeformed  state  suggests  a 
more  refined  mesh  may  alleviate  the  discontinuity  of  slope  near  the  central  beam.  To  this  end,  a 
convergence  s  tudy  w  as  co  nducted  tod  etermine  h  ow  m  any  1  inear  elements  ar  e  r  equired  t  o 
sufficiently  conduct  an  accurate  frequency  analysis.  The  results  of  this  convergence  study  are 
shown  below,  in  Figure  35.  Although  the  frequency  values  vary  with  mesh  size,  the  modeshapes 
remain  similar,  e  ven  w  ith  a  s  m  any  a  s  750,000  membrane  e  lements  with  756,504  nodes  (by 
comparison.  Figure  34  shows  1,950  e  lements).  In  f  act,  t  he  onl  y  change  i  n  t  he  e  igenvalue 
structure  is  the  thickness  of  the  membrane  “peak”  near  the  beam.  As  the  membrane  element  size 
approaches  an  infinitesimal  value,  the  width  of  the  peak  approaches  the  thickness  of  the  beam. 
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Membrane  Elements 

Figure  35:  Convergence  study  of  a  finite  element  model  with  three  quadratic  wire 
beams  connected  by  two  linear  membrane  sections. 


It  is  a  pparent  f  rom  e  xamining  th  e  results  o  f  th  e  above  convergence  s  tudy  t  hat  linear 
membrane  elements  are  insufficient  for  modeling  the  Manduca  Sexta  wing.  Even  in  a  simple 
case  involving  three  wire  beams  and  two  membrane  sections,  nearly  one  million  elements  were 
required  be  fore  a  be  havior  r  esembling  c  onvergence  w  as  a  chieved.  F  urthermore,  uni  ike  t  he 
behavior  obs  erved  e  xperimentally,  th  e  lin  ear  m  embrane  e  lements  d  o  n  ot  a  How  c  ontinuity  o  f 
displacement  ( or  s  lope)  normal  t  o  t  he  e  lement  pi  ane.  Either  a  q  uadratic  m  embrane,  o  r  a  n  ew 
element  entirely,  is  required. 

A  second  analysis  was  conducted  using  the  model  shown  in  Figure  33,  this  time  with 
eight-noded  quadratic,  quadrilateral  membrane  (M3D8)  elements.  The  M3D8  element  resembles 
the  linear  membrane  of  Figure  32,  with  added  midside  nodes;  like  the  linear  membrane  element, 
each  node  has  only  the  two  in-plane  degrees  of  freedom.  The  first  modeshape  is  shown  below, 
in  Figure  36.  In  t  his  cas  e,  t  he  r  esults  a  re  s  omewhat  s  imilar  to  th  ose  obtained  u  sing  lin  ear 
membrane  e  lements,  w  ith  tw  o  d  istinctions.  F  irst,  the  a  dded  m  idside  n  odes  i  n  t  he  qua  dratic 
element  pr  oduce  i  ntra-element  v  ibrations  a  t  th  e  n  atural  f  requencies  in  the  v  ertical  d  irection. 
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Second,  in  keeping  with  their  second-order  shape  functions,  the  quadratic  membrane  elements 
are  able  to  affect  the  displacement  of  adjacent  membranes;  this  advantage,  however,  should  not 
be  ove  rstated  s  ince  t  he  m  odel  r  esults  s  how  a  ny  di  splacement  i  nfluence  i  s  ne  gligible 
approximately  three  elements  from  the  central  beam. 


Figure  36:  Finite  element  results  of  three  quadratic,  wire  beams  connected  by  two 
quadratic  membrane  sections.  Color  bands  denote  displacement. 


Figure  37  shows  the  same  model  as  in  Figure  36,  but  with  a  more  discretized  mesh.  The 
effects  of  mesh  refinement  using  quadratic  membrane  elements  are  similar  to  those  of  the  linear 
elements.  Decreased  el  ement  s  ize,  co  rresponding  t  o  an  i  ncrease  i  n  t  he  o  verall  n  umber  o  f 
elements,  r  esults  in  an  exaggerated  m  odeshape  t  hat  r  esembles  a  D  irac  d  elta  f  unction  at  th  e 
middle  beam.  Again,  the  effects  of  the  beam  displacement  are  unfelt  at  a  distance  of  only  three 
membrane  elements  from  the  beam. 


Figure  37:  Finite  element  results  of  three  quadratic,  wire  beams  connected  by  two  quadratic 
membrane  sections,  with  a  refined  mesh.  Color  bands  denote  displacement. 
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Repeating  t  he  approach  u  sed  w  ith  t  he  1  inear  m  embrane  el  ements,  a  co  nvergence  s  tudy 
was  conducted  to  determine  whether  the  mesh  produces  a  stable  result.  As  before,  if  the  model 
results  are  independent  of  mesh  size,  one  may  deduce  that  increasing  the  number  of  elements 
will  ha  ve  1  ittle  (  or  no)  e  ffect  on  t  he  be  havior  of  t  he  s  tructure.  Thus,  i  f  t  he  pr  oblems  of 
displacement  only  near  the  central  beam  persist,  even  though  the  frequency  analysis  is  stable,  it 
becomes  evident  that  membrane  elements  are  ill-suited  for  modeling  the  Manduca  Sexta  wing. 
The  results  of  the  convergence  study  involving  quadratic  membranes  are  shown  below,  in  Figure 
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Membrane  Elements 

Figure  38:  Convergence  study  of  a  finite  element  model  with  three  quadratic  wire 
beams  connected  by  two  quadratic  membrane  sections. 

Highlighting  the  computational  usefulness  of  quadratic  formulations,  this  model  achieved 
convergence  after  only  403  elements  (1109  nodes),  at  ap proximately  3  Hertz  below  the  linear 
membrane  s  olution.  Despite  t  he  us  efulness  of  t  his  w  ell-known  r  esult,  t  he  b  ehavior  of  t  he 
converged  model  is  unchanged  from  that  of  Figure  37.  Membrane  elements  -  linear  or  quadratic 
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-  are  not  able  to  deseribe  the  motion  observed  in  Manduca  Sexta  wings,  but  why?  T  o  answer 
this  question,  one  must  revisit  the  first  prineipals  of  membrane  element  formulation. 

3.2.2.  Membrane  Element  Formulation 

Membrane  elements  are  plane-stress  formulations,  that  is,  the  only  non-zero  stresses  lie 
in-plane  with  the  midsection  of  the  element.  Figure  39  shows  a  membrane  element  described  by 
the  Xii  —  X22  plane,  in  which  all  stresses  in  the  5z-direction  are  assumed  to  be  zero  (z  =  1,  2  3). 
There  is  a  vast  difference  between  a  structure  in  plane-stress  and  a  structure  assumed  to  be  in 
plane-stress,  though  the  difference  wording  is  minor.  Structures  in  plane-stress  are  sufficiently 
thin,  s  uch  t  hat  onl  y  s  tress  tractions  nor  mal  t  0  t  he  m  idplane  a  re  allowed  t  0  be  nonz  ero. 
Computationally,  elements  assumed  to  be  in  plane  stress  leave  components  in  the  5z-positions  of 
stress  matrix  uncalculated  and  filled  with  zeros  [44].  In  ot her  w ords,  a ny  out-of-plane  stress 
present  in  a  me  mbrane  element  w  ill  n  ot  b  e  c  alculated;  if  stress  v  alues  a  re  i  ncluded  as  in  itial 
conditions,  they  will  be  overwritten  by  zeros. 


^33 


Figure  39:  Membrane  element  described  by  the  Xn-X22  plane. 

To  illustrate  this  point,  a  model  is  shown  below  in  Figure  40,  consisting  of  a  rectangular 
part  w  ith  a  ve  rtical  1  oad  i  n  t  he  c  enter.  T  he  pa  rt  i  s  bou  nded  on  a  11  f  our  s  ides  by  an  en  castre 
condition,  and  meshed  with  quadratic  membrane  elements.  For  convenience,  aluminum  material 
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properties  are  used,  and  the  load  in  the  center  is  a  concentrated  point-load  of  1,000  1  bf  (4.448 
kN).  The  results  of  this  analysis  are  given  below,  in  Figure  41 . 


Figure  40:  Rectangular  structure  comprised  of  membrane  elements,  with  a  central  load. 


The  r  esults  mig  ht  a  stonish  o  ne  u  nfamiliar  w  ith  th  e  u  nderlying  ma  thematical  b  asis  o  f 
membrane  elements:  a  vertical  load  of  1,000  1  bf  results  in  zero  von  M  ises  stress  at  all  points. 
Clearly  these  results  are  erroneous,  but  only  because  the  initial  conditions  violate  the  plane-stress 
scenario.  The  analysis  results  in  a  zero  stress,  zero  displacement  condition,  because  the  solver 
set  all  loading  terms  normal  to  the  midsection  equal  to  zero  -  eliminating  the  1,000  Ibf  load. 


Figure  41 :  Results  of  a  vertical  point  loads  using  only  membrane 
elements.  Blue  denotes  zero  von  Mises  stress. 


Returning  t  o  t  he  f  requency  a  nalysis  of  t  he  above  beam-membrane  m  odel  with  th  is 
understanding  helps  clarify  why  such  odd  be  havior  is  observed.  Consider  the  central  beam,  its 
adjacent  linear  membrane  element,  and  another  membrane  element  adjacent  to  the  first,  as  shown 
in  Figure  42.  In  this  scenario,  the  beam  is  deflected  vertically,  as  in  Figure  36.  Due  to  the  tie 
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constraint,  membrane  one  is  foreed  to  accept  a  rigid-body  rotation  and  displace  vertically  at  the 
right  e  nd,  m  aintaining  continuity  o  f  di  splacement  w  ith  t  he  be  am.  By  definition,  m  embrane 
elements  a  re  o  nly  a  ble  to  tr  ansmit  in  -plane  d  isplacement;  t  hey  can  not  accep  t  out  -of  p  lane 
(transverse)  strain  or  b  ending  t  erms.  Due  t  o  t  his  f  ormulation,  m  embrane  t  wo  i  n  Figure  42 
maintains  continuity  of  in-plane  (horizontal)  displaeement  with  membrane  one,  without  rotation 
(due  1 0  1  ack  oft  ransverse  s  hear).  T  he  r  esult  of  t  hese  combined  e  ffects  i  s  t  he  di  scontinuous 
“kink”  seen  in  Figure  36.  As  the  user  refines  the  mesh,  the  angle  between  membrane  elements 
one  a  nd  t  wo  i  nereases,  a  pproaching  t  he  1  imiting  e  ondition  of  e  lement  one  be  ing  nor  mal  t  o 
element  two.  T  his  results  in  the  lack  of  convergence  seen  using  linear  membrane  elements  in 
Figure  35. 


Figure  42:  Schematic  o  f  a  d  eflected  w  ire  b  earn 
connected  to  a  linear  membrane  element, 
which  is  connected  to  another  linear 
membrane  element. 


Although  the  above  commentary  describes  the  linear  membrane  elements,  a  similar  effect 
is  s  een  i  n  t  he  q  uadratic  m  embrane  elements.  The  ke  y  d  ifference  i  n  t  heir  r  esponses  is  their 
respective  shape  functions.  As  their  names  imply,  linear  serendipity  elements  have  one  node  at 
each  c  omer,  and  di  splace  i  n  a  1  inear  fashion  along  t  heir  edges,  ass  hown  i  n  Figure  34. 
Conversely,  qua  dratie  s  erendipity  elements  ha  ve  m  idside  node  s  be  tween  t  he  c  orners,  a  nd 
displace  according  to  a  quadratic  interpolation  function.  These  interpolation  functions  are  given 
in  Ref  [44]. 
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Another  notable  difference  in  the  responses  of  the  linear  and  quadratic  elements  is  the 
displacement  o  f  m  embranes  that  ar  e  not  ad j  acent  t  o  t  he  w  ire  b  earn.  Figure  43  shows  s  imilar 
results  as  Figure  37,  but  with  fewer  quadratic  membrane  elements,  to  aid  visibility.  Due  to  the 
requirement  of  continuity  of  displacement  at  the  nodes,  the  membranes  adjacent  to  the  center 
beam  cause  elements  not  connected  to  the  beam  to  displace  vertically.  Furthermore,  since  the 
membrane  elements  have  no  bending  stiffness,  the  slope  is  discontinuous  at  the  nodes. 


Figure  43:  First  modeshape  of  a  three  beam  stmeture  separated  by  quadratie  membrane 
elements.  Color  bands  denote  J5-displaeement. 


Despite  their  higher-order  interpolation  functions,  the  quadratic  elements  suffer  from  the 
same  pr  oblem  s  een  i  n  t  he  1  inear  m  embrane  e  lements:  t  he  di  scontinuity  of  s  lope  at  t  he  node  s 
results  in  a  lo  ss  o  f  tr  ansmitted  s  tress.  Thus,  th  e  majority  o  f  th  e  me  mbrane  s  ection  remains 
unaffected  by  the  vibration  of  the  central  beam.  From  these  studies,  it  is  clear  that  a  thin  element 
capable  o  f  c  arrying  b  ending  an  d  s  hear  i  s  r  equired  tod  escribe  t  he  “m  embrane”  t  issue  o  f  t  he 
Manduca  Sexta  wing. 


3.2.3.  Modeling  with  Shell  Elements 

In  1  ight  of  t  he  de  ficiencies  of  t  he  m  embrane  e  lement  f  or  t  he  purpose  of  m  odeling  a 
Manduca  Sexta  wing,  a  s  hell  el  ement  i  s  s  elected  f  or  t  he  n  ext  s  cries  o  f  an  alyses.  There  ar  e 
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essentially  three  types  of  shell  elements  used  in  Abaqus.  The  first  type  is  a  classical  Kirchhoff 
shell,  in  which  bending  is  restricted  so  that  vectors  normal  to  the  midpoint  of  the  shell  surface 
remain  normal  after  deformation.  These  elements  are  appropriate  for  sufficiently  thin  structures 
[44],  T  he  s  econd  t  ype  of  s  hell  f  alls  i  n  t  he  M  indlin  pi  ate  f  amily,  i  n  which  be  nding  i  s  not 
restricted  b  y  t  he  nor  mal  ve  ctor  f  ound  i  n  Kirchhoff  elements.  S  hell  e  lements  o  f  th  e  M  indlin 
formulation  are  appropriate  for  thick  plate  scenarios  [14].  The  third  type  of  shell  element  in  the 
Abaqus  library  is  the  Ko iter- Sanders  shell.  The  Koiter-Sanders  shell  is  a  general-purpose  shell, 
in  w  hich  t  hickness  i  s  a  llowed  t  o  va  ry  dur  ing  deformation  due  t  o  P  oisson  e  ffects  [1].  The 
general-purpose  shell  is  a  four-noded  serendipity  element  resembling  the  membrane  of  Figure 
32,  but  with  all  six  degrees  of  freedom  active.  The  Abaqus  Theory  Manual  states  the  following 
regarding  the  Koiter-Sanders  (“general-purpose”)  shell  elements  [1]; 

The  general-purpose  elements  provide  robust  and  accurate  solutions  in  all 
loading  conditions  for  thin  and  thick  shell  problems.  Thickness  change  as  a 
function  of  in-plane  deformation  is  allowed  in  their  formulation.  They  do  not 
suffer  from  transverse  shear  locking,  nor  do  they  have  any  unconstrained 
hourglass  modes.... No  hourglass  control  is  required  for  the  axisymmetric 
general-purpose  shells,  nor  in  the  bending  and  membrane  response  of  the  fully 
integrated  element  S4.  The  membrane  kinematics  of  S4  are  based  on  an  assumed- 
strain  formulation  that  provides  accurate  solutions  for  in-plane  bending 
behavior  ....These  elements  use  simplified  methods  for  strain  calculation  and 
hourglass  control  and  offer  significant  advantages  in  computational  speed. 

With  this  knowledge,  the  previous  exercise  in  which  three  beams  were  separated  by  two 

membrane  s  ections  w  as  r  epeated.  T  his  t  ime,  t  he  m  embrane  s  ections  w  ere  r  eplaced  b  y  shell 

sections  of  the  recommended  Koiter-Sanders  formulation.  In  light  of  the  advantage  offered  by 

general-purpose  elements,  only  four-noded,  doubly-curved  shells  with  reduced  integration  and 

hourglass  control  (S4R)  were  used;  other  shell  types  were  excluded  from  this  analysis.  A  total  of 

950  elements  (1,020  nodes)  were  used  in  this  example. 
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The  f  irst  s  ix  modeshapes  of  t  his  s  tructure  resembled  pi  ate-like  m  odes,  a  s  t  hough  t  he 
beams  were  not  present.  T  he  seventh  modeshape,  however,  was  of  a  configuration  useful  for 
comparison  to  the  first  modeshape  of  the  beam-membrane  analysis  (shown  below,  in  Figure  44). 
In  t  his  c  ase,  t  he  c  entral  be  am  i  s  flapping  out  of  pha  se  w  ith  t  he  t  wo  e  nd  be  ams,  but  t  he  s  hell 
elements  a  re  not  e  xhibiting  t  he  di  scontinuities  of  di  splacement  a  nd  s  lope  obs  erved  i  n  t  he 
membrane  elements. 


Figure  44:  Seventh  modeshape  of  a  three  beam  structure  separated  by  general-purpose  shell  elements. 

The  results  shown  in  Figure  44  qualitatively  resemble  how  one  might  expect  the  beam- 
shell  structure  to  deform.  Still,  simply  passing  the  first-glance  “sanity-check”  is  an  insufficient 
means  of  de termining  t he  a ppropriateness  of  the  general-purpose,  K oiter-Sanders  el  ement  for 
modeling  a  Manduca  Sexta  forewing.  With  no  physical  beam-shell  structure  against  which  one 
might  c  ompare  the  finite  e  lement  results,  th  e  e  fficacy  o  f  th  ese  s  hells  remains  in  q  uestion. 
Fortunately  there  is  another  applicable  structure  with  known  frequency  data  and  modeshapes. 

Recall  the  results  of  the  paper  wing  frequency  tests,  given  in  Chapter  I  of  this  document. 
A  fl  at  (n  on-cambered),  paper  w  ing  o  f  s  imilar  g  eometry  t  o  t  he  Manduca  Sexta  forewing  wa  s 
scanned  using  a  laser  vibrometer.  From  these  scans,  natural  frequencies  and  modeshapes  were 
obtained  f  or  c  omparison  t  o  t  he  a  ctual  Manduca  Sexta  wing.  While  th  ese  m  odeshapes  p  roved 
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different  than  those  of  the  hawkmoth  wing  at  the  third  and  higher  modes,  they  nevertheless  serve 
as  an  adequate  basis  for  finite  element  model  comparison. 

3.2.4.  Modeling  a  Paper  Wing 

Built  into  Abaqus  is  a  feature  in  which  the  user  can  import  an  image,  and  trace  over  it  to 
form  sketches.  These  sketches  can  later  be  turned  into  parts.  Using  this  method,  a  model  of  the 
paper  wing  was  constructed  with  the  same  dimensions  as  the  physical,  paper  wing.  Figure  45 
shows  the  completed  outline  of  the  paper  wing.  The  wing  used  during  modeling  is  truncated  at 
the  root  so  that  the  clamped  portion  of  the  wing  (near  the  root)  is  not  modeled.  Instead,  a  fixed 
boundary  condition  is  enforced  at  the  left  edge  of  the  wing,  denoted  by  the  “V”  in  Figure  45  . 


Figure  45:  Geometry  tracing  of  paper  wing  used  in  frequency  experiments. 


The  w  ing  m  odel  w  as  m  eshed  us  ing  t  he  four-noded,  doubl  y-curved  s  hells  w  ith  r  educed 
integration  and  hourglass  control  (S4R)  -  the  same  Koiter-Sanders  shell  elements  used  in  the 
previous  analysis.  A  converged  mesh  of  929  of  these  general-purpose  elements  (1,016  nodes) 
was  used.  Material  properties  for  the  wing  (made  of  copier  paper)  were  obtained  by  calculating 


For  comparison,  see  Figure  25. 
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the  de  nsity  ba  sed  on  k  nown  w  eight  a  nd  vol  ume  m  easurements,  e  stimating  P  oisson’s  r  atio 
(assumed  a  value  of  0.3),  and  ehanging  Young’s  modulus  until  reasonable  frequeney  values  were 
obtained^^.  The  results  of  t  his  a  nalysis  a  re  given  below,  i  n  Table  9.  As  i  s  e  vident  f  rom 
examining  t  he  va  lues,  e  ven  a  poor  ly  t  uned  m  odel  ba  sed  on  e  stimated  m  aterial  pr  operties 
produees  reasonably  aeeurate  results. 

Table  9.  C  omparison  of  natural  frequencies  of  a  paper 
wing  i  n  v  acuum,  a  nd  a  n  F  E  s  olution  us  ing 
general-purpose  shell  elements  (S4R). 

Experimental  Wing,  Hz  FE  Result,  Hz 
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22 

0)2 

96 

109 

"3 

162 

159 

U)4 

472 

466 

The  m  odeshapes  o  ft  he  f  inite  el  ement  m  odel  ar  e  al  so  r  epresentative  oft  hose 
experimentally  observed  in  paper  wing.  Reeall  from  Chapter  I,  the  first  four  modeshapes  of  the 
paper  w  ing  ar  e,  i  n  order,  f  irst  be  nding  ( or  flapping),  t  orsion,  s  eeond  be  nding,  a  nd  f  inally  t  he 
“saddle”  mode  with  tip-defleetion.  The  results  of  the  finite  element  model  of  the  paper  wing  are 
given  below.  Figure  46  shows  the  first  modeshape,  a  bending  mode;  Figure  47  shows  a  torsional 
mode.  Figure  48  gives  the  eharaeteristie  seeond-bending  mode,  and  the  modeshape  shown  Figure 
49  resembles  the  saddle-mode  with  tip  displaeement  noted  in  the  paper  wing. 

The  r  esults  of  t  he  finite  e  lement  model  -  both  t  he  f  requeneies  a  nd  m  odeshapes 
(displaeements)  -  appear  to  suffieiently  deseribe  the  paper  wing  tested  in  the  vaeuum  ehamber 
(modeshapes  s  hown  Chapter  I ,  pa ges  39  through  41).  F  urthermore,  t he  r  esults  of  the  modal 
analysis  support  the  statement  made  in  the  Abaqus  Theory  Manual,  in  whieh  the  Koiter-Sanders 
elements  were  praised  for  their  aeeuraey  and  robustness. 

Material  properties  used  where:  p  =  8.6x10'^  g/mm^,  v  =  0.3,  E  =  7.5  MPa 
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Figure  46:  First  modeshape  (22  FIz)  of  a  flat,  FE  wing,  shown  from  an 
isometric  view  (a)  and  a  tip  view  (b). 


(a)  (b) 

Figure  47:  Second  modeshape  (109  Flz)  of  a  flat,  FE  wing,  shown 
from  an  isometric  view  (a)  and  a  tip  view  (b). 


(a)  (b) 

Figure  48:  Third  modeshape  (159  Flz)  of  a  flat,  FE  wing,  shown  from 
an  isometric  view  (a)  and  a  tip  view  (b). 


(a)  (b) 

Figure  49:  Fourth  modeshape  (466  Flz)  of  a  flat,  FE  wing,  shown 
from  an  isometric  view  (a)  and  a  tip  view  (b). 
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3.3.  Camber  Effects 


Previous  sections  focused  on  the  effects  of  various  element  types  when  modeling  a  flat 
wing.  Based  on  the  results,  the  quadratic  beam  and  Koiter-Sanders  shell  were  obvious  choices. 
Despite  promising  results  when  comparing  the  finite  element  wing  to  the  paper  wing,  neither 
provides  a  m  odal  r  esponse  similar  to  th  e  a  ctual  w  ing.  Indeed,  at  t  his  s  tage  i  n  t  he  s  tructural 
dynamic  a  nalysis  of  t  he  Manduca  Sexta  wing,  it  is  u  nclear  to  w  hat  extent  the  r  esponse  i  s 
dominated  by  camber  or  the  presence  of  veins.  To  simplify  the  process  of  sorting  camber  and 
vein  stiffness  effects,  the  complexity  of  the  finite  element  model  will  be  incrementally  increased. 
Now,  we  focus  on  camber. 

Using  an  approach  similar  to  that  of  Section  3.2.4,  a  model  of  a  homogeneous,  isotropic, 
cambered  wing  was  constructed.  However,  in  order  to  introduce  camber  without  pre-stressing 
the  wing,  a  slightly  different  approach  to  model  creation  was  taken.  This  time,  a  curved  surface 
of  constant  radius  was  created  in  Abaqus.  Using  the  same  photo  sketching  technique  as  before, 
the  planform  of  the  Manduca  Sexta  forewing  was  sketched  onto  a  plane  above  the  curved  shell 
and  projected  downward,  as  shown  in  Figure  50.  Finally,  the  material  outside  wing  planform 
sketched  on  the  curved  surface  was  removed,  leaving  the  cambered  wing  of  Figure  51. 

One  of  t  he  di  sadvantages  of  using  a  planform  vi  ew  of  a  cambered  w  ing  to  t  race  the 
venation  structure  is  the  loss  of  wing  topography  data.  Despite  any  camber  in  the  wing,  a  p  Ian 
view  of  the  venation  structure  is  planar  {i.e.,  does  not  include  three-dimensional  topography). 
Thus,  since  the  Manduca  Sexta  wing  is  a  three-dimensional  structure,  a  two-dimensional  image 
simply  cannot  adequately  describe  the  wing  curvature.  Creating  a  curved  shell  and  sketching  the 
venation  pa  item  f  rom  t  he  C  T  i  mages  as  s  hown  i  n  Figure  50  helps  to  mitig  ate  th  is  lo  ss  in 
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topography  data.  Essentially,  the  goal  of  this  method  is  to  use  a  two-dimensional  photo  and  an 
average  value  of  eamber  to  restore  as  mueh  of  the  original  wing  geometry  as  possible. 

The  radius  of  eurvature  of  the  shell  onto  whieh  the  planform  of  the  wing  was  projeeted  is 
eonstant.  In  other  words,  all  parts  of  the  wing  undergo  the  same  camber.  Conversely,  the  actual 
Manduca  Sexta  wing  i  s  not  of  c  onstant  c  urvature,  w  ith  c  amber  be  ing  do  minant  ne  ar  t  he  r  oot. 
The  camber  value  used  in  the  model  was  achieved  by  examining  images  of  the  cross-sections 
from  the  CT  scans.  An  arithmetic  mean  of  the  angle  created  between  the  leading  edge  and  the 
wing  centroid  was  taken  and  used  to  construct  a  three-point  arc.  This  arc  was  extruded  and 
became  the  shell  onto  which  the  wing  planform  was  sketched. 


Figure  50:  Projected  Manduca  Sexta  wing  pi  anform  us  ed  dur  ing 


cambered  part  cr  cation.  The  p  art  w  as  sketched  o  n  t  he 
plane  above  the  curved  surface. 


This  value  was  determined  to  be  15°,  or  13%  of  the  chord.  An  explanation  of  wing  camber  is  found  in  Appendix 
D  -  Wing  Camber. 
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(a)  (b) 

Figure  5 1 :  Cambered  Manduca  Sexta  wing  model.  S  hown  a  re  t  he 
isometric  view  (a)  and  tip  view  (b) 


The  w  ing  s  hown  i  n  Figure  5 1  was  m  eshed  w  ith  t  he  four-noded,  doubl  y-curved  s  hells 
(S4R)  used  previously.  Nine  hundred  elements  covered  the  cambered  wing  during  the  frequency 
analysis.  In  order  to  make  a  direct  comparison  to  the  flat  wing  model,  the  material  properties  of 
paper  were  used  during  the  cambered  analysis.  For  comparison,  the  frequency  results  for  the  flat 
and  cambered  cases  are  given  below,  in  Table  10.  By  examining  the  results  of  both  cases,  one 
immediately  finds  that  the  curvature  of  the  wing  dramatically  increases  its  stiffness. 

Table  10.  Comparison  of  natural  frequencies  of  flat 
and  c  ambered  w  ing  FE  solutions,  using 
general-purpose  shell  elements  (S4R). 

Flat  Wing,  Hz  Cambered  Wing,  Hz 


22 

159 

(x)2 

109 

173 

0)3 

159 

797 

0)4 

466 

1057 

The  modeshapes  of  the  cambered  wing  analysis  are  given  in  the  figures  below.  The  first 

modeshape,  shown  in  Figure  52,  differs  slightly  from  that  of  the  flat  wing.  Although  the  wing  is 

still  displacing  in  a  f  ashion  r  eminiscent  o  f  classical  f  irst  b  ending,  a  la  rge  p  ortion  o  f  th  e 

movement  is  now  seen  in  along  the  trailing  edge  instead  of  at  the  tip.  T  he  second  modeshape, 

however,  appears  to  b  e  similar  in  ma  ny  respects  to  th  e  d  eformations  o  f  the  fiat  w  ing.  H  ere, 
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torsion  i  s  the  pr  imary  d  river,  t  hough  t  he  a  xis  o  f  r  otation  a  Iso  s  eems  t  o  ha  ve  s  hifted  s  lightly 
towards  the  trailing  edge. 


(a)  (b) 


Figure  52:  First  modeshape  (159  FIz)  of  a  cambered,  FE  wing,  shown 
from  an  isometric  view  (a)  and  a  tip  view  (b). 


(a)  (b) 

Figure  53:  Second  modeshape  (173  Hz)  of  a  cambered,  FE  wing, 
shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 

The  third  modeshape,  shown  below  in  Figure  54,  is  by  far  the  most  striking  evidence  of 
the  ef  feet  o  f  c  ambering  t  he  w  ing.  Here,  t  he  t  bird  m  ode  i  s  of  a  s  addle  s  hape,  w  ith  s  imilar 
deflections  t  o  t  he  t  bird  m  ode  of  t  he  a  ctual  Manduca  Sexta  forewing.  T  his  i  s  a  n  obvi  ous 
departure  from  the  flat  wing  analyzed  previously,  in  which  an  extra  second-bending  mode  was 
observed,  pus  hing  t  he  s  addle  m  ode  t  o  t  he  f  ourth  pos  ition.  The  ap  pearance  o  f  t  he  s  addle 
modeshape  (corresponding  to  the  third  mode)  indicates  that  camber  is  an  important  feature  in  the 
frequency  response  of  the  Manduca  Sexta  wing.  Moreover,  the  large  midside  deflections  at  the 
leading  and  trailing  edges  of  the  homogeneous  shell  model  highlight  the  structural  importance  of 
the  stiffness  imparted  by  the  veins.  By  contrast,  the  hawkmoth  wing  shows  little  midside  action 
in  the  third  eigenvector,  likely  due  to  the  added  vein  stiffness. 
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(a)  (b) 

Figure  54:  Third  modeshape  (797  Hz)  of  a  cambered,  FE  wing,  shown 
from  an  isometric  view  (a)  and  a  tip  view  (b). 


The  p  receding  s  tudy  o  f  camber  effects  o  n  n  atural  f  requency  was  r  epeated  f  or  s  everal 
values  between  zero  and  15°  (13%  of  chord)^'.  The  same  homogeneous,  isotropic  wing  meshed 
with  general-purpose  shells  w as  u sed.  T  he  results  are  plotted  below,  in  Figure  55.  The  first 
natural  f  requency  i  ncreases  1  inearly  with  c  amber,  w  hile  t  he  s  econd  m  ode,  a  Iso  i  ncreasing,  i  s 
substantially  less  affected  by  a  change  in  camber.  The  third  mode  has  the  unique  distinction  of 
changing  modeshapes  during  this  analysis.  At  low  camber  values  (less  than  4.5%  chord),  the 
third  m  odeshape  i  s  t  he  second  be  nding  m  ode  o  bserved  pr  eviously  i  n  flat  w  ings.  Above  t  he 
cutoff  value  of  4.5%  of  chord,  the  camber-induced  stiffness  eliminates  the  second  bending  mode. 
Here,  t  he  s  addle  m  ode  experiences  a  n  egative  frequency  s  hift,  an  d  becomes  t  he  third  m  ode 
observed  in  the  Manduca  Sexta  wing.  It  is  notable  that  the  second  bending  and  saddle  modes  do 
not  experience  a  mode  reversal  at  high  camber  values.  Instead,  the  second  bending  mode  simply 
disappears  -  or  changes  shape  so  dramatically  that  it  can  no  longer  be  characterized  as  “second 
bending.”  Regardless  of  the  modeshape  presented,  the  third  mode  is  strongly  affected  by  even 
slight  changes  in  camber. 


The  cam  her  i  n  t  he  Manduca  Sexta  wing  va  ries  f  rom  0  %  c  t  o  approximately  2  0%c,  with  an  a  verage  o  f 
approximately  13  %c. 
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Wing  Camber ,  %c 

Figure  55:  Effects  of  camber  on  natural  frequency  in  a  homogeneous  Manduca  Sexta  wing. 


3.4.  Vein  Effects 

The  Manduca  Sexta  wing  is  partially  characterized  by  an  intricate  arrangement  of  veins, 
primarily  arranged  i  n  t  he  s  panwise  d  irection.  At  th  is  p  oint,  it  s  eems  th  at  th  e  mo  deshapes  are 
primarily  dom  inated  b  y  planform  i  n  t  he  f  irst  t  wo  m  odes,  a  nd  c  amber  i  n  t  he  t  bird  a  nd  hi  gher 
modes.  S  till,  the  associated  modeshapes  show  significant  midside  motion  not  observed  in  the 
hawkmoth  wings  -  a  likely  indicator  that  the  veins  are  structurally  significant.  F  rom  a  future 
development  standpoint,  a  significant  cost  savings  may  be  realized  by  increasing  the  diameter  of 
the  ve  ins.  R  ecreating  t  he  s  mall-diameter  ve  ins  obs  erved  i  n  t  he  Manduca  Sexta  wing  w  ould 
involve  manufacturing  wing  s  pars  on  t  he  m  icro-scale,  thus  requiring  great  precision  dur  ing 
manufacturing.  With  th  is  v  iew  in  min  d,  in  creasing  th  e  vein  d  iameter  i  s  d  esirable;  doi  ng  s  o 
lowers  production  costs  and  increases  the  wing  volume,  thus  reducing  intra-vein  stress.  These 
considerations  invite  a  study  of  the  structural  dynamic  effects  of  increasing  vein  diameter. 
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The  pr  esence  of  t  hese  ve  ins  causes  one  t  o  ask  several  que  stions  regarding  t  heir 
importance.  Dot  hey  s  erve  i  n  a  pr  imarily  c  irculatory  capacity?  A  sked  a  nother  w  ay,  t  o  w  hat 
extent  do  they  influence  the  structural  response  of  the  wing?  How  significant  is  the  thickness 
(diameter)  of  the  vein  structure?  To  answer  these  questions,  a  flat  model  of  the  Manduca  Sexta 
wing  was  constructed  in  Abaqus.  Although  it  is  true  that  previous  model  results  indicated  that 
camber  a  ffects  th  e  r  esponse  o  f  th  e  s  tructure,  it  is  e  qually  imp  ortant  to  is  olate  t  he  i  mpact  of 
adding  veins  from  that  of  adding  camber. 

3.4.1.  Vein  Contributions  to  Modal  Behavior 

The  flat  wing  model  was  devised  using  the  same  photo-tracing  method  as  the  paper  wing 
model.  T  he  b  earns  ( veins)  ar  e  p  rismatic,  e  ach  h  aving  a  u  nique  cross-section.  Cross  s  ectional 
properties  of  each  vein  were  averaged  using  the  CT  measurements.  T  he  membrane  (denoting 
tissue,  not  an  element  type)  is  also  of  constant  thickness  (0.12  mm)  at  all  points  .  The  three- 
noded  qua dratic  beam  (B32)  and  m embrane  s  ections  (four-noded  general  purpose  s  hell,  type 
S4R)  were  t  ied  t  ogether,  en  forcing  co  ntinuity  o  f  d  isplacement  and  r  otation  a  mong  be  am  a  nd 
shell  elements.  Figure  56  shows  the  assembly  used  during  analysis. 


Figure  56:  Assembly  of  vein  (quadratic  beam)  and  membrane  (Koiter-Sanders  shell)  sections. 


Very  little  variation  in  membrane  thickness  was  observed  during  the  CT  scans. 
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The  inclusion  of  veins  raises  an  important  question  regarding  the  material  properties  of 
the  w  ing  s  tructures.  To  da  te,  no  w  ork  ha  s  be  en  publ  ished  on  t  he  m  aterial  pr  operties  of  t  he 
Manduca  Sexta;  however,  foilwing  their  work  with  the  Tibicen  Canicularis  (cicada),  Mengesha, 
Vallance,  Barraja  and  Mittal  published  material  property  values  for  that  species  [25].  Table  1 1 
gives  the  reported  values  for  the  cicada  wing,  obtained  through  a  nano-indentation  procedure. 
Note  th  at  th  e  th  ickness  of  th  e  c  icada  me  mbrane  is  c  lose  to  th  at  o  f  th  e  Manduca  Sexta  -  a 
heartening  result. 


Table  11.  Material  p  roperties  o  f  a  Tibicen  Canicularis 
(cicada)  r  eported  b  y  Me  ngesha,  V  allance, 
Barraja,  and  Mittal. 


Mass  Density,  g/cm^ 
Young’s  Modulus,  GPa 
Poisson’s  Ratio 
Thickness,  /rm 


Veins 

Membrane 

2.30 

2.30 

1.90 

3.70 

0.495 

0.495 

12.20 

The  Y oung’s  modulus  of  the  membrane  is  greater  than  that  of  the  veins  in  the  cicada 
wing.  This  result  is  not  surprising,  since  many  entomologists  have  noted  that  there  are  two  basic 
classes  of  insect  wings:  s  tiff  an  d  flexible.  T  he  Tibicen  Canicularis  (cicada)  is  of  the  order 
Hemiptera,  and  has  a  rigid  wing  [25],  while  the  Manduca  Sexta  has  a  flexible  wing,  as  do  most 
insects  in  the  order  Lepidoptera  [40].  As  a  r  esult,  the  initial  hawkmoth  model  c  onstructed  in 
Abaqus  uses  the  material  properties  of  the  Tibicen  Canicularis,  but  in  reverse:  t  he  hawkmoth 
veins  are  assumed  to  be  stiffer  than  its  membrane. 

The  quadratic  beam,  Koiter-Sanders  shell  setup  provided  an  accurate  solution  in  Section 
3.2.  B  uilding  on  t  his  approach,  the  wing  shown  in  Figure  56  was  meshed  with  three-noded 
quadratic  b  earns  (  B32),  a  nd  four-noded  general-purpose  shells  (  S4R).  A  Lanczos  m  ethod 
eigenvalue  solver  was  used  to  compute  the  modes  and  modeshapes  of  the  first  three  modes.  The 
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natural  frequencies  obtained  from  this  analysis  are  given  below  in  Table  12.  The  robustness  of 
the  selected  elements  was  evident  from  a  convergence  study  that  showed  virtually  no  change 

23 

with  mesh  refinement  . 


Table  12.  Natural  frequencies  of  a  flat  Manduca  Sexta  wing  model 

Frequency,  Hz 


0)1 

41 

0)2 

91 

OJ3 

225 

0)4 

308 

The  modeshapes  of  the  flat  Manduca  Sexta  wing  model  are  shown  in  the  figures  below. 
By  now,  the  first  two  modeshapes  are  familiar;  the  planform  dictates  first  bending  and  torsion  at 
the  first  and  second  modes,  respectively.  The  third  modeshape  shows  some  deviation  from  the 
fiat  paper  wing  and  its  model.  Instead  of  the  deformations  resembling  pure  second-bending  seen 
previously,  the  beams  influence  the  modeshape  to  add  significant  flattening  near  the  tip .  S  till, 
the  t  bird  m  odeshape  h  as  m  any  aspects  t  hat  ap  pear  s  imilar  t  o  s  econd-bending,  a  nd  i  s  not  s  een 
experimentally  in  the  Manduca  Sexta  wing. 


Figure  57: 


(a)  (b) 

First  modeshape  (41  Hz)  of  a  flat,  FE  wing  with  added  veins,  shown  from  an 
isometric  view  (a)  and  a  tip  view  (b).  Color  bands  denote  displacement. 
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Largest  change  of  any  mode  was  1.5  Hz  when  the  mesh  was  refined  from  348  elements  to  4,121. 
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(a)  (b) 

Figure  58:  Second  modeshape  (92  Hz)  of  a  flat,  FE  wing  with  added  veins,  shown  from 
an  isometric  view  (a)  and  a  tip  view  (b).  Color  bands  denote  displacement. 


Figure  59:  Third  modeshape  (225  Hz)  of  a  flat,  FE  wing  with  added  veins,  shown  from 
an  isometric  view  (a)  and  a  tip  view  (b).  Color  bands  denote  displacement. 

The  fourth  modeshape  is  shown  below,  in  Figure  60,  and  is  similar  to  the  saddle  mode 
seen  previously.  Adding  the  vein  (beam)  structure  has  shifted  the  central  axis  toward  the  trailing 
edge,  and  they  are  now  nearly  parallel.  Accompanying  the  axis  shift  is  a  more  prominent  tip 
deflection  than  previously  observed.  Also  notable  is  the  lack  of  midsection  motion  in  the  leading 
and  trailing  edges,  resembling  the  Manduca  Sexta  more  than  previous  models.  The  midsection 
stiffness  and  the  presence  of  the  third  (extra)  mode  not  seen  in  the  cambered  model  are  indicative 
that  both  venation  and  camber  are  important  structural  characteristics. 
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(a)  (b) 

Figure  60:  Fourth  modeshape  (308  FIz)  of  a  flat,  FE  wing  with  added 
veins,  shown  from  an  isometrie  view  (a)  and  a  tip  view  (b). 


3.4.2.  Effects  of  Vein  Thickness  (Diameter) 

The  preceding  finite  element  case  studies  have  provided  valuable  insight  into  the  level  of 
sophistication  ne  eded  tor  eproduce  a  s  ynthetic  Manduca  Sexta  wing.  From  th  e  r  esults,  it  is 
apparent  t  hat  t  he  w  ing  cannot  c  onsist  of  a  s  ingle,  horn  ogeneous  p  art.  Moreover,  to  a  void  a 
pseudo-homogenous  r  espouse,  t  he  v  eins  m  ust  ha  ve  a  Young’s  m  odulus  g  reater  t  han  the 
membrane,  due  to  their  volume  fractions.  Likewise,  the  finite  element  models  indicate  camber 
is  required  to  prevent  spurious  modes.  From  a  manufacturing  standpoint,  it  is  important  to  gage 
how  s  ensitive  t  he  w  ing  i  s  t  o  c  hanges  i  n  ve  in  thickness  (  diameter).  T  olerance  co  ntrol  an  d 
uniqueness  of  vein  geometry  will  almost  certainly  drive  costs  higher.  Additionally,  should  it  be 
proven  that  the  wing  is  highly  sensitive  to  changes  in  vein  geometry,  techniques  offering  greater 
resolution  than  CT  scans  will  likely  be  required. 

Using  the  flat  wing  model  shown  in  Figure  56,  a  study  of  the  effects  of  vein  thickness  on 
natural  frequency  w  as  c  onducted.  In  t  his  s  tudy,  a  v  ein  t  hickness  f  actor  w  as  de  fined  as  a 
coefficient  that  acted  uniformly  on  all  veins  {e.g.,  a  thickness  factor  of  2  doubled  the  diameter  of 
all  ve  ins).  Throughout  this  s  tudy,  a  r  efmed,  c  onverged  m  esh  w  as  h  eld  unc  hanged,  us  ing 
quadratic  b  earns  an  d  general-purpose  s  hell  elements  to  m  odel  the  ve  ins  a  nd  m  embranes, 
respectively.  Based  up  on  t  he  pr  evious  a  ssertion  t  hat  t  he  Manduca  Sexta  has  a  f  lexible  w  ing 
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structure,  th  e  ma  terial  properties  o  f  Tibicen  Canicularis  veins  a  nd  m  embrane  were  a  gain 
reversed.  The  results  of  the  vein  thickness  study  are  plotted  in  Figure  61.  Examining  the  results 
of  t  he  f  requency  vs  .  v  ein  t  hickness  pi  ot  r  eveals  a  c  ommon  t  heme:  t  he  f  irst  t  wo  na  tural 
frequencies  i  ncrease  n  early  1  inearly  w  ith  v  ein  s  tiffness  ( via  d  iameter),  a  nd  at  t  he  s  ame  r  ate; 
conversely,  t  he  t  bird  m  ode  i  s  s  trongly  a  ffected  b  y  v  ein  t  hickness  ( thus  s  tiffness).  The  t  bird 
modeshape  is  a  second  bending  mode  at  low  stiffness  values,  and  the  fourth  mode  is  a  saddle 
mode.  Increasing  the  vein  diameter  beyond  12  0%  of  the  measured  values  causes  the  second 
bending  m  ode  t  o  di  sappear;  w  hen  t  his  oc  curs,  the  saddle  m  ode  ex  periences  a  r  eduction  i  n 
frequency  and  becomes  the  third  mode. 


Figure  61:  Natural  frequency  as  a  function  of  vein  t  hickness  fora  flat 
wing,  us ing  r  eversed  Tibicen  Canicularis  material  properties. 
A  V ein  t hickness  o f  1  is  equivalent  to  t he  av erage  m ensured 
value  for  each  individual  vein.  Average  vein  radii  vary  from 
0.07  mm  to  0.6  mm. 
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The  V  alidity  o  f  th  e  h  ypothesis  th  at  th  e  ma  terial  p  roperties  o  f  Lepidopterae^"^  are  t  he 
reverse  of  those  of  the  Hemipterae  remains  in  question.  T  o  test  this  position,  the  previous 
study  of  vein  thickness  effects  on  natural  frequency  was  repeated  using  the  material  properties 
from  the  Tibicen  Canicularis  (cicada),  reported  by  Mengesha,  Vallance,  Barraja  and  Mittal  [25]. 
The  results  of  this  analysis  are  plotted  below,  in  Figure  62.  The  stunning  conclusion  is  that  the 
natural  frequencies  are  nearly  identical  when  the  material  properties  are  reversed  (even  though 
the  membrane  is  nearly  twice  as  stiff  as  the  veins  in  the  cicada).  The  same  behavior  as  before  is 
observed:  t  here  are  s  mall  va  nations  i  n  t  he  f  irst  a  nd  s  econd  n  atural  f  requencies  w  ith  v  ein 
thickness,  and  a  strong  correlation  of  the  third  and  fourth  modes  with  vein  thickness.  As  seen 
previously,  increasing  the  vein  thickness  beyond  a  certain  point  results  in  the  elimination  of  the 
second  bending  mode. 


Figure  62:  Natural  frequency  as  a  function  of  vein  t  hickness  fora  flat 
wing,  using  Tibicen  Canicularis  material  properties. 


Order  to  which  the  Manduca  Sexta  belongs 

Order  to  which  the  Tibicen  Canicularis  (cicada)  belongs 
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As  one  m  ight  e  xpect  f  rom  a  t  heoretical  s  tandpoint,  t  he  pr  eeeding  i  nvestigations  of 
material  pr  operties,  c  amber,  a  nd  ve  in  di  mensions  r  eveal  a  pos  sible  non  -uniqueness  of  m  odal 
solutions.  Indeed,  structural  stiffness,  El,  is  a  combination  of  material  property  and  moment  of 
inertia  e  ffects,  i  n  which  a  n  i  nadvertent  unde  restimation  of  one  c  ontributing  f  actor  m  ay  be 
compensated  for  b  y  a  corresponding  ov  erestimation  of  a  nother.  W  ith  Y  oung’s  m  odulus, 
Poisson’s  ratio,  and  density  estimated  (but  unknown),  and  geometric  measurements  limited  by 
the  resolution  of  the  CT  scans,  producing  a  validated  structural  model  of  a  Manduca  Sexta  may 
be  possible,  albeit  non-unique. 
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IV.  Results  and  Discussion 


A  finite  el  ement  model  of  a  Manduca  Sexta  forewing  w  as  c  onstructed  i  n  A  baqus  for 
eomparison  to  experimental  frequency  data.  Throughout  the  model  creation  process,  there  were 
two  competing  requirements:  ma  ximizing  geometric  simplifications  to  ensure  tractability,  and 
minimizing  unnecessary  deviation  from  the  physical  structure  to  preserve  the  accuracy  of  results. 
The  trade  studies  of  Chapter  III  were  heavily  relied  upon  to  produce  a  model  that  most  accurately 
reflected  the  Manduca  Sexta  forewing,  and  to  help  understand  its  shortcomings.  In  light  of  the 
preceding  i nvestigation  of  appropriate  element  s  election,  four-noded,  ge  neral-purpose  ( Koiter- 
Sanders)  shells  were  used  to  model  the  wing  membrane  tissue,  and  three-noded  quadratic  beams 
served  to  represent  the  vein  structure. 

This  chapter  is  organized  chronologically  beginning  with  model  construction,  to  include 
material  properties.  Frequency  analysis  and  convergence  study  results  are  shown,  followed  by  a 
discussion  of  the  model  limitations.  Recommendations  are  made  regarding  possible  methods  to 
improve  the  modeling  process.  With  the  model  limitations  in  mind,  a  scaling  factor  is  developed 
that  normalizes  the  Abaqus  and  Polytech  scanning  laser  outputs,  and  a  stress  analysis  of  the  first 
two  modes  is  conducted. 

4.1.  Model  Construction 

The  finite  e  lement  studies  i  n  Chapter  I II  reveal  t  hat  t  he  n  atural  f  requencies  o  f  t  he 
Manduca  Sexta  wing  are  strongly  correlated  to  camber  and  vein  thickness  (geometry),  but  only 
weakly  co  rrelated  tom  aterial  p  roperties  w  ithin  t  he  r  ange  of  2  -5  GPa.  C  onsequently,  both 
camber  and  vein  thickness  variation  as  a  function  of  span  were  included  in  the  finite  element 
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model.  It  was  also  assumed  (without  verification  )  that  the  material  properties  of  the  Manduca 
Sexta  are  reasonably  close  to  those  of  the  Tibicen  Canicularis. 

An  average  cam  her  value  was  c  omputed  ba  sed  on  C  T  m  easurements,  a  nd  f  ormed  t  he 
foundation  of  t  he  m  odel.  U  sing  t  he  a  pproach  de  scribed  i  n  Chapter  III ,  a  pr  ismatic  s  hell  o  f 
constant  c urvature  was  extruded,  c rearing  a  rectangular  pi  anform  w ith  constant  thickness.  A 
planform  sketch  of  the  Manduca  Sexta  wing  was  then  projected  onto  the  surface  of  the  extruded 
shell,  as  shown  in  Figure  63.  Next,  an  extruded  cut  was  performed  to  eliminate  material  outside 
the  wing  planform,  resulting  in  a  similar  wing  to  the  one  shown  previously  (see  Figure  56). 


Figure  63:  Planform  sketch  of  a  Manduca  Sexta  wing,  projected 
on  to  a  curved  surface  to  simulate  an  average  camber. 


Preliminary  nano-indentation  studies  by  O’Hara  place  the  Young’s  modulus  of  dry,  aged  Manduca  Sexta  wings 
between  2  and  5  GPa  [31].  Material  properties  for  the  veins  have  yet  to  be  determined. 
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The  ve  nation  s  hown  i  n  Figure  63  eonsists  onl  y  of  p  artition  s  ketehes,  not  pa  rts,  that 
effeets  only  the  mesh  eonstruction.  Unfortunately,  Abaqus/CAE  currently  lacks  the  capability  to 
convert  partition  geometry  to  wire  parts  .  To  circumvent  this  problem,  a  fine  mesh  of  general- 
purpose  shell  elements  (S4R)  was  created  on  the  part  shown  in  Figure  63.  Using  the  script  given 
in  Appendix  A  -  CoordGen,  the  coordinates  of  the  nodes  corresponding  to  the  locations  of  the 
vein  partitions  were  extracted  from  the  Abaqus  input  files.  These  coordinates  were  imported  into 
SolidWorks  to  create  splines  that  were  then  exported  as  wire  parts.  With  each  of  the  veins  saved 
as  wire  parts  in  a  standardized  format  (IGES),  they  were  imported  into  Abaqus  and  mated  with 
the  partition  sketches  from  which  they  were  created.  Consistent  with  the  vein-membrane  models 
of  Chapter  III,  the  beams  and  their  adjacent  shells  are  held  together  using  the  “tie”  constraint, 
forcing  continuity  of  displacement  at  the  beam  nodes. 

Beam  sections  (profiles)  were  created  using  the  CT  cross-section  measurements.  E  ach 
vein  (beam)  cross-section  changes  at  different  spanwise  locations,  and  tapers  at  its  own,  unique 
rate.  To  reflect  t  his  be  havior  i  n  t  he  finite  e  lement  model,  ve  ins  w  ere  s  ectioned  a  t  i  rregular 
intervals  based  on  changes  in  the  CT  measurements.  B  etween  each  of  the  section  locations,  an 
average  value  of  vein  thickness  was  computed,  and  a  new  section  was  created  when  the  vein 
thickness  varied  by  more  than  10%  from  the  average  section  value.  Since  Abaqus/CAE  lacks  the 
capability  t  o  cr  cate  t  apered  s  ections,  b  earn  cr  oss-section  c  hanges  oc  curred  di  scontinuously. 
Eigure  64,  below,  shows  the  Manduca  Sexta  wing  model  color-coded  by  cross-section. 

The  b  earns  (veins)  were  tie  d  a  t  th  e  mid  section  of  th  e  s  hells  (membrane  tis  sue)  in  th  e 
model  of  Eigure  64  for  two  reasons.  Recall  Eigure  19,  the  CT  image  of  a  wing  cross-section  near 
the  root.  Erom  this  figure  it  is  evident  that  some  of  the  veins  pass  through  the  midsection  of  the 
membrane  s  tructure,  w  hile  ot  hers  do  not  .  T  hose  ve  ins  t  hat  do  a  ppear  t  o  ha  ng  be  low  t  he 
The  author  has  requested  this  feature,  and  Simulia  has  agreed  to  include  it  in  a  future  release  of  Abaqus/CAE  [30]. 
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membrane  do  s  o  near  the  root,  and  are  gradually  engulfed  by  the  membrane  as  the  vein  tapers 
with  span.  Second,  all  veins,  at  all  points,  are  surrounded  by  a  thin  coating  of  membrane  in  the 
actual  Manduca  Sexta  wing.  When  modeling  the  membrane  path  (over  or  through)  the  veins,  it 
was  assumed  that  the  thin  layers  of  membrane  above  and  below  the  veins  were  of  roughly  the 
same  mass  moment  of  inertia  as  an  idealized  membrane  passing  through  the  midsection. 


Figure  64:  Planform  view  of  a  finite  element  model  of  a  Manduca 
Sexta  forewing,  color  coded  by  cross-section  dimension. 

In  light  of  the  previous  studies  asserting  that  the  Manduca  Sexta  wing  is  far  more  flexible 
than  the  Tibicen  Canicularis  wing,  and  O  ’Hara’s  preliminary  findings  regarding  the  membrane 
modulus  [31],  the  material  properties  used  mirrored  the  reversed  cicada  properties.  With  the 
exception  of  the  modulus  of  elasticity  of  the  veins,  all  properties  are  the  same  as  the  Tibicen 
Canicularis',  here,  Y oung’s  modulus  was  increased  s lightly  (7.5%),  from  3.7  G  Pa  to  4.0  GPa. 
Table  13  shows  the  material  properties  used  during  model  construction. 
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Table  13.  Assumed  values  of  material  properties  fora 
Manduca  Sexta  (hawkmoth)  wing. 

Veins  Membrane 


Density,  g/cwd 

2.30 

2.30 

Young’s  Modulus,  GPa 

4.00 

1.90 

Poisson’s  Ratio 

0.495 

0.495 

Thickness,  /rm 

12 

4.2.  Mass  Properties 

During  the  experimental  phase  of  his  researeh  Norris  demonstrated  that  the  camber  and 
frequency  response  of  the  Manduca  Sexta  wing  change  dramatically  after  approximately  three 
hours,  with  mass  being  less  dramatically  effected  [29].  Asa  result,  preserving  the  freshness  of 
the  wing  during  CT  and  frequency  investigations  was  of  the  upmost  priority.  Estimation  of  the 
wing  mass  was  conducted  using  fresh-wing  data  previously  acquired  by  Norris. 

In  the  initial  phase  of  his  work  at  AFIT,  Norris  only  conducted  vibration  experiments  in 
air.  D  uring  t  his  t  ime,  one  of  hi  s  primary  obj  ectives  w  as  t  o  de  velop  a  statistical  d  atabase  o  f 
Manduca  Sexta  wing  properties.  Using  his  data,  the  average  wing  mass  was  0.0416  grams,  with 
a  standard  deviation  of  0.0037  grams.  In  comparison,  Abaqus  determined  the  mass  of  the  finite 
element  model  to  be  0.040  grams,  well  within  one  standard  deviation  of  the  average  Manduca 
Sexta  wing  m  ass.  This  r  esult  1  ends  f  urther  c  redence  t  o  t  he  a  ssumption  t  hat  t  he  m  aterial 
properties  of  the  Tibicen  Canicularis  (cicada)  are  likely  close  to  those  of  the  Manduca  Sexta. 

4. 3.  Frequency  Analysis  Results 

The  wing  model  shown  in  Figure  64  was  meshed  with  three-noded  quadratic  beams  and 
four-noded,  doubl  y-curved,  s  hells,  us  ing  a  r  educed  i  ntegration  f  ormulation  of  1 1  t  hrough  t  he 
thickness  integration  points.  The  first  three  natural  frequencies  of  this  model  (converged  values) 
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are  presented  in  Table  14,  and  compared  to  the  experimentally  obtained  values  of  the  Manduca 
Sexta  forewing. 


Table  14.  Natural  frequency  results  of  a  Manduca  Sexta  FE 
model,  and  comparison  to  experimental  values. 


Experimental,  Hz 

FE  Model,  Hz 

Minimum  Difference 

OJl 

86  ±2 

84.6 

0.0% 

(x)2 

106  ±2 

106.1 

0.0% 

"3 

155  ±2 

317.7 

102.4% 

Table  15.  Comparison  of  frequency  ratios  of  Manduca  Sexta  and  FE  wings. 

Experimental,  Hz 

FE  Model,  Hz 

Minimum  Difference 

1.23  ±0.05 

1.25 

0.0% 

tog/nji  1.80  ±0.07 

3.91 

100.8% 

It  is  immediately  obvious  that  the  first  two  natural  frequencies  of  the  Manduca  Sexta 
wing  and  the  finite  element  model  match  very  well,  as  does  their  ratio.  The  finite  element  model 
does  not  adequately  match  the  third  mode,  however.  The  predicted  third  natural  frequency  and 
its  ratio  to  th  e  f  ist  mo  de  d  iffers  from  e  xperimental  v  alues  by  a  factor  of  roughly  two.  As 
previous  s  tudies  ha  ve  s  hown,  t  he  t  bird  m  ode  i  s  t  he  m  ost  s  ensitive  (of  t  he  t  hree  m  odes 
considered)  to  s  tiffness  e  ffects  r  esulting  from  camber,  c  ross-sectional  s  hape,  and  m  aterial 
property  variation. 

Not  surprisingly,  the  first  two  modeshapes  of  the  finite  element  model  correlate  well  with 
the  experimental  results  seen  in  the  Manduca  Sexta  wing.  The  first  modeshape  of  the  finite 
element  wing  is  plotted  below,  in  Figure  65  .  This  modeshape  is  the  familiar  first  bending 


For  readability’s  sake,  the  most  refined  (converged)  mesh  is  not  s  hown  in  the  figures.  T  heconverged  mesh 
contains  s  o  many  e  lements  t  hat  t  he  c  olor  ba  nds  de  noting  di  splacement  a  re  no  1  onger  v  isible  d  ue  t  o  t  he  bl  ack 
outline  of  the  elements. 
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mode,  w  ith  c  olor  ba  nds  de  noting  m  agnitude  of  displacement.  Ass  een  e  xperimentally  i  n  t  he 
Manduca  Sexta  wing,  the  wing  displacement  is  a  function  of  span,  with  little  variation  in  the 
chordwise  direction. 


(a)  (b) 


Figure  65:  First  modeshape  (84.6  Flz)  of  FE  Manduca  Sexta  wing 

model,  shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 

The  s  econd  m  odeshape  is  a  t  orsional  m  ode,  a  nd  i  s  pi  otted  i  n  Figure  66.  Unlike  t  he 
torsional  modes  seen  in  Chapter  III,  the  second  modeshape  of  this  model  does  not  affect  the  wing 
in  a  rigid,  plate-like  fashion.  That  is,  instead  of  the  entire  wing  (or  a  majority  of  it)  rotating  as  a 
single,  planar  unit  about  a  central  axis,  the  wing  of  Figure  66  shows  significant  curvature  in  the 
tip  view.  This  deflection  is  consistent  with  the  second  modeshape  of  the  Manduca  Sexta  wing, 
shown  in  Figure  23.  In  both  cases,  the  central  axis,  about  which  the  rotation  occurs,  is  offset 
slightly  forward  of  the  half-chord  line.  M  oreover,  the  finite  element  model  accurately  predicts 
that  the  vertical  displacement  of  the  trailing  edge  has  a  slightly  greater  magnitude  than  that  of  the 
leading  edge. 


(a)  (b) 


Figure  66:  Second  modeshape  (106. 1  Hz)  of  FE  Manduca  Sexta  wing 
model,  shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 
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The  third  modeshape  is  a  saddle-mode  and  is  shown  below,  in  Figure  67.  Although  the 
frequency  value  predicted  by  the  finite  element  model  does  not  coincide  with  the  155  Hz  value 
observed  i  n  t  he  ha  wkmoth  w  ing,  i  ts  m  odeshape  i  s  s  tunningly  s  imilar.  Like  t  he  s  addle  m  ode 
shown  in  Figure  24,  the  finite  element  model  predicts  a  smooth  deflection  of  nearly  constant 
radius  f  rom  t  he  t  ip  vi  ew.  Likewise,  t  he  chordwise  de  flection  i  s  a  ngled,  sot  hat  t  he  a  xis  of 
rotation  is  slightly  off  center,  giving  the  trailing  edge  the  greatest  magnitude  of  deflection  in  the 
vertical  direction.  This  effect  causes  the  trailing  edge  slope  to  be  greater  than  that  of  the  leading 
edge;  as  a  result,  the  root  is  visible  in  the  upper  left  comer  of  Figure  67  (b)  -  mimicking  the 
Manduca  Sexta ’s  response  in  Figure  24  (b). 


(a)  (b) 


Figure  67:  Third  modeshape  (3 17.7  Hz)  of  FE  Manduca  Sexta  wing 

model,  shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 

4.4.  Convergence  Studies 

Convergence  studies  were  carried  out  to  ensure  accuracy  of  the  finite  element  solutions. 
Mesh  sizes  between  295  and  80,633  elements  were  investigated,  and  the  results  of  these  studies 
are  plotted  in  the  figures  below.  For  each  mode,  the  trend  observed  in  the  convergence  plots  is 
the  same:  as  the  mesh  size  is  increased,  the  solution  obtained  by  the  eigenvalue  solver  decreases 
to  some  asymptotic  value.  In  this  analysis,  the  mesh  was  considered  sufficiently  refined  when  a 
large  change  in  the  number  of  elements  produced  a  frequency  shift  of  less  than  one  percent  in 
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any  mode.  During  the  eonvergenee  studies,  the  number  of  beams  and  shells  was  eontrolled  by 
assigning  a  g  lobal  s  eed  s  ize,  t  bus  e  ontrolling  t  he  1  ength  of  a  ny  s  ingle  e  lement  a  long  its  e  dge. 
Curvature  eontrol  was  enforeed  to  ensure  that  highly  eurved  surfaees  were  adequately  meshed. 

A  solution  satisfying  the  above  eonvergenee  eriterion  was  found  with  a  mesh  eonsisting 
of  three-noded  quadratie  beam  elements  (B32)  and  four-noded  general-purpose  shell  elements. 
The  eonverged  mesh  used  a  total  of  80,633  elements  (83,664  nodes).  The  seeond  most  refined 
mesh  contained  13,969  elements  (15,330  nodes),  and  had  a  maximum  frequency  shift  of  1.0%  (in 
the  third  mode)  from  the  converged  solution.  Since  the  most  refined  mesh  contained  more  than 
six  t  imes  t  he  n  umber  o  f  el  ements  as  t  he  s  econd  m  ost  r  elined  m  esh  but  1  acked  s  ignilicant 
variation  i  n  f  requency  values,  t  he  s  olution  w  as  c  onsidered  s  ufficiently  converged.  From  a 
computational  standpoint,  an  efficient  solution  was  obtained  with  3,850  elements  (4,570  nodes), 
at  w  hich  poi  nt  a  11  m  odes  w  ere  w  ithin  1 .5%  of  t  he  c  onverged  va  lue  obt  ained  w  ith  80,633 
elements.  Likewise,  the  robustness  of  the  selected  elements  is  evident  in  the  solution  obtained 
with  only  636  elements  (896  nodes);  here,  the  largest  difference  with  the  converged  solution  is 
2.8%,  observed  in  the  second  mode.  Table  16  indicates  the  number  of  elements  and  nodes  used 
during  the  convergence  study,  and  how  many  of  these  elements  were  beams  (B32)  and  shells 
(S4R). 
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Table  16.  Number  of  elements  and  nodes  (and  of  which 
type)  used  during  convergence  studies. 


Total  Elements 

Nodes 

Beams  (B32) 

Shells  (S4R) 

251 

389 

116 

135 

316 

481 

138 

178 

428 

627 

167 

261 

636 

896 

221 

415 

1,199 

1,582 

332 

867 

3,850 

4,570 

658 

3,192 

5,754 

4,774 

823 

4,931 

13,969 

15,330 

1,314 

12,655 

80,633 

83,664 

3,285 

77,348 
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Figure  68:  Convergence  study  results  for  the  first  mode  of  a  Manduca 
Sexta  finite  element  model. 
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Figure  69:  Convergence  st  udy  r  esults  f  or  the  seco  nd  m  ode  of  a 
Manduca  Sexta  finite  element  model. 


Figure  70:  Convergence  s  tudy  r  esults  f  or  the  t  bird  m  ode  of  a 
Manduca  Sexta  finite  element  model. 
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4. 5.  Model  Lint  itations 


Before  undertaking  a  n  umerieal  analysis  of  the  Manduca  Sexta  wing,  one  of  the  stated 
objeetives  w as  t o  de  velop  a  finite  e  lement  model  eap  able  o  f  ad equately  r  epresenting  the  first 
three  natural  frequeneies  and  their  as soeiated  modeshapes.  The  above  results  indicate  that  the 
finite  element  model  produced  acceptable  results  for  the  first  two  frequencies  and  modeshapes. 
While  t  he  t  bird  m  odeshape  a  ppears  qu  alitatively  similar  t  o  t  he  s  addle  mode  obs  erved  i  n  t  he 
Manduca  Sexta  wing  experiment,  the  frequency  values  are  different  by  approximately  a  factor  of 
two.  T  here  are  several  possible  reasons  why  this  deviation  might  exist:  low  resolution  during 
component  di  mensioning,  i  nadequate  c  haracterization  o  f  w  ing  c  amber,  and  poor  e  stimation  of 
material  properties.  Each  of  these  possibilities  is  explored  below. 

4. 5. 1.  Component  Dimensioning 

The  di  mensions  of  t  he  ve  in  a  nd  m  embrane  s  tructures  w  ere  obt  ained  vi  a  c  omputed 
tomography  (CT)  measurements.  During  this  process  it  was  observed  that  the  vein  dimensions 
(and  to  a  lesser  extent  the  membrane  thickness)  area  function  of  span.  As  the  cross-section 
location  moves  from  root  to  tip  the  veins  taper,  and  the  low  resolution  of  the  CT  imager  prohibits 
accurate  measurement.  Figure  71  compares  images  of  the  first  three  veins  (leading  edge  is  to  the 
left)  at  the  wing  root,  and  approximately  the  half-span  location  (20  mm  from  the  root).  Due  to 
the  change  in  vein  dimensions,  the  differences  in  image  quality  are  striking.  Gone  from  the  half¬ 
span  i  mage  i  s  t  he  cl  ear  i  nner-diameter,  an  d  t  he  separation  of  t  he  first  tw  o  v  eins  is 
indistinguishable. 
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Figure  71:  Comparison  of  first  three  veins  at  the  wing  root  (a)  and  20  mm  (approximately 
half-span)  from  the  root  (b).  In  both  cases,  the  leading  edge  is  to  the  image  left. 

Figure  70  (  b)  hi  ghlights  t  he  e  ffect  of  de  ereasing  t  he  ve  in  di  ameter  a  t  eonstant  i  mage 
resolution:  the  relative  importanee  of  a  single  pixel  is  dramatically  elevated.  The  pixel  size  is 
the  limit  ing  v  alue  o  f  p  recision  in  d  igital  ima  ge  me  asurement,  th  us  a  ttempting  to  d  imension 
structures  with  length  scales  of  the  same  order  of  magnitude  as  the  pixel  size  greatly  increases 
the  probability  of  measurement  errors.  The  trend  observed  in  Figure  70  continues  as  the  cross- 
section  1  ocation  a  dvances  t  oward  t  he  t  ip.  A  s  t  he  ve  ins  t  aper,  t  hey  b  ecome  1  ess  vi  sible,  a  nd 
finally  appear  to  blend  into  the  membrane  prior  to  their  terminal  locations. 

The  effects  of  varying  vein  dimensions  on  frequency  content  were  studied  in  Chapter  III. 
These  p  arametric  s  tudies  d  emonstrated  t  hat  a  11  m  odes  are  s  trongly  affected  b  y  even  s  mall 
changes  in  vein  diameter.  As  discussed  in  Chapter  III,  the  first  two  modes  are  linear  functions  of 
vein  thickness,  while  the  third  mode  can  be  approximated  by  a  polynomial  that  changes  rapidly 
with  d  iameter  v  ariation.  A  n  otable  r  esult  f  rom  th  is  p  arametric  s  tudy  is  th  at  th  e  f  requency 
separation  between  the  first  two  modes  appears  fixed,  regardless  of  vein  diameter. 

Comparing  the  results  of  the  vein  thickness  (diameter)  studies  to  the  results  of  the  finite 
element  model  of  the  Manduca  Sexta  wing  suggests  the  outer  diameter  of  the  veins  may  have 
been  overestimated.  The  higher-quality  CT  images  near  the  wing  root  reveal  that  the  veins  are 
surrounded  by  a  coating  of  membrane  tissue.  As  the  cross-section  location  moves  from  root  to 
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tip,  the  relative  thiekness  of  the  membrane  eoating  inereases  due  to  a  eorresponding  deerease  in 
vein  thickness. 

It  is  likely  that  low  CT  image  quality  resulted  in  measuring  the  total  diameter  of  what 
appeared  t  o  be  t  he  ve  in  s  tructure,  instead  of  di  mensioning  onl  y  t  he  ve  in  and  e  xcluding  t  he 
surrounding  membrane  (see  Figure  72,  below).  Since  the  assumed  value  of  Young’s  modulus  of 
the  veins  (beams)  was  more  than  twice  that  of  the  membrane  tissue  (shells),  measuring  the  total 
diameter  s  urrounding  t  he  ve  in  would  result  in  a  significant  s  tiffening  ef  fect^^.  Conversely, 
because  the  veins  and  membrane  are  assumed  to  be  of  the  same  density,  over-measuring  the  vein 
diameter  would  not  increase  the  mass  of  the  model  -  making  measurement  errors  more  difficult 
to  detect.  Thus,  if  a  more  accurate  model  is  to  be  produced,  a  geometry  measurement  technique 
that  offers  higher  resolution  than  the  CT  imager  is  necessary. 


Figure  72:  Schematic  of  membrane  tissue  surrounding  the  vein.  Low  CT  resolution  results 
in  measuring  the  outer-diameter  of  both  stmctures,  instead  of  the  vein  only. 


This  extent  to  which  this  over-dimensioning  of  vein  diameter  affects  frequency  results  increases  with  the  number 
of  poorly-dimensioned  veins.  If  the  CT  resolution  were  such  that  only  a  few  vein  segments  had  questionable 
dimensions,  this  effect  might  be  negligible;  this  is  not  the  case  here. 
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4.5.2.  Modeling  Camber 


Camber  is  an  important  feature  of  many  wings,  and  the  Manduca  Sexta  forewing  is  no 
exeeption.  While  o  ften  e  xamined  for  reasons  relating  t  o  fluid  dynamies,  t  he  finite  e  lement 
results  of  Chapter  III  demonstrate  that  eamber  also  has  an  important  role  in  struetural  dynamies. 
Adding  eurvature  to  the  Manduca  Sexta  wing  adds  stiffness,  whieh  in  turn  inereases  its  natural 
frequeneies.  The  approaeh  taken  during  the  finite  element  modeling  proeess  was  to  model  the 
wing  with  a  shell  of  eonstant  radius  and  thickness.  The  selected  camber  was  an  average  value 
based  on  C  T  measurements,  and  was  measured  using  the  NACA  wing  standard;  the  maximum 
distance  between  the  chord  line  and  the  mean  camber  line  is  the  camber  value. 

The  curvature  of  a  Manduca  Sexta  wing  is  far  more  complex  than  that  of  wings  typical  of 
man-made  d  esigns.  Instead  of  t  he  s  mooth  w  ing  s  urface  o  f  c  onstant  c  oncavity  seen  i  n  m  ost 
commercial  w  ing  d  esigns,  t  he  Manduca  Sexta  wing  i  s  of  pos  hive  c  oncavity  n  ear  t  he  1  eading 
edge,  and  ne gative  c  oncavity  ov  er  the  majority  of  t he  w ing  s urface,  a s  shown  in  Figure  19. 
Additionally,  the  camber  is  a  function  of  span  in  the  Manduca  Sexta  wing.  Future  modeling 
efforts  should  include  this  geometric  sophistication  if  high-frequency  modes  are  to  be  matched. 

4. 5. 3.  Material  Properties 

No  material  properties  for  the  Manduca  Sexta  (or  monarch  butterfly,  a  similar  species  of 
the  order  Lepidoptera)  are  currently  published,  leaving  one  forced  to  estimate  Young’s  modulus, 
Poisson’s  ratio,  and  density.  As  a  result,  the  material  properties  used  in  the  finite  element  model 
were  based  on  those  of  the  Tibicen  Canicularis  (cicada),  as  reported  b  y  Mengesha,  V  allance, 
Barraja  and  Mittal.  The  cicada  is  of  the  order  Hemiptera,  a  cl  ass  of  species  observed  to  have 
stiff,  semi-rigid  wings;  conversely,  the  Manduca  Sexta  (hawkmoth)  is  of  the  order  Lepidoptera, 
and  is  known  to  have  more  flexible  wings  than  Hemipterae  [25].  Based  on  this  knowledge,  the 
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assumption  t  hat  t  he  Y  oung’s  m  odulus  of  t  he  Manduca  Sexta  wing  c  omponents  is  lik  ely  th  e 
reverse  of  the  Tibicen  Canicularis  was  made. 

Stiffness  is  a  funetion  of  Young’s  modulus,  geometry,  and  mass-density.  While  the  mass 
of  the  finite  element  model  was  representative  of  what  one  might  expeet  from  a  typical  Manduca 
Sexta  wing,  it  is  possible  that  the  agreement  with  the  first  two  experimental  modes  was  the  result 
of  the  non-uniqueness  of  the  stiffness  value.  In  other  words,  the  vein  and  membrane  densities  of 
the  hawkmoth  might  be  different  than  those  of  the  cicada.  If  this  is  the  case,  assuming  the  same 
values  as  the  cicada  forced  the  moduli  of  elasticity  to  also  reflect  the  values  reported  in  Ref  [25]. 
Thus,  it  may  be  said  that  the  ratio  of  Y oung’s  modulus  to  mass  density,  E/ju,  is  an  important 
characteristic  i  n  t  he  f  requency  r  esponse  o  f  t  he  Manduca  Sexta  forewing.  W  hile  t  he  m  odd 
produces  good  results  for  the  first  two  modes,  further  investigation  into  the  material  properties  is 
required  t  o  de  termine  w  hether  t  he  va  lues  of  E  and  ji  are  co  rrect,  o  r  w  hether  t  heir  r  atio  w  as 
arbitrarily  satisfied. 

The  f  act  t  hat  t  he  n  atural  f  requencies  o  f  t  he  finite  element  model  a  gree  with  t  he 
experimental  results  only  in  the  first  two  modes  likely  indicates  some  total  stiffness  error  in  the 
model.  Indeed,  a  fundamental  concept  in  vibration  theory  is  the  effect  of  density  and  stiffness  on 
frequency.  A  more  accurate  finite  element  model  of  the  Manduca  Sexta  will  likely  require  one 
to  obtain  material  properties  for  the  hawkmoth. 

4. 6.  Modeshape  Stresses 

It  is  imp  ortant  to  q  uantify  th  e  s  tresses  th  e  Manduca  Sexta  wing  un  dergoes  dur  ing 
operation.  A  Ithough  the  magnitude  of  velocity  and  displacement  at  the  modes  may  be  greater 
than  typically  experienced  in  normal  operation,  the  wing  is  evidently  capable  of  sustaining  such 
large  oscillations.  From  a  design  life  and  fatigue  analysis  standpoint,  a  stress  analysis  of  the 
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observed  operational  defleetion  shapes  will  likely  play  an  important  role  in  future  MAV  design. 
As  discussed  in  the  preceding  sections,  the  final  finite  element  model  produced  accurate  results 
for  the  first  two  modes,  but  was  unable  to  accurately  predict  the  third  natural  frequency  value. 
With  this  in  mind,  a  von  Mises  stress  analysis  of  the  first  two  operational  deflections  shapes  was 
conducted;  the  third  modeshape  was  omitted  from  this  analysis  due  to  the  lack  of  confidence  in 
the  results. 

4. 6. 1.  Scaling  of  Results 

An  important  result  from  the  experimental  work  detailed  in  Chapter  II  was  the  structural 
linearity  observed  in  the  Manduca  Sexta  wing.  By  definition  of  a  linear  system,  the  Manduca 
Sexta  wing  satisfies  additivity  and  homogeneity  of  response  to  an  input.  The  modeshapes  in 
Abaqus  are  normalized  to  report  some  relative  value  of  displacement  at  the  nodes  for  a  mode  of 
vibration,  however  this  value  is  different  than  the  displacement  reported  by  the  Polytec  scanning 
laser  vi  brometer.  A  Ithough  a  m  odeshape  i  s  magnitude  in  dependent  in  its  ma  thematical 
formulation,  the  Manduca  Sexta  forewing  was  underwent  oscillations  of  a  known  magnitude.  By 
comparing  t  he  wing  d  efiections  r  ecorded  dur  ing  t  he  e  xperimental  p  hase  t  o  t  he  r  dative 
displacements  predicted  by  Abaqus,  one  can  determine  what  stresses  were  present  in  the  wing  at 
the  time  of  its  vibration. 

In  0  rder  t  o  co  rrelate  t  he  an  alytical  an  d  experimental  d  isplacements  (  and  s  tresses),  a 
proportional  scaling  factor,  1,  must  be  developed.  By  the  property  of  homogeneity,  this  scaling 
factor  will  serve  as  a  coefficient  to  correlate  stresses  reported  in  Abaqus  to  those  experienced  in 
the  w  ing.  This  s  caling  method  w  ill  o  nly  b  e  u  sed  f or  a  nalyzing  s  tress  in  th  e  first  tw  o  mo  des 
where  the  frequency  results  of  the  finite  element  model  and  the  experimental  wing  are  close. 
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Since  t  he  t  hird  m  ode  i  s  not  a  dequately  r  epresented  i  n  the  finite  e  lement  model,  it  w  ill  b  e 


exeluded  from  stress  analysis  using  the  sealing  faetor,  A. 

Invoking  the  property  of  superposition,  the  scaling  factor.  A,  is  simply  the  ratio  of  the 
wing  tip  displaeements  from  the  Polytee  output  and  Abaqus  model: 


^PSV 

^Abq 


(7) 


where  5psv  is  the  wing  tip  displaeement  measured  by  the  Polytee  seanning  laser  vibrometer,  and 
dAbq  is  the  wing  tip  displaeement  reported  by  Abaqus.  Using  Eq.  (7),  sealing  faetors  for  the  first 
bending  mode,  Ai,  and  for  the  torsional  mode,  A2,  are  eomputed. 

The  tip  d  isplaeement  in  th  e  f  irst  mo  de  r  eported  b  y  th  e  Polyteeh  s  eanning  vi  brometer 
(PSV)  software  was  27.29  xlO'^  meters,  while  Abaqus  reported  a  first  mode  displaeement  of  3.09 
xlO'  meters.  Inserting  these  values  into  Eq.  (7), 

27.29x10“®  m 

Ai  = - 

3.09x10  ^m 
Ai  =  8.92x10-3 

Similarly,  Abaqus  reported  the  tip  displacement  in  the  seeond  mode  as  1.98  xl O'  meters, 
while  t  he  1  aservi  brometer  r  eported  a  s  eeond  m  ode  di  splaeemento  f3.32  xlO'^  meters. 
Substituting  these  values  into  Eq.  (7)  gives. 


3.32x10-®  m 
A,  = - 

1.98x10  3  m  ^9^ 

A2  =  1.68x10-3 

With  the  s  ealing  f  aetors  (Ai  and  A2)  obtained,  the  m  odel  s  tress  v  alues  ean  her  elated  1 0  t  he 
stresses  in  the  Manduca  Sexta  wing. 


92 


4. 6. 2.  Wing  Stress  Analysis  -  First  Bending  Mode 


A  von  Mises  stress  analysis  of  the  first  operational  deflection  shape  (first  bending,  84.6 
Hz)w  as  c  onductedus  ingt  he  finite  e  lement  model  described  pr  eviously^^.  Since  t  he 
proportional  limit  of  the  materials  used  in  the  wing  are  unknown,  this  analysis  assumes  purely 
elastic  behavior. 

The  maximum  (scaled  by  X\)  von  Mises  stress  is  258.4  KPa,  and  occurs  at  the  trailing 
edge  near  the  clamp.  The  scaled  minimum  stress  is  567  P  a  (0.08  psi  -  essentially  zero  stress), 
and  occurs  at  the  wing  tip.  Stress  contour  plots  given  in  Abaqus  indicate  that  the  veins  appear  to 
be  the  primary  load-bearing  members,  and  have  an  average  stress  of  approximately  three  orders 
of  magnitude  greater  than  that  of  the  membrane.  The  most  striking  feature  is  the  nearly  uniform 
low  stress  in  the  membrane;  only  portions  near  the  clamped  boundary  -  and  a  small  section  at 
approximately  one  -third  o  f  th  e  w  ing  s  pan  -  are  s  ignificant  s  tress  c  oncentrations  v  isible. 
Likewise,  it  is  apparent  that  the  stress  in  the  veins  is  the  highest  in  approximately  the  first  third 
of  the  s  pan.  It  w  as  o  bserved  t  hat  v  ein  i  ntersections  i  n  w  hich  ne  w  ve  ins  a  re  c  reated  t  end  t  o 
increase  the  von  Mises  stress  to  approximately  65  KPa  (scaled  value),  while  intersections  that 
combine  existing  veins  tend  to  have  little  effect  on  stress  concentrations. 

4. 6. 3.  Wing  Stress  Analysis  -  Torsional  Mode 

The  stress  analysis  was  repeated  for  the  second  (torsional)  m ode,  this  time  scaling  the 
results  b  y  th  e  to  rsional  f  actor,  2,2.  The  m  aximum  von  M  ises  s  tress  obs  erved  i  n  t  he  t  orsional 
eigenvector  is  28.91  KPa,  and  occurs  in  the  leading  edge  vein  at  approximately  80%  of  the  span. 
The  minimum  stress  is  located  in  the  membrane  near  the  leading  edge  and  was  found  to  be  277 
Pa  -  approximately  half  that  of  the  first  mode.  In  general,  the  second  modeshape  causes  a  lower 

Contour  plots  of  this  analysis  are  given  in  Appendix  E  -  Stress  Plots. 
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average  s  tate  o  f  s  tress  than  t  he  f  irst  m  ode,  with  less  variation  be  tween  t  he  minimum  and 
maximum  stress  values. 

A  more  even  stress  distribution  is  present  in  the  torsional  mode  than  in  the  first  bending 
mode,  but  many  of  the  trends  are  similar.  The  maximum  stresses  are  loeated  near  the  clamp  and 
at  vein  intersections,  and  generally  decrease  with  span.  In  a  noticeable  departure  from  this  trend, 
the  t  hickening  of  t  he  first  ve  in  s  lightly  b  eyond  the  ha  If-span  r  esults  in  st  ress  reduction;  t  his 
reduction  appears  t  o  also  c  ause  a  1  ocal  s  tress  r  eduction  i  n  t  he  ve  ins  be  hind  it.  A  s  t  he  c  ross- 
section  of  the  leading  edge  vein  narrows  near  the  tip,  the  stress  increases  to  24.88  KPa  (scaled 
value),  reflecting  the  large  leading-edge  displacement  in  the  torsional  mode.  Unlike  the  behavior 
observed  i  n  t  he  f  irst  m  ode,  t  he  ve  in  or  iented  i  n  t  he  c  hordwise  di  rection  pi  ays  a  dr  amatic 
stiffening  role  in  the  torsional  mode.  Indeed,  its  stresses  are  much  higher  than  observed  in  the 
bending  mode,  with  an  average  von  Mises  stress  of  9.22  KPa. 
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V.  Conclusions  and  Recommendations 


The  focus  of  this  research  was  to  conduct  a  system  identification  of  a  Manduca  Sexta 
forewing,  a  nd  develop  a  finite  el  ement  model  capable  of  accurately  p  redicting  the  o  bserved 
modes.  D  uring  this  process,  parametric  studies  were  conducted  examining  the  effects  of  vein 
diameter  and  camber  on  the  natural  frequencies  of  vibration.  C  onvergence  studies  and  a  stress 
analysis  w  as  c  onducted  by  nor  malizing  t  he  di  splacement  of  t  he  finite  e  lement  model  tot  hat 
obtained  b  y  t  he  1  aser  vibrometer.  The  co  nclusions  r  cached  d  uring  each  p  hase  oft  hese 
investigations  are  detailed  below. 

5.1.  Analysis  of  Boundary  Conditions 

A  foam-lined  clamp  was  used  to  secure  the  Manduca  Sexta  wing  during  the  experimental 
portion  of  the  system  identification.  Before  frequency  extraction  was  conducted  on  the  Manduca 
Sexta  wing,  a  p  aper  b  earn  w  as  p  laced  i  nside  t  he  cl  amp  and  an  alyzed  in  va  cuum  ( <  4  T  orr) 
conditions.  The  ratio  of  modes,  02/wi  is  6.40  and  ©3/02  is  2.813.  The  analytical  solution  for  the 
modes  of  a  p  rismatic,  homogeneous,  i  sotropic  be  am  w  ith  a  clamped  bounda  ry  c  ondition 
(cantilever)  were  developed  in  Chapter  II,  and  deviated  from  the  experimental  values  by  2.08% 
for  03/02  and  0.45%  for  03/02.  Thus,  the  foam-lined  clamp  used  during  the  frequency  analysis 
of  the  hawkmoth  wing  may  be  considered  rigidly  clamped. 

5.2.  Wing  Preparation  and  CT 

Wing  s  election  an  d  preparation  w  ere  car  efully  co  ntrolled  t  0  ensure  r  epeatability  0  f 
results.  O  nly  healthy,  adult,  female  Manduca  Sextae  were  used.  Before  wing  liberation,  each 
hawkmoth  w  as  p  laced  i  nto  a  r  efrigerator  (  34°F)  a  nd  cooled  for  3  0  m  inutes  t  0  i  nduce 
hypothermia.  T  he  moths  were  checked  for  responsiveness  before  wing  liberation.  Liberation 
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was  accomplished  using  a  scalpel,  and  the  wings  were  earefuhy  eut  adjaeent  to  the  thorax.  Onee 
removed,  the  wings  were  plaeed  into  a  foam-lined  elamp  for  analysis. 

Wing  ge  ometry  a  nd  di  mensions  w  ere  m  easured  us  ing  a  e  omputed  t  omography  (CT) 
imager.  T  he  C  T  s  oftware  er  eated  a  d  igital,  t  hree-dimensional  m  odel  us  ing  a  s  eries  of  X-Ray 
images.  C  T  ima  ging  1  asted  a  pproximately  20  min  utes,  a  nd  th  e  mo  del  w  as  s  aved  for  la  ter 
analysis.  Vein  dimensions,  eamber,  and  membrane  thiekness  were  measured,  and  the  membrane 
was  f  ound  t  o  be  of  a  g  enerahy  e  onstant  0.12  mm  thiekness.  Signifieant  e  stimation  o  f  v  ein 
diameter  was  neeessary,  due  to  the  low  resolution  of  the  CT  images.  In  order  to  ereate  a  model 
representative  o  f  t  he  Manduca  Sexta  wing,  it  is  r  eeommended  th  at  a  n  ima  ging  te  ehnique  o  f 
higher  resolution  than  CT  be  used. 

5.3.  System  Identification  of  a  Manduca  Sexta  Forewing 

A  Manduca  Sexta  forewing  was  plaeed  inside  the  foam-lined  elamp  deseribed  above,  and 
seeured  to  a  flapping  motor.  The  motor  exeited  the  wing  only  in  the  direetion  normal  to  the  wing 
planform  and  at  low  voltage  (0.2  V)  to  produee  small  displaeements.  A  pseudo-random  input 
was  used  from  frequeneies  approaehing  zero  to  1  kHz.  The  low  voltage  in  the  motor  produeed  a 
maximum  displaeement  of  27.29  pm  at  the  tip  during  the  first  bending  mode,  and  a 
eorresponding  maximum  veloeity  of  282.7  pm/s;  s  imilarly,  th  e  tip  d  isplaeement  d  uring  th  e 
torsional  mode  was  3.32  pm,  and  the  maximum  tip  veloeity  was  215.0  pm/s.  Frequeney  values 
and  m  odeshapes  w  ere  r  eeorded  us  ing  a  P  olytee®  s  eanning  1  aser  vi  bro meter.  T  he  f  irst  t  hree 
natural  f  requeneies  of  t  he  w  ing  i  n  vaeuum  w  ere  86  H  z,  106  H  z,  a  nd  155  H  z,  w  ith  a  n  e  rror 
margin  of  ±  2  H  z  at  ah  modes;  the  eorresponding  modeshapes  were  first  bending,  torsion,  and 
saddle,  respeetively. 
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Tests  were  conducted  in  air  and  vacuum  to  gage  the  extent  of  aeroelastic  effects  on  the 
Manduca  Sexta  wing.  Due  to  drying  concerns,  a  wing  cannot  be  used  for  air,  vacuum,  and  CT 
experiments  within  the  3-hour  freshness  criterion  set  by  Norris  [32],  T  o  circumvent  this  aging 
problem,  a  fresh  wing  was  used  for  comparison  of  modes  in  air  and  vacuum.  Vibration  tests  in 
air  r  esulted  i  n  f  requency  values  of  58.75  H  z,  75,  H  z,  a  nd  95  H  z  for  t  he  fi  rst  t  hree  m  odes. 
Frequency  testing  in  vacuum  (<  4  T  orr)  resulted  in  the  following  first  three  modes:  80  H  z,  98 
Hz,  and  118  Hz.  The  conclusion  from  these  experiments  is  that  aeroelasticity  plays  some  role  in 
the  f  requency  r  esponse  of  a  Manduca  Sexta  wing;  c  ontributions  f  rom  a  ir  de  crease  first  t  hree 
modes  by  an  average  of  23.17%. 

5. 4.  Linearity  of  Wing  Structure 

The  spectral  coherence  from  the  Manduca  Sexta  system  identification  is  plotted  in  Figure 
20.  O  nly  the  first  three  modes  of  the  forewing  were  considered  during  this  investigation,  the 
lowest  of  which  oc  curred  at  86  H  z,  while  the  third  mode  appeared  at  155  Hz.  T  he  average 
coherence  within  the  range  of  30  H z  to  200  H z  was  0.983,  w  ith  a  standard  deviation  of  0.014. 
The  high  coherence  value  observed  in  the  Manduca  Sexta  forewing  indicates  that  the  behavior  of 
the  wing  structure  is  linear  -  that  is,  satisfies  homogeneity  and  additivity  -  within  the  regime 
tested.  This  important  result  allows  scaling  of  stress  analyses  using  normalization  factors,  and 
Xi,  corresponding  to  the  ratios  of  the  tip  displacements  observed  in  experiment  to  the  first  and 
second  modes,  respectively. 

It  should  be  noted  that  although  a  1  inear  r esponse  w as  observed  during  low  amplitude 
flapping,  on  e  c  annot  c  onclusively  c  laim  that  the  w  ing  b  ehaves  lin  early  for  a  11  f  lapping 
conditions.  W  ithout  a  dditional  e  xperimentation,  i  t  is  unknow  n  w  hether  t  he  Manduca  Sexta 
wingbeat  produces  nonlinear  structural  effects. 
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5.5.  Inadequacy  of  Flat,  Homogeneous  Wings 


A  s  ystem  id  entification  o  f  a  paper  w  ing  o  f  t  he  s  ame  p  lanform  as  a  Manduca  Sexta 
forewing  was  also  eondueted  using  a  laser  vibrometer.  The  first  three  modeshapes  of  the  aetual 
Manduca  Sexta  wing  were  first  bending,  torsion  and  saddle;  the  first  four  modeshapes  of  the 
paper  wing  were  first  bending,  torsion,  seeond  bending,  and  saddle  (with  negative  tip-defleetion). 
The  ratios  of  the  Manduca  Sexta  wing  modes  in  vacuum  (<  4  Torr)  were  02/wi  =  1.23  and  03/01 
=  1.80;  the  ratios  of  the  paper  wing  modes  (in  vacuum)  were  02/01  =  4.80,  and  03/01  =  8.10. 
The  i  nclusion  of  t  he  s  purious  s  econd-bending  mode  a  nd  t  he  c  onsiderable  di  fferences  i  n  t  he 
frequency  ratios  highlight  the  inability  for  a  flat,  homogeneous,  isotropic  wing  to  characterize  a 
hawkmoth  wing.  Based  on  these  findings,  it  is  strongly  recommended  that  future  researchers  not 
conduct  a  nalogous  Manduca  Sexta  (or  a  ny  Lepidoptera)  w  ing  s  tudies  us  ing  f  lat,  one  -piece 
wings. 

5. 6,  Modeling  with  Membrane  Elements 

Finite  element  studies  were  conducted  to  examine  what  type  of  element  best  describes  the 
membrane  t  issue  found  i  n  t  he  Manduca  Sexta  wing.  A  m  odel  consisting  of  t  hree  b  earns 
separated  b  y  t  wo  m  embrane  s  ections  w  as  u  sed.  The  el  ements  selected  were  the  eight-noded, 
quadratic,  quadrilateral  membrane  (M3D8),  and  the  four-noded  linear  quadrilateral  membrane 
(M3D4).  Since  membrane  elements  only  carry  in-plane  stresses,  kinking  was  observed  at  the 
membrane-beam  i  nteraction  f  or  a  11  m  esh  s  izes.  Convergence  s  tudies  w  ere  c  onducted  s  bowing 
that  even  with  nearly  one  million  membrane  elements,  the  discontinuity  of  slope  at  the  beams 
could  not  be  r  esolved.  T  he  c  onclusion  f  rom  t  hese  s  tudies  i  s  t  hat  t  he  m  embrane  t  issue  i  n  t  he 
Manduca  Sexta  contains  be  nding  a  nd  t  ransverse  s  hear  c  omponents  not  carried  b  y  m  embrane 
elements. 
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5. 7.  Modeling  with  Shell  Elements 


The  model  used  in  the  investigation  of  membrane  elements  was  modified  to  eonsist  of 
three  beams  separated  by  two  shell  seetions.  Based  on  the  doeumentation  supplied  with  Abaqus, 
the  four-noded,  doubi  y-eurved,  r  edueed-integration,  s  hells  (S4R)  of  t  he  K  oiter-Sanders 
formulation  were  seleeted.  The  results  of  even  a  simple  model  (5,814  degrees  of  freedom)  were 
remarkably  better  than  those  observed  using  membrane  elements.  The  modeshapes  eontained  no 
disoontinuities  of  slope  at  the  beams.  Asa  result  of  these  studies,  shell  elements  were  seleeted 
for  future  modeling  of  the  Manduca  Sexta  wing. 

5.8.  Effects  of  Camber  on  Natural  Erequency 

A  homogeneous,  isotropic  shell  model  of  the  same  planform  as  the  Manduca  Sexta  wing 
was  constructed  in  Abaqus.  A  parametric  study  of  the  effect  of  camber  on  natural  frequency  was 
conducted.  Camber  values  between  0  %c  and  1 3  %c  (15”  angle  between  the  leading  edge  and 
maximum  mean  camber  line  location)  were  considered,  using  a  converged  mesh  of  300  elements 
(954  nodes,  5724  degrees  of  freedom).  All  modes  except  the  torsional  mode  were  observed  to  be 
strongly  dependent  on  c  amber-induced  stiffness.  T  he  spurious  second-bending  mode  observed 
in  the  flat  paper  wing  disappeared  at  a  wing  camber  of  4.5  %c,  at  which  point  the  saddle  mode 
became  the  third  mode. 

5.9.  Effects  of  Vein  Diameter  on  Natural  Erequency 

A  parametric  study  of  the  effects  of  vein  thickness  on  natural  frequency  was  conducted 
using  a  flat  wing  of  the  same  planform  as  the  Manduca  Sexta  forewing.  Ten  veins  were  modeled 
as  p  rismatic  b  earns  u  sing  average  cr  oss-sectional  da  ta  obt  ained  from  C  T  i  mages.  A  v  ein 
thickness  factor  was  applied  as  a  s  caling  coefficient  for  all  veins  and  the  first  four  modes  were 
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recorded.  Scaling  factors  between  0.5  and  2  were  used,  and  data  were  collected  at  semi  regular 
intervals  with  a  thickness  factor  step  size  of  approximately  0.4.  At  a  vein  thickness  factor  of  1.2 
{i.e. ,  t  he  ve  ins  cross  s  ections  w  ere  20%  1  arger  t  han  obs  erved  f  rom  C  T  images),  t  he  s  purious 
second-bending  mode  disappeared,  and  the  saddle  mode  became  the  third  mode. 

All  ve  ins  w  ere  s  trongly  c  orrelated  t  o  ve  in  t  hickness.  The  f  irst  be  nding  m  ode  a  nd 
torsional  mode  increased  at  approximately  the  same  rate  with  vein  thickness,  though  the  saddle 
mode  di  verged  f  rom  t  he  f  irst  be  nding  a  nd  t  orsion  a  s  v  ein  d  iameter  i  ncreased.  T  he  s  econd 
bending  m  ode  c  hanged  b  y  t  he  greatest  f  requency  v  alue  (  129  H  z)  i  n  t  he  dom  ain  of  t  he 
investigation,  while  the  first  bending  mode  increased  by  the  greatest  factor  (4.42).  There  are  two 
conclusions  from  this  study:  first,  adequately  modeling  the  natural  frequencies  of  a  Manduca 
Sexta  wing  is  strongly  d ependant  on  obtaining  correct  vein  dimensions;  second,  the  frequency 
difference  between  the  first  bending  and  torsional  modes  is  relatively  fixed,  though  the  frequency 
difference  between  the  first  and  third  modes  diverges  rapidly  with  increasing  vein  diameter. 

5.10,  Material  Property  Considerations 

The  ma  terial  p  roperties  o  f  th  e  Manduca  Sexta  are  unknow  n  a  s  of  t  his  doc  ument’s 
publication.  Published  values  for  the  Tibicen  Canicularis  (cicada)  [25]  and  preliminary  findings 
by  O’Hara  [31]  place  the  material  properties  of  the  vein  and  membrane  structures  between  1.9-5 
GPa.  Mengesha,  V allance,  Barraja  and  Mittal  report  the  veins  and  membrane  having  Y oung’s 
moduli  of  1.90  G  Pa  and  3.70  G  Pa,  r  espectively  [25].  O  ’Hara’s  r  esearch  i  ndicates  t  hat  t  he 
hawkmoth’s  material  values  may  be  the  reverse  of  the  cicada’s,  with  a  membrane  modulus  of  2 
GPa  [31].  Although  the  values  reported  in  Refs.  [25]  and  [31]  vary  by  a  factor  of  approximately 
two,  they  are  within  the  same  order  of  magnitude. 
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The  vein  diameter  studies  of  Chapter  III  indieate  that  varying  Y oung’s  modulus  within 
the  regime  of  2-5  GPa  plays  a  relatively  small  role  in  ehanging  the  natural  frequencies  of  the 
Manduca  Sexta  wing.  Examining  the  vein  diameter  plots  at  a  vein  thickness  factor  of  one  gives 
the  following  frequency  values  using  the  material  properties  from  the  cicada  wing:  ©i  =  33  Hz, 
©2  =  89  Hz,  ©3  =  206  Hz,  and  ©4  =  297  Hz.  In  comparison,  reversing  the  moduli  of  the  veins  and 
membranes  gives  the  following  frequency  values:  ©1  =  41  Hz,  ©2  =  91  Hz,  ©3  =  226  Hz,  and  ©4 
=  308  Hz.  While  these  changes  are  significant  in  the  first  mode  (a  change  of  19.5%),  modulus 
variation  n  evertheless  h  as  a  1  ess  dramatic  effect  than  al  tering  cam  ber  or  v  ein  d  iameter  b  y  a 
similar  magnitude  (a  factor  of  two,  in  this  case). 

Density  values  for  the  Manduca  Sexta  are  also  unknown.  Unlike  the  elastic  moduli,  the 
mass  densities  of  the  veins  and  membrane  have  yet  to  be  evaluated  -  even  preliminarily.  The 
density  V alues  used  during  the  modeling  process  were  taken  from  the  cicada  (2.30  g/cm  )  and 
produced  a  model  with  mass  representative  of  a  Manduca  Sexta  (0.04  g).  During  this  analysis, 
an  inherent  assumption  is  that  the  densities  of  the  membrane  and  veins  are  equal;  this  assumption 
may  be  t  he  r  oot  c  ause  of  r  dative  i  nvariance  off  requency  w  ith  c  hanging  Y  oung’s  m  oduli. 
Borrowing  f  rom  b  earn  t  heory,  t  he  ratio  of  e  lastic  m  odulus  1 0  m  ass  de  nsity,  E/ju,  is  t  he  onl  y 
influence  of  material  property  values  on  natural  frequency.  The  others  -  length  and  moment  of 
inertia  -  are  set  by  geometry.  In  the  case  of  the  Manduca  Sexta,  it  may  be  that  the  ratio  of  E/ju 
plays  a  critical  role  that  is  understated  by  the  equal  vein  and  membrane  density  assumption.  In 
short,  before  a  more  accurate  model  can  be  produced.  Young’s  modulus  and  density  values  must 
be  i  nvestigated  f  urther.  W  ithout  doi  ng  s  0,  it  i  s  i  mpossible  to  ve  rify  t  he  c  onclusions  of  the 
material  pr  operty  i  nfiuence  on  f  requency  given  i  n  Chapter  I II,  b  ecause  t  hey  a  re  ba  sed  upon 
unvalidated  underlying  assumptions  regarding  density. 
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5.11.  Modeling  and  Model  Limitations 


A  finite  element  model  of  a  Manduca  Sexta  wing  was  produeed  to  accurately  predict  the 
first  tw  0  natural  f  requencies  an  d  m  odeshapes  as  seen  in  th  e  e  xperimental  d  ata.  T  he  th  ird 
modeshape  is  q  ualitatively  s  imilar  to  th  at  o  bserved  in  th  e  Manduca  Sexta  wing,  t  hough  t  he 
frequency  v  alues  di  ffer  b  y  102%  .  T  he  1  aser  vi  brometry  experiments  pr  edicted  m  odes  i  n  t  he 
Manduca  Sexta  wing  at  86  Hz,  106  Hz,  and  155  Hz,  with  an  error  margin  of  ±  2  Hz  for  each 
mode.  T  he  finite  element  model  gives  the  first  bending  mode  at  84.6  H  z,  a  torsional  mode  at 
106.1  Hz,  and  a  saddle  mode  at  317.7  Hz.  These  values  were  obtained  using  a  converged  model 
with  a  total  of  80,633  quadratic  beam  and  general-purpose  shell  elements,  comprised  of  83,664 
nodes.  T  he  agreement  of  the  first  two  natural  frequencies  of  the  experimental  wing  and  finite 
element  model  suggest  that  the  model  is  useful  for  characterizing  frequencies  up  to  and  including 
the  s  econd  m  ode  onl  y.  Further  r  esearch  i  s  r  equired  t  o  de  termine  t  he  i  mportance  of  t  he 
experimental  and  numerical  disagreement  seen  in  the  third  and  higher  natural  frequencies. 

Comparing  t  he  m  odd  t  o  t  he  pa  rametric  s  tudies  of  c  amber  a  nd  ve  in  d  iameter  v  ersus 
natural  frequency  suggest  more  p  recise  mo  deling  is  required.  W  hile  it  is  p  ossible  to  mo  del  a 
more  accurate  distribution  of  camber  using  the  CT  images,  the  resolution  is  not  sufficient  to 
accurately  m  easure  v  ein  p  rofiles.  I  t  i  s  s  trongly  r  ecommended  t  hat  hi  gher  r  esolution  i  maging 
techniques  be  investigated  if  accuracy  beyond  the  second  mode  is  required. 

5.12.  Stress  Analysis  of  the  First  Bending  Mode 

A  stress  analysis  was  conducted  on  the  first  bending  mode  (84.6  Hz)  of  the  finite  element 
model.  A  s  caling  factor,  Ai,  w  as  de  veloped  ton  ormalize  t he  di  splacement  r  esults  r  eported  i n 
Abaqus  to  those  of  the  wing  examined  using  a  laser  v ibrometer.  T  he  scaling  factor,  li,  was 
defined  as  t  he  r  atio  o  f  the  t  ip  di  splacement  r  ecorded  e  xperimentally  t  o  t  he  tip  d  isplacement 


102 


reported  by  Abaqus:  Ai  =  5psv/5Abq,  and  was  found  to  be  8.92  xlO'  for  the  first  bending  mode. 
Sinee  the  Manduca  Sexta  wing  is  behaves  linearly,  the  sealing  faetor  approaeh  is  valid. 

The  maximum  stress  in  the  first  bending  mode  oeeurs  in  the  trailing  edge  at  the  elamped 
boundary.  The  maximum  (l-sealed  output)  von  Mises  stress  experieneed  in  the  first  mode  was 
258.4  KPa.  The  minimum  von  Mises  stress  was  eomputed  to  be  567  Pa,  and  oeeurred  at  the 
wing  tip.  N  otable  was  the  eonsistently  low  stress  observed  in  the  membrane,  while  the  veins 
earned  the  bulk  of  the  stresses  in  this  mode.  V  ein  (beam)  stress  was  observed  to  deerease  with 
span  in  a  11  me  mbers,  with  s  tress  e  oneentrations  a  rising  a  t  th  e  v  ein  in  terseetions.  The  e  yelie 
operation  a  nd  hi  gh  s  tress  e  oneentrations  ne  ar  t  he  w  ing  r  oot  s  uggest  f  atigue  f  ailure  m  ay  be  a 
design  ehallenge  in  future  MAVs;  appropriate  materials  and  geometry  must  be  earefully  seleeted 
to  seeure  the  wing-body  interfaee. 

5.13,  Stress  Analysis  of  the  Torsional  Mode 

Using  the  approaeh  deseribed  above  to  determine  the  sealing  faetor,  X,  a  stress  analysis  of 
the  torsional  mode  was  eondueted.  The  sealing  faetor  eorresponding  to  the  torsional  mode,  X2, 
was  e  omputed  1 0  be  1 .  68  x  10'  .  The  m  aximum  (  sealed)  v  on  M  ises  s  tress  i  n  t  he  t  orsional 
eigenveetor  was  28.91  KPa,  loeated  in  the  leading  edge  at  approximately  80%  of  the  span.  The 
minimum  von  Mises  stress  was  277  Pa,  loeated  in  the  membrane  just  aft  of  the  leading  edge, 
near  the  14  span  loeation.  The  stress  variation  in  the  torsional  mode  was  signifieantly  less  than 
that  of  the  first  bending  mode.  While  the  veins  still  earried  the  majority  of  the  load,  a  far  greater 
stress  di stribution  w as  found  in  the  membrane.  S  tress  was  observed  to  be  a  funetion  of  span, 
with  eoneentrations  in  the  ehordwise  vein  not  observed  in  the  previous  mode. 
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5.14.  Recommendations  for  Future  Work 


With  the  above  eonelusions  in  mind,  it  is  recommended  that  future  researchers  work  to 
reduce  the  number  of  model  approximations  and  conduct  further  parametric  analyses.  From  an 
experimental  pe  rspective,  f  uture  w  ork  s  hould  f  ocus  on  obt  aining  high-fidelity  g  eometry 
measurements.  S  econd,  the  material  properties  of  the  Manduca  Sexta  wing  components  should 
be  investigated.  Third,  in  light  of  the  observed  frequency  shift  when  comparing  air  and  vacuum 
measurements,  one  s  hould  unde  rtake  a  c  omprehensive,  noni  nvasive  e  xperimental  s  tudy  o  f  t  he 
aeroelastic  effects  on  the  Manduca  Sexta  wing. 

Focusing  on  possible  numerical  research,  a  model  that  represents  camber  as  a  function  of 
wing  coordinates  should  be  created,  with  special  attention  paid  to  the  bi-camber  observed  using 
the  CT  imager.  F  inally,  to  aid  future  MAV  designers,  a  parametric  study  of  the  wing  structure 
and  its  components  should  be  conducted.  At  this  point,  it  is  evident  that  the  venation  observed  in 
the  hawkmoth  wing  has  a  dramatic  effect  on  frequency  response;  it  is  unclear,  however,  whether 
the  specific  venation  pattern  is  optimized  for  flight  or  circulation  -  or  neither. 
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Appendix  A  -  CoordGen 


The  following  software  uses  input  files  generated  in  Abaqus  to  write  text  files  eontaining 
the  t  hree-dimensional  e  oordinates  o  f  geometrie  s  ets.  T  hese  s  ets  ar  e  ereated  b  y  t  he  u  ser  i  n 
Abaqus  during  the  modeling  proeess.  The  purpose  of  CoordGen  is  to  address  Abaqus’s  inability 
to  internally  sketch  wire  beams  in  three  dimensions.  CoordGen  is  a  program  written  in  Python 
2.6,  and  feeds  coordinates  for  splines  into  SolidWorks;  these  splines  can  be  imported  back  into 
Abaqus  as  wire  beams  in  a  common  format  (IGES,  STEP,  etc.).  With  permission,  the  code  is 
reprinted  below. 


#  Copyright  (c)  2010  Jason  R.  Sims,  jsims02@gmail.com 

#  Permission  is  hereby  granted,  free  of  charge,  to  any  person 

#  obtaining  a  copy  of  this  software  and  associated  documentation 

#  files  (the  "Software"),  to  deal  in  the  Software  without 

#  restriction,  including  without  limitation  the  rights  to  use, 

#  copy,  modify,  merge,  publish,  distribute,  sublicense,  and/or  sell 

#  copies  of  the  Software,  and  to  permit  persons  to  whom  the 

#  Software  is  furnished  to  do  so,  subject  to  the  following 

#  conditions: 

#  The  above  copyright  notice  and  this  permission  notice  shall  be 

#  included  in  all  copies  or  substantial  portions  of  the  Software. 

#  THE  SOFTWARE  IS  PROVIDED  "AS  IS",  WITHOUT  WARRANTY  OF  ANY  KIND, 

#  EXPRESS  OR  IMPLIED,  INCLUDING  BUT  NOT  LIMITED  TO  THE  WARRANTIES 

#  OF  MERCHANTABILITY,  FITNESS  FOR  A  PARTICULAR  PURPOSE  AND 

#  NONINFRINGEMENT.  IN  NO  EVENT  SHALL  THE  AUTHORS  OR  COPYRIGHT 

#  HOLDERS  BE  LIABLE  FOR  ANY  CLAIM,  DAMAGES  OR  OTHER  LIABILITY, 

#  WHETHER  IN  AN  ACTION  OF  CONTRACT,  TORT  OR  OTHERWISE,  ARISING 

#  FROM,  OUT  OF  OR  IN  CONNECTION  WITH  THE  SOFTWARE  OR  THE  USE  OR 

#  OTHER  DEALINGS  IN  THE  SOFTWARE. 


#  Following  are  objects  and  their  properties: 

#  Instance 

#  - 

#  *  name  -  string 

#  *  nodes  -  list  of  Nodes 

#  *  nsets  -  list  of  Nsets 

# 
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Node 


# 

# 

# 

# 

# 

# 

# 

# 

# 

# 

# 

# 


*  X  -  x-coord  (float) 

*  y  -  float 

*  z  -  float 

*  id  -  id  in  Node  list  (int) 

Nset 

*  name  -  string 

*  instance  name  -  string 

*  nodes  -  list  of  ints  representing  Node  ids 


from  optparse  import  OptionParser 


class  Context (object) : 
is  done  =  False 

def  add  line (self,  line) : 

if  self.  class  .end  pred(line) : 

self. is  done  =  True 
else : 

self,  add  line (line) 

class  Instance (Context ) : 

def  init  (self,  line) : 

Context.  init  (self) 
self. nodes  =  list() 
self.nsets  =  list() 

self,  parsing  nodes  =  False 

self,  parse  name (line) 

def  add  line (self,  line) : 
if  self,  parsing  nodes: 

if  Node. end  pred(line) : 

self,  parsing  nodes  =  False 
else : 

self . nodes . append (Node (line)  ) 
elif  Node. start  pred(line) : 

self,  parsing  nodes  =  True 

def  parse  name (self,  line) : 

for  res  in  filter ( lambda  1:  "name"  in  1,  line . split (",") ) 
self. name  =  res.strip()  .split("  =  ")  [1] 

@classmethod 

def  start_pred (els ,  line): 

return  "*Instance"  in  line 
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@classmethod 

def  end  pred(self,  line) : 

return  "*End  Instance"  in  line 

class  Node: 

def  init  (self,  line) : 

self. id,  self.x,  self.y,  self.z  =  [float(s)  for  s  in 
line . split  ( " , " ) ] 

@classmethod 

def  start_pred (els ,  line): 
return  "*Node"  in  line 

@classmethod 

def  end  pred(self,  line) : 
return  in  line 

class  Nset (Context ) : 

def  init  (self,  line) : 

Context.  init  (self) 
self. nodes  =  list() 
self,  parse  name (line) 

def  add  line (self,  line) : 

self. nodes  +=  [int(n)  for  n  in  line . split (",") ] 

def  parse  name (self,  line) : 

for  res  in  filter ( lambda  1:  "nset"  in  1,  line . split (",")) : 

self. name  =  res.strip()  .split("  =  ")  [1] 
for  res  in  filter ( lambda  1:  "instance"  in  1,  line . split (",")) : 
self . instance  name  =  res. strip  ()  .split ("  =  ")  [1] 

@classmethod 

def  start_pred (els ,  line): 
return  "*Nset"  in  line 

@classmethod 

def  end  pred(self,  line) : 
return  in  line 


def  map_nsets_to_instances (elts)  : 

def  partition ( (insts ,  nsets)  ,  elt)  : 
if  type (elt)  is  Instance: 

insts . append (elt ) 
else : 

nsets . append (elt ) 
return  insts,  nsets 

insts,  nsets  =  reduce (partition,  elts,  [[],  []]) 

for  nset  in  nsets: 

for  inst  in  insts: 
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if  nset . instance  name  ==  inst.name: 
inst.nsets. append (nset ) 

return  insts 

def  gen  lines (file) : 

fd  =  open (file,  "r") 
for  line  in  fd: 

yield  line 
f d . close  ( ) 

def  process_file (f ile) : 

def  do_process ( (context,  acc) ,  line) : 
if  context  is  not  None: 

context. add  line  (line) 
if  context. is  done: 

return  (None,  acc) 
return  (context,  acc) 
for  kls  in  [Instance,  Nset] : 
if  kls. start  pred(line): 
context  =  kls (line) 
acc. append (context) 
return  (context,  acc) 

return  map_nsets_to_instances (reduce (do_process,  gen_lines (file) , 

(None, 

list  0 ))[!]) 


def  write  instance (inst)  : 
for  nset  in  inst.nsets: 

fd  =  open ( "%s . txt"  %  nset. name,  "w") 
for  n  in  nset. nodes: 

node  =  inst . nodes [n-1 ] 

fd. write ("%s  %s  %s\n"  %  (node.x,  node.y,  node.z)) 
f d . close ( ) 

if  name  ==  "  main  ": 

parser  =  OptionParser (usage="coordgen . py  [files]") 
(options,  args)  =  parser . parse_args ( ) 

map (write_instance ,  sum (map (process_file,  args),  [])) 
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Appendix  B  -  CT  Images 


These  i  mages  are  s  upplemental  tot  hose  s  upplied  i  n  Chapter  1 1.  Figure  73  shows  t  he 
unaltered  version  of  Figure  19.  Although  visible  at  high-zoom  levels,  the  inner  portions  of  the 
leading  and  trailing  edge  veins  are  diffieult  to  distinguish.  U nmagnified  versions  of  this  image 
do  not  adequately  reproduce  the  grey  discoloration  present  to  the  CT  operator.  Asa  result,  the 
areas  inside  the  inner  diameters  were  converted  to  white,  with  a  3%  coloration  threshold.  In 
other  words,  if  the  color  of  the  pixels  adjacent  to  the  inner  portion  was  more  than  3%  different 
than  the  grey  inside  the  vein,  it  kept  its  original  coloring.  By  applying  this  color  change  with  a 
small  threshold,  the  reader  is  better  able  to  see  the  circulatory,  inside  portion  of  the  vein.  Finally, 
the  noise  at  the  left  of  the  image  was  removed  in  Figure  19. 


Figure  73:  Unaltered  version  of  Figure  19. 

The  next  figure  shows  the  exact  locations  of  the  cross-sectional  measurements  described 
in  Chapter  1 1.  The  s  olid,  v  ertical  1  ines  r  epresent  t  he  digital  “slices”  taken  t  o  eh  aracterize  t  he 
geometry  of  the  wing.  Finally,  Figure  75  shows  the  stereo-lithographic  planform  view  of  the 
Manduca  Sexta  forewing,  as  exported  from  the  CT. 
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Figure  74:  Location  of  cross-sectional  measurements,  denoted  by  vertical  lines. 


Figure  75:  Stereo-lithographic  planform  view  of  the  Manduca  Sexta  wing. 
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Appendix  C  -  Higher  Order  Modes 


The  following  images  show  the  modeshapes  of  the  hawkmoth’s  fourth  and  fifth  modes. 
These  modes  were  collected  in  the  laboratory,  but  omitted  from  analysis.  Since  every  continuous 
structure  has  an  infinite  number  of  natural  frequencies  and  associated  mode  shapes,  one  is  forced 
to  narrow  an  analysis  to  a  set  number  of  modes.  In  this  case,  the  first  three  modes  were  deemed 
sufficient  for  a  preliminary  finite  element  comparison. 


Figure  76:  Fourth  modeshape  (187  Hz)  of  a  Manduca  Sexta  wing, 
shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 


Figure  77:  Fifth  modeshape  (23 1  Hz)  of  a  Manduca  Sexta  wing, 
shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 


The  following  figure  shows  the  fifth  modeshape  of  the  paper  wing  tested  for  comparison 
to  th  e  h  awk  mo  th.  T  his  mo  de  h as  th  e  s  ame  qualitative  c  haracteristics  a s  th e  Manduca  Sexta 


ill 


saddle  mode,  with  the  exception  of  the  circular,  plate-like  vibration  in  the  center  (see  image  (a) 
in  Figure  78),  below. 


(b) 


Figure  78:  Fifth  modeshape  (692  Hz)  of  a  paper  Manduca  Sexta  wing, 
shown  from  an  isometric  view  (a)  and  a  tip  view  (b). 
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Appendix  D  -  Wing  Camber 


Wing  camber  is  defined  as  the  maximum  distance  between  the  mean  camber  line  and  the 
chord  line  [41],  T  ypically,  the  camber  is  expressed  as  a  p  ercentage  of  chord  (%c)  -  the  linear 
distance  b  etween  t  he  1  eading  an  d  t  railing  ed  ges.  T  he  m  ean  cam  ber  1  ine  ( technically  a  cu  rve, 
though  c  ailed  a  “  line”  b  y  convention)  i  s  t  he  c  ollection  of  poi  nts  t  hat  i  s  ha  If- way  b  etween  t  he 
upper  and  lower  surfaces  of  the  wing  -  hence  mean  camber  line).  Figure  79  shows  a  side  view 
of  a  conventional  airfoil,  with  the  mean  camber  line,  chord  line,  thickness  and  camber  labeled. 

Mean  Camber 


Figure  79:  Wing  geometry  image,  eourtesy  NASA  Glenn  Research  Center. 

The  Manduca  Sexta  wing  has  a  thin  membrane  structure,  such  that  the  mean  camber  line 
is  essentially  the  midpoint  of  the  membrane  curvature.  As  discussed  in  Chapter  IV,  there  is  a 
difference  b  etween  t  he  camber  m  odeled  i  n  t  he  finite  e  lement  studies,  a  nd  t  he  c  amber  of  t  he 
actual  Manduca  Sexta  wing.  The  finite  element  model  features  a  single  shell  of  constant  radius, 
to  w  hich  b  earns  a  re  tie  d  a  t  th  e  mid  point,  ass  hown  i  n  Figure  80.  T  he  a  ctual  w  ing,  how  ever, 
features  a  double  camber,  with  both  positive  and  negative  concavity,  as  shown  in  Figure  81. 
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Mean  camber  line 


Chord,  c 


Figure  80:  Schematic  of  modeled  FE  wing  at  the  root. 


Figure  8 1 :  Schematic  of  a  Manduca  Sexta  wing  at  the  root. 
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Appendix  E  -  Stress  Plots 


The  images  that  follow  are  contour  plots  corresponding  to  the  stress  analyses  of  Chapter 
IV.  Figure  82  shows  the  von  M  ises  stress  contour  of  the  first  bending  mode.  Figure  83  also 
shows  the  stress  in  the  first  bending  mode,  with  the  shells  made  transparent  for  easier  viewing  of 
the  beams. 


Min 


Max 

Figure  82:  First  mode  von  Mises  stress  contour  for  a  finite  element  model  of  a  Manduca  Sexta  forewing. 


Figure  83:  First  mode  stress  in  a  finite  element  model  showing  Manduca  Sexta  wing  veins. 
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Figure  84  and  Figure  85  show  the  von  M  ises  stress  contour  for  the  second  mode  (torsion).  As 
before,  t  he  s  hells  ar  e  p  resent  i  n  Figure  85,  b  ut  a  re  tr  ansparent  to  b  etter  illu  strate  th  e  s  tress 
distribution  in  the  veins. 


Figure  84:  Second  mode  stress  in  a  finite  element  model  of  a  Manduca  Sexta  forewing. 


Figure  85:  Second  mode  stress  in  a  finite  element  model  showing  Manduca  Sexta  wing  veins. 
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